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Abstract 

As  legacy  neutron  irradiation  facilities  are  shut  down  due  to  security  and  financial 
constrictions,  a  growing  need  arises  for  alternatives  that  can  provide  the  same  or 
similar  radiation  environments  using  methods  and/or  fuels  that  do  not  pose  the  same 
risks.  For  this  reason,  facilities  that  provide  much  lower  neutron  fluxes  and  different 
spectral  shapes  are  being  leveraged  over  longer  irradiation  intervals  to  meet  this  need. 
However,  the  question  arises  as  to  whether  the  use  of  these  types  of  facilities  provide 
a  valid  comparison  to  the  legacy  systems’  results. 

To  this  end,  a  model  using  a  system  of  coupled,  non-linear  ordinary  differential 
equations  has  been  developed  to  track  defect  species  in  silicon  for  short  pulse  neutron 
irradiations.  This  model  has  been  used  to  predict  current  gain  degradation  in  silicon 
BJTs  for  various  neutron  fluences.  These  predictions  have  been  compared  against 
experimental  data  collected  at  two  neutron  irradiation  facilities  with  different  time 
and  spectral  prohles.  The  damage  constant  during  irradiation  has  been  determined, 
and  it  is  different  for  both  facilities.  However,  the  time  prohle  is  found  to  have  no 
effect  in  the  region  tested  in  this  work.  Now  that  this  analysis  has  been  done,  these 
types  of  facilities  can  be  used  for  radiation  vulnerability  analysis  studies  for  use  in 
short  pulse  neutron  damage  studies. 
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INVESTIGATING  TIME  AND  SPEGTRAL  DEPENDENGE  IN  NEUTRON 
RADIATION  ENVIRONMENTS  FOR  SEMIGONDUGTOR  DAMAGE  STUDIES 

I.  Introduction 


1.1  Overview 

In  the  radiation  effects  on  electronics  community,  the  overarching  goal  of  every 
theoretical  analysis,  simulation,  and  experiment  is  to  understand  the  response  of 
electronic  devices  and  components  when  subjected  to  a  type  and  level  of  radiation. 
Neutron  radiation  sources  such  as  fast  burst  reactors  at  the  Sandia  Pulsed  Reactor 
(SPR-III)  [1],  the  Aberdeen  Proving  Ground  (APG)  [2],  and  White  Sands  Missile 
Range  (WSMR)  [3]  have  been  the  gold  standard  for  pulsed  neutron  radiation  damage 
studies.  The  harsh,  but  well-quantihed,  neutron  environments  at  these  facilities  were 
used  to  gain  government  and  industry  approval  for  components  used  in  various  sys¬ 
tems,  such  as  missiles  and  satellites.  However,  some  of  these  facilities  have  been  shut 
down  due  to  a  growing  interest  in  the  security  of  Special  Nuclear  Material  (SNM). 

While  this  happens,  researchers  search  for  alternatives  that  can  provide  the  same 
or  similar  radiation  environments  using  methods  and/or  fuels  that  do  not  pose  a 
safety  or  security  risk.  Examples  of  these  alternatives  include  thermal  neutron  reac¬ 
tors  such  as  the  Annular  Gore  Research  Reactor  (AGRR),  neutron  surrogate  facilities, 
and  accelerator-based  neutron  sources.  Each  of  these  systems  has  benehts  and  disad¬ 
vantages,  but  there  is  one  trait  that  they  all  share.  This  trait  is  a  lower  time-averaged 
neutron  flux  requiring  longer  irradiation  times  to  reach  the  total  fluence  required  by 
scientists  making  comparisons  to  legacy  data  collected  from  fast  burst  reactors. 
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Conventional  wisdom  has  held  that  a  spectrum  of  neutrons  can  be  collapsed  into 
an  equivalent  number  of  mono-energetic  neutrons.  This  is  done  using  the  concept  of 
damage  Kinetic  Energy  Released  per  unit  MAss  (kerma),  or  the  amount  of  kinetic 
energy  deposited  in  matter  [4],  which  can  then  be  used  to  calculate  the  expected 
damage  in  Si  devices  for  a  given  fluence.  In  the  past,  this  has  been  verihed  for 
fluences  up  to  and  exceeding  10^^  n/cm^  at  facilities  such  as  the  SPR-III.  However,  as 
these  fast  burst  reactors  are  decommissioned  or  no  longer  widely  used,  the  question 
arises  as  to  whether  the  use  of  alternative  facilities  with  longer  irradiation  times  is  a 
valid  comparison  to  the  legacy  fast  burst  systems. 

It  is  the  overarching  goal  of  this  research  to  investigate  the  effect  of  temporal  and 
spectral  variance  on  Bipolar  Junction  Transistor  (BJT)  performance  for  facilities  that 
require  a  signihcant  amount  of  irradiation  time  to  reach  desired  neutron  fluences,  both 
constant  and  continuously-pulsed  sources.  In  this  work,  only  Si  n^pn-tjpe  BJTs  are 
considered,  since  these  devices  have  a  signihcant  data  history  and  are  well-understood. 

Legacy  Test  Systems. 

While  short  pulse  neutron  sources  have  been  used  for  decades  and  numerous  legacy 
data  sets  exist  in  databases  around  the  world,  a  growing  interest  in  the  control  and 
security  of  SNM  has  resulted  in  shutting  down  a  number  of  these  crucial  facilities.  The 
SPR-III  was  decommissioned  in  2007  amid  a  hurry  of  last-minute  experiments  [5,6]. 
The  U.S.  Army’s  fast  burst  reactor  at  APG  was  long-ago  decommissioned  due  to 
budget  and  security  concerns.  The  United  States’  only  remaining  fast  burst  reactor 
is  situated  at  the  WSMR  in  Arizona.  It  continues  to  operate,  however  it  also  faces 
budgetary  and  security  concerns.  It  is  for  this  reason  that  the  community  continues 
to  look  for  alternatives  to  facilities  fueled  by  SNM. 

The  American  Society  for  Testing  and  Materials  (ASTM)  provides  a  standard 
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method  to  test  BJT  response  to  neutrons  [7].  In  it,  the  prescribed  requirements 
include  very  few  hard-and-fast  stipulations.  However,  it  suggests  that  the  source 
provide  total  fluence  of  particles  in  very  short  bursts,  often  in  100  ps  or  less.  With 
fast  burst  reactors,  it  was  an  inherent  physical  safeguard  that  the  reaction  hnished  in 
this  time  frame  or  less,  which  was  one  reason  why  the  standard  was  written  this  way. 
However,  with  alternative  methods  currently  proposed,  it  can  be  difficult  to  reach 
relevant  neutron  fluence  levels  in  such  a  short  time  interval. 

Neutron  Source  Alternatives. 

Alternative  neutron  sources  for  radiation  testing  are  many  and  varied.  Fission 
reactors  are  sources  that  are  typically  driven  by  traditional  nuclear  reactors  fueled 
with  low-  to  medium-enrichment  uranium  or  mixed-oxide  fuel.  Examples  of  this 
source  include  the  Annular  Core  Research  Reactor  (ACRR)  [8]  as  well  as  university 
research  reactors.  Often,  they  can  be  operated  in  both  pulsed  or  steady-state  mode. 
The  pulse  is  caused  by  a  large  reactivity  insertion  instigated  by  quickly  withdrawing 
a  reactor  control  rod  from  the  core,  typically  using  pneumatic  or  hydraulic  methods. 
The  excess  reactivity  causes  the  reactor  to  go  prompt  supercritical  for  a  short  time 
before  the  control  rod  falls  back  into  the  reactor,  stifling  the  reaction. 

The  advantage  of  this  alternative  is  that  the  fuel,  while  still  nuclear  material,  is  of 
much  lower  enrichment  than  for  bare  metal  reactors  such  as  the  SPR-HI.  Additionally, 
while  the  pulse  width  of  the  critical  transient  is  short,  neutron  fluences  are  very  high. 
Unfortunately,  since  the  coolant  and  moderator  is  usually  light  water,  the  neutron 
spectrum  is  heavily  shifted  to  the  thermal  region.  Because  of  the  shift  to  the  thermal 
region,  facilities  such  as  the  ACRR  are  typically  used  for  components  that  will  be 
primarily  used  in  terrestrial  applications,  since  most  damaging  neutron  irradiation  in 
the  oxygen-rich  environment  near  the  earth’s  surface  has  been  greatly  thermalized 
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and  the  pulse  subsequently  stretched  out  in  time  before  reaching  the  component  in 
question.  However,  for  space-based  systems,  this  is  not  the  case,  as  much  of  the 
neutron  radiation  (either  cosmic  or  nuclear  weapon  induced)  remains  hard  and  the 
radiation  pulse  short.  Also,  the  gamma  flux  from  hssion  and  fission  products  in 
thermal  reactor  systems  is  relatively  high,  which  adds  additional  complexity  to  the 
analysis  of  experimental  results. 

Charged  particle  accelerators  bombard  a  Device  under  Test  (DuT)  with  ions  and 
equate  the  resultant  damage  to  a  fluence  of  neutrons  via  a  transfer  function  [9].  An 
example  of  this  system  is  the  Ion  Beam  Lab  (IBL)  at  Sandia  National  Laboratory 
(SNL).  Since  the  ions  are  charged,  packets  of  them  can  be  accelerated  through  an 
electromagnetic  held  and  directed  to  a  DuT.  An  advantage  to  this  method  is  the 
ability  to  control,  manipulate,  and  focus  the  incident  particles  to  the  target,  which 
can  lead  to  high  huence  and  short  pulse  widths.  Also,  since  the  charged  particles 
deposit  most  of  their  energy  at  the  end  of  their  track,  the  device  being  tested  can 
be  probed  for  sensitivity  in  different  areas.  However,  a  disadvantage  is  that,  because 
the  experimenter  is  not  using  neutrons,  an  approximation  must  be  made  for  damage 
equivalence.  Additionally,  since  charged  particles  deposit  almost  all  of  their  energy  at 
the  end  of  their  track  in  a  Bragg  peak,  the  damage  that  is  caused  is  highly-localized, 
which  is  contrary  to  the  homogeneity  of  damage  caused  by  neutrons.  This  may  be 
helpful,  if  the  experimenter  wishes  to  probe  a  very  specihc  region  of  a  device;  however, 
it  does  not  simulate  the  damage  done  by  neutrons  exactly  and  must  be  taken  into 
account. 

Compact  accelerator-driven  neutron  sources  encompass  a  wide  range  of  irradiation 
facilities.  Neutron  generators  use  very  small  linear  accelerators  to  fuse  deuterium, 
tritium,  or  a  combination  of  the  two  [10].  Each  combination  provides  a  different 
total  flux  and  spectrum  of  neutron  energies.  An  idealized  plot  of  these  spectra  are 
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plotted  in  Figure  1.  Typically,  these  sources  run  at  very  low  flux  and  are  used  for 
applications  more  pertinent  to  medicine  or  more  delicate  experiments.  More  recently, 
higher  flux  neutron  generators  are  being  developed  [10],  but  these  technologies  are 
still  in  their  infancy  and  not  fully  realized. 
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Figure  1.  The  normalized  neutron  spectrum  of  three  different  types  of  neutron  gen¬ 
erators  shows  that  each  provides  a  different  spectral  shape  and  may  be  relevant  for 
different  applications  [11]. 

Another  example  of  an  accelerator-driven  source  is  the  one  found  at  the  Indiana 
University  (lU)  Center  for  the  Exploration  of  Energy  and  Matter  (CEEM)  Neutron 
Radiation  Effects  Facility  (NREF).  Here,  neutrons  are  created  via  the  inelastic  col¬ 
lision  of  protons  in  beryllium.  The  protons  are  accelerated  into  a  beryllium  target 
using  a  linear  accelerator.  The  resulting  neutrons  depend  strongly  on  the  incident 
proton  energy  and  can  be  very  energetic.  However,  because  the  amount  of  energy 
deposited  in  the  beryllium  target  is  so  great,  the  linear  accelerator  typically  has  a 
short  duty  time.  For  example,  at  the  NREF  target  at  CEEM,  the  standard  operating 
parameters  have  a  pulse  width  of  400  ps  with  a  repetition  rate  of  20  Hz.  This  leads  to 
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a  duty  cycle  of  approximately  0.8%,  which  is  driven  by  the  structure  of  the  beryllium 
target.  Advantages  to  this  source  include  the  ability  to  turn  the  source  on  and  off 
quickly,  relatively  low  gamma  contamination  (about  8-10  neutrons  for  every  gamma), 
and  a  hard  spectrum.  However,  a  disadvantage  is  the  relatively  low  time-averaged 
neutron  flux,  leading  to  long  times  of  continuous  operation  to  reach  high  fluences. 
If  the  community  is  to  use  this  type  of  neutron  production  source  for  displacement 
damage  studies,  the  effect  of  producing  damage  over  long  periods  of  pulsed  irradiation 
must  be  well  understood.  It  is  the  goal  of  this  research  to  further  this  understanding. 

1.2  Testing  Si  Bipolar  Transistors 

Si  BJTs  are  ubiquitous  semiconductor  devices  used  in  various  applications  in¬ 
cluding  power  transmission,  amplihcation,  and  high-freqnency  oscillation  [12].  There 
exists  a  signihcant  amonnt  of  data  collected  over  many  years  regarding  operation  of 
BJTs  in  ideal  and  in  operational  environments.  For  this  reason,  the  response  of  these 
devices  is  well-known  and  highly  characterized. 

Transistors,  specihcally  the  Si  2N2222A  n^pn  BJT,  are  being  nsed  instead  of  a 
novel  type  of  semicondnctor  device  snch  as  III-V  type  transistors  or  carbon  nanotnbe- 
based  devices  for  a  variety  of  reasons.  The  primary  reason  has  already  been  men¬ 
tioned;  that  is,  the  devices  are  very  well-characterized  in  a  large  nnmber  of  environ¬ 
mental  conditions.  Additionally,  because  these  devices  are  so  copious,  it  is  easy  to 
procure  a  large  lot  size  for  experimental  considerations.  By  using  the  same  lot  for  all 
experimental  procednres,  the  variability  inherent  in  lot-to-lot  variation  can  be  mini¬ 
mized  [13].  While  the  parts  tested  in  this  dissertation  research  were  obtained  from  a 
military-standard  procnred  lot,  transistor  variability  is  still  a  consideration  to  make 
when  developing  experimental  procednres.  An  additional  reason  for  using  Si  is  that  it 
is  an  indirect  bandgap  semicondnctor  material,  meaning  that  it  does  not  reqnire  Liq- 
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uid  Nitrogen  (LiN)  cooling  or  other  environmental  considerations  to  minimize  system 
electronic  noise.  For  all  of  these  reasons,  Si  n^pn  transistors  were  used  in  this  work. 
In  essence,  due  to  the  vast  knowledge  and  maturity  of  fabrication  techniques,  these 
devices  allow  the  researcher  to  explore  variations  in  effect  of  different  environments. 

1.3  Motivation 

The  sponsor  of  this  research  (supporting  the  author,  with  gratitude)  is  the  Naval 
Surface  Warfare  Center  (NSWC),  Crane  Detachment.  One  of  the  key  missions  of  the 
organization  is  to  analyze  the  performance  of  critical  electronic  systems  and  compo¬ 
nents  for  use  in  the  harsh  radiation  environment  produced  by  nuclear  weapons,  and 
to  verify  that  operational  performance  measures  are  met.  Due  to  the  complexity  of 
the  electronic  systems,  as  well  as  the  variability  and  complexity  of  the  environments, 
this  inherently  requires  a  substantial  test  base  covering  multiple  types  of  radiation 
over  a  wide  range  of  time  and  energy  scales. 

Sources  for  testing  and  analysis  are  varied,  expensive  and  specialized  for  the  nu¬ 
clear  weapons  environment,  due  to  the  extreme  fluences  and  short  time  of  interaction. 
Testing  sources  can  be  categorized  by  producing  ionization  effects  (e.g.  photons  and 
electrons)  and  displacements  effects  (neutrons  and  ions).  These  categorizations  are 
mainly  due  to  the  different  physical  mechanisms  of  the  radiation  interactions  and 
the  ultimate  effect  on  charge  transport  in  electronics;  ionizing  effects  can  be  tran¬ 
sient,  as  well  as  accumulative,  and  neutron  effects  are  generally  permanent.  The  test 
environments  are  simulated  based  upon  the  analysis  of  these  effects.  To  even  more 
hdelity,  testing  and  analysis  are  further  separated  by  endo-  and  exo-atmospheric  en¬ 
vironments.  This  separation  is  largely  due  to  the  effect  of  the  atmosphere  on  the  en¬ 
ergy  spectrum  of  the  radiation,  especially  regarding  neutrons.  For  endo-atmospheric 
weapons  sources,  the  neutrons  are  strongly  down  scattered  from  the  initial  hssion 
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spectrum  and  time-broadened  due  to  elastic  collisions  in  the  air  at  practical  dis¬ 
tances.  Thus,  a  test  environment  must  account  for  this  variation  in  damage  effects, 
and  it  is  what  makes  the  testing  and  validation  of  DoD  systems  much  more  rigorous 
that  those  for  space  operations  alone. 

Many  sources  exist  for  testing  and  analysis  of  pulsed  ionizing  effects,  primarily 
through  use  of  controlled  pulsed  electron  systems  (e.g  MEDUSA,  Utah  [14]).  For 
facilities  like  this,  the  delivery  of  high  fluxes  and  short  pulses  are  facilitated  by  the 
electron  charge,  which  allows  for  electronic  and  magnetic  field  steering  and  focus¬ 
ing.  NSWC,  Crane,  is  developing  an  advanced  dose  rate  facility  at  Indiana  Univer¬ 
sity  which  would  provide  a  local  source  for  conducting  ionization  irradiation  testing 
and  analysis;  the  Advanced  Electron-Photon  Facility  (ALPHA)  at  Indiana  Univer¬ 
sity  [15].  The  system  offers  high  fluxes  and  target  irradiation  sizes  that  not  only 
meet  the  NSWC  mission  requirements,  but  also  DoD  and  DoE  dose  rate  survivability 
requirements  with  a  narrow  energy  band  window  to  allow  for  precision  in  ionization 
radiation  effects  analysis.  Aside  from  the  operational  parameters,  the  near  co-location 
with  NSWC,  Crane,  will  likely  increase  interoperability  and  lower  overall  cost  per  shot 
for  testing. 

Sources  for  endo-atmospheric  testing  and  analysis  of  short  pulse  neutron  effects 
on  electronics  are  primarily  accomplished  at  the  Sandia  National  Laboratory  Annual 
Core  Research  Reactor  (ACRR).  This  system  is  suitable  for  such  testing  due  to  its 
long  pulse  period,  flux,  and  moderated  spectrum.  Sandia  has  many  years  worth  of 
testing  and  analysis  using  this  system  for  pulsed  neutron  analysis.  However,  the 
ACRR  is  not  a  panacea,  given  that  the  proximity  of  the  devices  to  the  reactor  does 
not  allow  for  neutron  only  testing,  and  actions  similar  to  those  used  in  this  research 
must  be  accomplished  to  separate  ionizing  and  non-ionizing  effects. 

Since  the  closing  of  the  Sandia  Pulsed  Reactor  HI  [6]  in  September  2007,  the 


Department  of  Defense  (DoD)  and  the  Department  of  Energy  (DoE)  have  put  much 
effort  into  developing  methods  for  simulating  the  exo-atmospheric  short-pulse  neutron 
environment  [16].  Although  there  have  been  proposals  for  development  of  large,  high 
flux,  high  energy  (up  to  14  MeV)  neutron  sources,  much  of  the  effort  has  gone  to 
developing  a  physics  based  experimental  approach  that  allows  for  the  use  of  modeling 
to  extrapolate  to  regions  of  interest  with  high  hdelity  and  conhdence.  During  that 
process,  NSWC  and  much  of  the  DoD  have  defaulted  to  the  use  of  White  Sands 
Missile  Range  Fast  Burst  Reactor  (FBR)  for  much  of  their  experimental  validation. 
Although  useful  for  validation,  it  has  similar  concerns  for  validation  and  testing  due 
to  flux,  timing,  and  size  restrictions,  as  well  as  the  travel  and  coordination  costs;  and 
the  future  of  the  facility  is  the  driving  concern. 

One  of  NSWCs  approaches  to  meet  their  exo-atmospheric  testing  needs  has  been 
oriented  on  development  of  a  neutron  source  that  is  relatively  local,  does  not  require 
special  handling  permits  for  special  nuclear  material,  and  can  meet  flux  and  timing 
requirements  for  validation  of  their  critical  components.  The  Neutron  Radiation 
Effects  Facility  (NREF)  at  Indiana  University  (described  in  more  detail  later  in  this 
document)  has  strong  promise  for  this  purpose,  and  has  many  positive  attributes 
similar  to  the  ALPHA.  The  spectrum  can  be  tailored  by  changing  the  proton  energy 
to  allow  for  better  hdelity  of  spectral  effects;  the  neutron  source  is  relatively  gamma 
free,  which  reduces  the  complexity  of  experiments;  the  pulses  have  a  fast  rise  and 
fall  time  allowing  for  temporal  variation  in  experiments;  the  system  does  not  require 
special  nuclear  materials;  and  it  is  located  close  to  NSWC.  It  suffers,  however,  from  a 
low  duty  cycle  (0.8%)  and  low  hux,  both  requiring  extensive  engineering  and  design 
to  improve.  Furthermore,  it  is  operated  in  a  continuous  multiple-pulse  mode,  possibly 
presenting  variations  in  physical  mechanisms  that  cannot  be  related  to  legacy  research 
due  in  part  to  variation  in  defect  formations  from  high  instantaneous  huxes,  and 
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short  term  annealing  effects  that  may  alter  the  outcome  during  the  off  time  of  the 
duty  cycle.  Experiments  at  this  facility  have  shown  that  a  detailed  understanding 
of  neutron  damage  in  Si  is  not  comprehensive  and  requires  additional  investigation. 
It  is  for  this  reason  that  this  research  is  undertaken,  so  that  decisions  and  efforts  to 
engineer  a  better  source  are  done  with  an  understanding  of  the  trade-offs  and  provide 
the  nation  with  the  simulation  capabilities  required  to  meet  mission  requirements. 

1.4  Problem  statement 

The  traditional  method  to  predict  the  response  of  bipolar  junction  transistors  that 
have  been  subjected  to  displacement  damage  from  neutrons  is  using  the  Messenger- 
Spratt  equation  [17].  However,  this  equation  is  usually  applied  to  short,  single-pulse 
irradiation  environments.  The  hrst  question  answered  in  this  dissertation  is:  does 
defect  annealing  during  irradiation  affect  the  end  result  and,  if  so,  can  it  be  corrected? 

Additionally,  the  neutron  spectrum  is  often  collapsed  down  into  an  equivalent 
mono-energetic  spectrum  with  a  corresponding  damage  constant.  The  second  ques¬ 
tion  asked  in  this  work  is  whether  the  neutron  spectrum  shape  will  affect  the  damage 
response  in  Si  BJTs. 

To  address  these  questions,  the  work  is  broken  up  into  three  pieces: 

1.  Develop  a  model  that  tracks  bulk  defect  species  concentrations  so  that  the 
current  gain  in  a  BJT  can  be  predicted  with  acceptable  accuracy  for  a  variety 
of  input  spectra  and  time  prohles. 

2.  Answer  the  question  of  whether  a  BJT  damaged  in  an  accelerator-driven  neutron 
source  facility  will  follow  the  Messenger-Spratt  response. 

3.  Collect  empirical  evidence  of  whether  or  not  the  spectral  shape  affects  the  re¬ 
sponse  by  using  different  facilities  to  irradiate  the  same  type  of  BJT. 
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II.  Theory 


In  this  portion  of  the  document,  the  underlying  theory  that  will  be  necessary 
to  understand  the  overarching  goals  of  the  research  will  be  presented.  As  such,  the 
main  focus  will  be  on  semiconductor  physics,  defect  formation  in  semiconductors, 
and  radiation  damage  mechanisms  in  silicon  devices.  In  addition,  a  description  of 
radiation  environments  and  electronics  hardness  testing  will  be  laid  out. 

2.1  Basic  Semiconductor  Physics 

Material  Classification. 

All  solid  materials  can  be  categorized  into  three  main  electrical  classes:  conduc¬ 
tors,  insulators,  and  semiconductors  [1].  These  materials  are  classihed  in  this  way 
based  on  a  number  of  electrical  properties,  one  of  which  is  the  conductivity.  Con¬ 
ductors  have  such  high  conductivity,  and  inherently  low  resistivity,  that  electrons 
can  easily  travel  through  the  material,  when  induced  to  do  so,  with  little  loss  of  en¬ 
ergy  [2].  Conversely,  insulators  have  extremely  low  conductivity  and,  thus,  electrons 
have  difficulty  moving  throughout  the  material.  In  the  middle,  the  conductivity  of 
semiconductors  is  strongly  dependent  on  available  carriers,  as  presented  by  doping 
concentration  or  applied  helds. 

Another  way  to  classify  crystalline  solids  is  by  the  band-gap  energy  between  the 
valence  and  conduction  bands,  which  is  presented  in  the  chart  in  Figure  2.  This 
gap  arises  because,  as  the  atoms  in  the  solid  form  bonds,  energy  bands  form  with  a 
forbidden  region  in  between.  Electrons  can  occupy  the  lower,  or  valence,  band  or, 
if  given  enough  energy,  they  may  be  found  in  the  higher,  or  conduction,  band.  The 
band-gap  energy.  Eg,  is  the  width  of  the  forbidden  region  between  these  two  bands. 

In  a  conductor,  the  conduction  band  may  be  partially  hlled  or  it  may  overlap 
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Figure  2.  A  schematic  of  the  energy  bands  for  a  conductor  (left)  with  a  partially  filled 
conduction  band  (top)  and  with  overlapping  valence  and  conduction  bands  (bottom), 
a  semiconductor  (middle),  and  an  insulator  (right).  (Figure  derived  from  [3].) 

with  the  valence  band,  and  electrons  can  easily  move  thronghont  the  solid  with  little 
resistance.  In  an  insnlator,  the  band-gap  energy  is  very  large,  so  electrons  need 
signihcant  energy  to  occupy  the  conduction  band.  A  semiconductor  occupies  the 
“goldilocks”  region;  the  band-gap  energy  is  small  enough  that  some  electrons  may 
be  excited  from  the  valence  band  to  the  conduction  band  by  thermal  energy  or  an 
applied  electric  held  but  not  so  small  that  all  do. 

When  an  electron  is  excited  into  the  conduction  band,  it  leaves  behind  a  hole  in 
the  valence  band.  Despite  physically  being  the  lack  of  something,  this  hole  is  treated 
as  a  particle  and  is  the  opposite  of  the  electron.  While  the  electron  (with  the  symbol 
e~)  is  considered  the  carrier  of  negative  charge,  the  hole  carries  positive  charge,  and  is 
often  delineated  by  the  symbol  h^ .  In  the  calculation  of  semiconductor  properties  and 
performance  parameters,  electrons  and  holes  are  the  primary  negative  and  positive 
charge  carriers,  respectively. 
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Doping  of  Semiconductor  Material. 


In  an  ideal  Si  semiconductor,  the  crystal  is  a  perfect  lattice  of  silicon  atoms, 
each  one  covalently  bonded  to  4  others.  This  is  called  the  intrinsic  semiconductor 
state.  In  this  state,  the  Si  carrier  concentration  is  approximately  10^°  cm“^  and  is 
approximately  equal  in  electrons  and  holes. 

The  introduction  of  impurities  into  the  semiconductor  changes  the  carrier  trans¬ 
port  characteristics.  The  impurities  must  lie  within  the  silicon  bandgap  in  order  to 
act  as  either  donors  or  acceptors.  For  a  donor,  the  impurity  electron  energy  typically 
lies  close  to  the  conduction  band  so  that,  with  very  little  thermal  energy,  almost 
100%  of  their  electrons  are  excited  into  the  conduction  band.  This  creates  a  signif¬ 
icant  increase  of  negative  charge  carriers  in  the  conduction  band.  Similarly,  for  an 
acceptor,  the  impurity  usually  lies  near  the  valance  band  and  accepts  electrons  from 
the  valence  band  with  very  little  thermal  energy  required,  which  subsequently  adds 
additional  holes  to  the  semiconductor  material.  Figured  displays  some  of  the  more 
common  impurities  added  to  silicon  semiconductors  as  well  as  their  position  within 
the  bandgap. 

When  at  thermal  equilibrium,  the  concentration  of  electrons  and  holes  in  an  n-type 
extrinsic  (doped)  semiconductor  is  described  by  the  equation 

nl  =  np=  {NcNv)  (1) 

where  Ui  is  the  intrinsic  carrier  concentration;  n  and  p  are  the  electron  and  hole 
concentrations,  respectively;  Nc  and  Ny  are  the  effective  density  of  states  in  the 
conduction  and  valence  bands,  respectively;  Eg  is  the  band-gap  energy  for  the  semi¬ 
conductor;  k  is  the  Boltzmann  constant;  and  T  is  the  temperature. 

When  a  semiconductor  material  is  doped,  its  intrinsic  qualities,  such  as  carrier 
mobility  and  diffusivity,  change. 
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Figure  3.  Measured  ionization  energies  for  select  acceptors  and  donors  in  Si  (in  eV). 
Below  the  Fermi  level,  impurities  are  acceptors  (unless  otherwise  noted);  above  the 
Fermi  level,  impurities  are  donors.  (Figure  derived  from  [3].) 


Semiconductor  Junctions. 

When  a  j)-type  region  is  brought  in  contact  with  an  n-type  region,  the  result  is 
called  a  junction.  At  this  junction,  if  it  is  abrupt,  a  region  of  space  charge  develops. 
In  his  1949  paper  [4],  Shockley  succinctly  describes  a  pn  junction.  In  this  work,  a 
simple  dopant  prohle  leads  to  two  separate  regions  (p-  and  n-type)  with  a  transition 
region  in  between. 

Due  to  diffusion  and  the  spatial  variation  of  charge  distribution,  electrons  in  the  n- 
type  material  drift  to  towards  the  p-type  just  as  holes  move  in  the  opposite  direction. 
The  width  of  this  region  is  highly  dependent  on  quantities  such  as  the  built-in  voltage, 
Vbi,  and  the  dopant  concentration  in  the  lighter-doped  region,  Nb-  The  depletion 
width,  W,  is  described  by 
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where  is  the  permittivity,  V  is  an  applied  external  voltage,  and  q  is  the  elementary 
charge  constant  (1.602  x  10“^®  C).  In  real  world  applications  of  pn-junctions,  the 
abrupt  junction  is  not  possible;  however,  it  is  often  used  as  an  approximation  to 
simplify  calculations. 

The  junction  is  a  key  inspection  point,  as  much  of  the  focus  on  displacement 
damage  in  a  silicon  bipolar  junction  transistor  is  put  on  the  emitter-base  junction. 
Defects  in  this  region  can  affect  the  diffusion  of  charge  carriers,  which  in  turn  affects 
the  device’s  performance.  This  will  be  discussed  in  greater  detail  further  in  this 
document. 


Semiconductor  Device  Types  &;  Applications. 

Elemental  semiconductors  such  as  germanium  and  silicon  are  most  often  used 
as  radiation  detectors.  Additionally,  germanium  is  used  extensively  as  a  gamma 
spectrometer.  However,  because  of  the  small  band-gap  energy,  it  must  be  operated 
at  low  temperatures  (~77K)  to  minimize  electronic  noise.  Because  of  the  small  band- 
gap  energy,  thermal  energy  at  room  temperature  provides  enough  energy  to  electrons 
in  the  valence  band  that  they  jump  to  the  conduction  band,  providing  a  distorted 
signal.  Silicon  strip  detectors  are  the  choice  of  high-energy  physics  projects  such  as 
the  Linear  Hadron  Collider  (LHC).  The  larger  band-gap  energy  and  the  fact  that 
silicon  is  an  indirect  semiconductor,  which  means  that  energy  and  momentum  must 
be  applied  to  an  electron  for  excitation  to  the  conduction  band,  leads  to  low-noise 
operation  at  room  temperature. 

More  often,  semiconductor  devices  are  created  with  junctions  between  p-  and  n- 
type  regions.  To  be  considered  one  or  the  other,  the  semiconductor  region  must  be 
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doped  with  donors  or  acceptors  so  that  there  are  more  electrons  in  the  conduction 
band  (n-type)  or  more  holes  in  the  valence  band  (p-type)  than  are  expected  in  an 
intrinsic  semiconductor.  Figured  displays  select  ionization  energies  for  both  donors 
and  acceptors  in  silicon. 

The  pn-junction  is  the  fundamental  building  block  for  the  BJT.  The  BJT  has  been 
used  extensively  in  high-speed  analog  and  digital  circuits,  predominately  as  a  current 
ampliher.  These  devices  typically  have  three  regions:  the  collector,  the  base,  and  the 
emitter;  and  these  regions  interact  through  the  pn  or  np  junctions.  In  the  case  of  a 
silicon  p'^np  transistor,  the  substrate  is  grown  as  lightly-doped  p-type  and  acts  as  the 
collector.  Next,  n-type  dopants  are  diffused  into  the  substrate  to  create  an  n-type 
base  region.  Finally,  a  heavily-doped  p-type  region  is  diffused  into  the  n-type  base 
to  create  the  emitter.  Metal  contacts  are  attached  at  each  region  to  facilitate  carrier 
flow.  A  schematic  of  an  idealized  n+pn  BJT  is  displayed  in  Figured.  For  an  p'^np 
BJT,  regions  of  n-type  doping  are  replaced  by  regions  of  p-type,  and  vice  versa. 


COLLECTOR 


Figure  4.  A  schematic  view  of  an  idealized  ri^pn  bipolar  junction  transistor.  The  silicon 
dioxide  is  used  to  passivate  the  surface  boundary.  (Figure  derived  from  [3].) 


In  a  BJT,  charge  carriers  are  classihed  as  either  the  majority  or  minority  carrier 
based  on  the  extrinsic  state  of  the  device  region.  So,  for  example,  in  the  emitter  region 
of  an  n'^pn  BJT,  which  is  n-type  and  doped  with  elements  such  as  P  or  arsenic,  the 
majority  carrier  is  the  electron  while  the  minority  carrier  is  the  hole.  However,  as  the 
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electron  travels  from  the  emitter  into  the  base,  which  is  j)-type  and  may  be  doped 
with  elements  such  as  boron  or  gallium,  it  becomes  the  minority  carrier. 

When  in  thermal  and  electronic  equilibrium,  i.e.  room  temperature  and  no  volt¬ 
age  bias  applied  to  the  device  leads,  the  band  structure  of  an  BJT  would  be 

represented  as  in  Figured.  There  is  band  bending  at  the  junctions  due  to  the  effect 
of  space  charge  in  the  depletion  region.  The  Fermi  level,  i?/  is  constant  and  its  re¬ 
lationship  with  respect  to  the  conduction  and  valence  band  is  based  on  the  carrier 
concentrations,  which  are  affected  by  dopant  levels  in  each  region. 
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Figure  5.  The  Fermi  level,  Ef,  is  continuous  across  the  entire  bandgap  schematic  of  an 
idealized  rFpn  BJT  in  thermal  and  electronic  equilibrium. 

If  a  potential  is  applied  to  various  regions  of  the  BJT,  the  Fermi  level  is  perturbed. 
An  example  of  the  iiFpn  BJT  from  Figured  with  positive  bias  applied  to  the  collector 
and  negative  bias  applied  to  the  emitter  is  displayed  in  Figured.  The  positive  po¬ 
tential  applied  to  the  n-type  collector  region  causes  the  overall  potential  to  decrease 
and  puts  the  base-collector  junction  in  reverse  bias.  Alternatively,  the  negative  bias 
applied  to  the  emitter  causes  the  opposite  effect,  causing  the  base-emitter  junction  to 
operate  in  forward  bias.  For  this  particular  configuration,  the  BJT  is  in  aetive  mode. 
A  table  of  the  four  BJT  modes  is  displayed  in  Table  1. 
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Table  1.  Modes  of  BJT  operation. 


Base-Emitter 
Junction  Bias 

Base-Collector 
Jnnction  Bias 

Applied 

Voltages 

Mode 

Forward 

Reverse 

E  <B  <C 

Forward-Active 

Forward 

Forward 

E  <B  >C 

Saturation 

Reverse 

Reverse 

E  >B  >C 

Cut-Off 

Reverse 

Forward 

E  >B  >C 

Reverse-Active 

Emitter  Base  Collector 


Figure  6.  The  BJT  bandgap  shifts  as  bias  is  applied.  In  this  case,  postive  voltage  is 
applied  at  the  emitter,  and  negative  voltage  is  applied  to  the  collector  while  the  base 
is  kept  at  ground. 


The  BJT,  specifically  the  current  amplihcation  BJT  with  the  designation  2n2222A, 
is  the  semiconductor  device  that  will  be  studied  extensively  in  this  research.  For  non¬ 
ionizing  radiation,  the  base  region  is  of  greatest  interest.  As  early  as  1964,  Goben  and 
his  colleagues  reported  a  neutron-induced  base  current  component  which  increases 
in  proportion  to  neutron  fluence  [5-7].  This  base  current  component  is  primarily 
responsible  for  gain  degradation  because  of  degradation  of  the  emitter  efficiency  [8]. 

Current  gain  is  one  of  the  most  important  device  properties  that  can  be  measured 
for  the  2n2222A  BJT.  When  in  active  mode,  electrons  from  the  emitter  are  injected 
into  the  base  since  the  emitter-base  junction  is  forward  biased.  If  the  base  width 
is  small  compared  with  the  diffusion  length  of  the  minority  carriers  (in  this  case, 
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electrons),  most  of  the  electrons  survive  and  make  it  to  the  collector  region  [9].  This 
is  often  done  intentionally  to  maximize  carrier  transport  properties.  The  current  gain, 
which  can  be  called  f3  or  hfe,  is  dehned  as  (3  =  hfe  =  Ic/h  where  Ic  is  the  collector 
current  and  h  is  the  base  current. 

To  help  further  dehne  the  current  gain.  Figure  7  displays  a  schematic  of  the  cur¬ 
rents  inside  a  BJT.  The  forward  biased  junction  at  the  emitter-base  interface  injects 
electrons  into  the  base  where  they  diffuse,  and  this  injected  current  of  electrons  is 
named  /e„.  Additionally,  there  is  a  current  of  holes  in  the  emitter,  Aj,  that,  when  con¬ 
sidered  together  with  the  electron  current,  makes  up  the  total  current  in  the  emitter. 
The  electrons  that  do  not  recombine,  surviving  to  the  collector,  make  up  the  current 
of  electrons  in  the  collector,  or 


EMITTER  BASE  COLLECTOR 


Figure  7.  Idealized  schematic  of  electron  and  hole  currents  in  an  n+pn  BJT. 

To  calculate  hf^,  the  base  transport  factor  (ap)  and  the  emitter  injection  efficiency 
(7)  are  necessary,  as  displayed  in  Equation  Set  3. 
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In  Equation  Set  3,  PnEo  and  UpB^  are  the  minority  carrier  concentrations  for  emitter 
and  base,  respectively;  Dpj^  and  Dng  are  the  diffusivity  of  the  minority  carriers  for 
emitter  and  base,  respectively;  L„  is  the  diffusion  length  of  electrons  in  the  base;  and 
Wh  and  We  are  the  physical  widths  of  emitter  and  base,  respectively.  The  diffusion 
length  is  a  measure  of  the  average  distance  that  an  electron  can  travel  from  birth  to 
recombination  and  is  dehned  as  L  =  D^gT,  where  r  is  the  minority  carrier  lifetime. 

The  current  gain  in  a  semiconductor  is  then  dehned  as 
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It  is  (3  that  will  be  the  main  focus  of  the  empirical  data  collection  efforts  in  this 
research.  Since  a  number  of  the  parameters  that  go  into  calculating  it  are  affected  by 
non-ionizing  radiation,  it  will  provide  the  best  metric  of  damage  in  the  2N2222A  Si 
BJTs. 


2.2  Semiconductor  Defects  in  Silicon 

At  the  most  basic  level,  the  primary  defect  mode  in  solid  state  semiconductors  is 
the  Frenkel  pair  [1] .  This  type  of  defect  is  created  when  enough  energy  is  transferred 
to  an  atom  in  the  crystal  lattice  of  the  material  so  that  it  is  relocated  from  its  initial 
location.  Often,  this  energy  is  deposited  in  the  crystal  lattice  via  incident  particle 
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collisions,  as  in  the  case  of  nentrons  colliding  with  silicon  atoms  in  a  single  crystal. 
This  relocation  canses  a  vacancy  defect  at  the  initial  location.  Additionally,  the 
displaced  atom  becomes  an  interstitial  defect  if  its  hnal  resting  place  is  not  at  the 
appropriate  lattice-aligned  location.  The  defects  thns  created  become  charge  carrier 
traps  that  may  act  as  recombination  sites  for  electrons  and  holes  if  they  lie  within 
the  bandgap.  In  essence,  these  defects  become  allowed  states  within  the  forbidden 
region  of  the  semicondnctor  bandgap,  increasing  the  possibility  of  recombination  as 
electrons  and  holes  drift  and  diffnse  throngh  the  bnlk  material.  Figure  8  is  a  schematic 
of  the  Frenkel  pair  creation  that  illustrates  this  process  using  neutron  irradiation  of 
the  crystal  lattice. 
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Figure  8.  Schematic  of  the  creation  of  a  Frenkel  pair  (in  this  case  using  low-energy 
neutron  bombardment  of  the  crystal  lattice).  [10] 


As  the  bombarding  particle  travels  into  the  material,  it  collides  with  a  stationary 
lattice  atom,  imparting  energy  and  dislodging  the  atom.  The  electron  may  be  ab¬ 
sorbed  or  may  continue  on  through  the  bulk,  potentially  colliding  and  causing  other 
dislocations.  The  original  atom  is  bumped  out  of  place  (diagonal  arrow  in  Figure  8) 
where  it  replaces  the  atom  at  that  site.  This  process  may  continue  many  times  until 
the  hnal  atom  comes  to  rest  in  between  lattice  sites  (hlled  in  circle  in  Figure  8).  The 
time  scale  for  this  process  is  short  and  is  often  less  than  10“’^  seconds  [11]. 

In  addition  to  Frenkel  pairs,  there  are  other  combinations  of  defects  that  are 
known  by  other  names.  Two  vacancies  adjacent  to  one  another  form  a  di-vacancy. 
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It  is  possible,  but  less  common,  for  the  same  phenomenon  to  occur  with  interstitials, 
forming  a  di-interstitial.  Also,  it  has  been  found  that  vacancies  and  interstitials 
in  silicon  act  as  drivers  for  various  diffusion  mechanisms;  they  can  interact  with 
impurities  in  the  crystal  to  form  complex  defects  [12],  One  common  example  is  the 
vacancy-phosphorus  pair,  known  as  the  E  center  in  irradiated  Si  [13].  This  is  quite 
common  in  n-type  semiconductor  material,  as  P  is  a  very  common  donor  atom.  The 
result  of  defect  formation  in  semiconductors  is  the  insertion  of  allowable  energies  in 
the  forbidden  bandgap. 

Table  2  presents  a  small  selection  of  the  more  common  defects  as  collected  by 
Myers  et  al.  [14].  In  the  table,  AE  is  the  activation  barrier  for  electron  and  hole 
emission,  a  is  the  capture  cross  section  for  electrons  and  holes.  Do  is  the  diffusion 
prefactor,  and  is  the  diffusion  activation  energy.  In  their  study,  Myers  used  this 
information  to  calculate  clustering  in  pulsed-neutron  irradiated  silicon.  The  list  of 
defects  selected  here  displays  some  of  the  most  common  defects,  but  it  is  a  very  small 
portion  of  the  total  list.  It  is  included  here  to  show  the  complexity  of  defect  formation. 


Table  2.  Select  defect  species  and  their  electronic  parameters.  Here  e~  denotes  the 
conduction  electron  and  denotes  the  hole,  V  is  the  vacancy,  the  subscript  /  denotes  an 
interstitial  atom,  and  charge  states  are  indicated  by  superscripts.  (Derived  from  [14]) 


AE‘ 

AE'‘ 

Species 

[eV 

[eV] 

[cm^] 

[cm^] 

Do 

E"* 

Si? 

0.62 

0.86 

3 

X  10-1^ 

3 

X  10-1® 

o 

1 

CO 

0.17 

Si+i 

0.47 

0.50 

3 

X  10-1^ 

3 

X  10-1® 

O 

1 

CO 

0.50 

Si7^ 

0.26 

1.01 

3 

X  10-1® 

3 

X  10-1^ 

10-® 

0.29 

1.07 

0.72 

3 

X  10-1® 

3 

X  10-1® 

1.3  X  10-3 

0.45 

vyo 

0.91 

0.71 

3 

X  10-1® 

3 

X  10-1® 

B? 

0.13 

0.75 

3 

X  10-1® 

3 

X  10-1® 

2.3  X  10-® 

0.53 

VP° 

0.85 

0.68 

3 

X  10-1® 

3 

X  10-1® 

Shi  et  al.  [15]  found  that  initial  defect  formation  is  over  on  a  nanosecond  time 
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scale.  In  their  work,  they  determined  that  after  about  10“^°  s,  recombination  and 
divacancy  production  are  independent  of  a  production  constant,  indicating  that  most 
of  the  interstitials  and  vacancies  have  annihilated.  It  is  for  this  reason  that  the  model 
that  was  developed  in  this  research  only  included  a  generic  generation  term  since  the 
initial  movement  of  silicon  ions  is  over  on  such  short  time  scales.  The  resulting  defects 
then  migrate  under  different  mechanisms,  and  this  behavior  is  the  focus  of  the  model. 

Defect  Clusters. 

The  most  common  damage  mode  in  a  fast-neutron  irradiation  is  the  defect  cluster. 
While  they  still  are  not  fully  understood,  there  is  extensive  experimental  data  proving 
their  existence  including  transmission  electron  microscopy  (TEM)  observations.  They 
have  been  observed  for  silicon  ion  bombardment  [16],  medium  [17]  and  heavy  [18,19] 
ions,  protons  [20],  and  electrons  [21]. 

In  this  mode,  a  signihcant  amount  of  kinetic  energy  is  deposited  in  the  Primary 
Knock-on  Atom  (PKA).  The  PKA  causes  secondary  and  tertiary  collisions  in  the 
surrounding  lattice,  depositing  and  dissipating  energy  in  the  form  of  thermal  energy 
[22].  In  Figure  9  is  displayed  a  representative  track  for  a  PKA  given  a  starting  kinetic 
energy  of  50  keV  at  the  origin.  There  are  a  few  very  energetic  collisions  that  give  rise 
to  other  high-energy  recoils,  but  low-energy  transfers  are  more  common.  The  lower- 
energy  transfers  often  create  Frenkel  pairs  that  coalesce  into  stable  defect  clusters  [15] . 
Each  track  path  ends  in  a  terminal  cluster  [23]. 

Gossick  [25]  proposed  a  physical  model  regarding  defect  clusters  that  involved  a 
disordered  volume  surrounded  by  a  depletion  region.  This  occurs  because  the  highly 
disordered  region  (defect  cluster)  captures  majority  carriers  from  the  surrounding 
bulk,  creating  opposing  space  charge.  Subsequent  work  from  Bertolotti  observed  this 
phenomenon  using  electron  microscopy  and  etching  [26]. 


25 


Silicon  buDc 
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Figure  9.  Recreation  of  a  typical  history  track  of  a  50  keV  recoil  atom  in  silicon 
produced  using  a  Monte  Carlo  simulation.  [24] 

Defect  formation  is  thought  to  be  energy  dependent,  meaning  that  the  number  of 
defects  and  the  number  of  clusters  that  form  depends  on  the  PKA  recoil  energy.  If  a 
neutron  deposits  barely  sufficient  energy  to  cause  a  Frenkel  pair,  called  the  threshold 
energy  for  displacement,  a  simple  defect  forms.  The  threshold  energy  for  silicon  is 
about  20  eV  [28] .  However,  for  a  fast  neutron  (one  with  higher  energy)  that  deposits 
its  energy  to  a  silicon  atom,  it  can  create  many  thousands  of  defects  which  cluster  as 
shown  in  Figure 9.  This  damage  anneals  in  the  short-term  (up  to  about  10^  seconds), 
as  well  as  long-term. 

Short  Term  Annealing. 

The  displacement  damage  caused  in  a  BJT  subjected  to  neutron  irradiation  under 
standard  temperature  (300K)  usually  finishes  short-term  annealing  on  the  order  of 
10^  seconds.  Before  this  time,  the  clusters  rearrange  due  to  both  thermal  effects 
and  injection  annealing.  It  is  important  to  note  that,  during  this  time,  the  damage 
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Figure  10.  Schematic  representation  of  defect  and  subcascade  formation  as  a  function 
of  PKA  energy.  The  y-axis  of  the  plot  is  logarithmic  and  representative  of  the  number 
of  interactions  [27]. 


constant  is  greater  than  for  when  the  device  is  stable. 

One  measure  of  damage  in  the  short  term,  i.e.  before  reaching  the  stable  damage 
condition,  is  the  annealing  factor  defined  by  Sander  and  Gregory  [29-32]  as 


AF{t) 


-  /So-' 


(5) 


In  Equation  5,  /3  is  the  device  gain,  is  the  pre- irradiation  gain,  and  /Soo  is  the  gain  “at 
inhnite  time”  which  is  usually  simulated  by  performing  an  80°C  thermal  annealing 
bake  for  2  hours.  Binder  et  ah  [33]  conjectured  that  the  short  term  annealing  is 
governed  by  a  second  order  recombination  of  mobile  and  immobile  defects  inside  a 
spherical  defect  cluster.  The  rapid  recombination  reaction  is  slowed  down  by  diffusion 
out  of  the  cluster. 

The  annealing  factor  is  affected  by  temperature  [34],  injection  level  [35],  and 
neutron  energy  levels.  The  short  term  annealing  process  is  thermally  activated,  with 
longer  decay  times  at  lower  temperatures  [36].  This  is  due  to  the  extra  energy  im¬ 
parted  to  the  lattice  at  higher  temperatures  and  increases  vibrational  energy,  which 
rearranges  the  defects  faster. 

In  p-type  silicon,  high  current  injection  causes  the  defects  to  anneal  out  of  the 
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device  faster.  Conversely,  in  n-type  silicon,  high  injection  causes  annealing  to  proceed 
more  slowly  [37].  The  difference  in  annealing  factor  due  to  material  type  is  attributed 
to  higher  vacancy  mobility  in  its  negative  charge  state  [24].  This  will  have  a  great 
affect  on  the  results  as  the  in  situ  measurements  will  require  the  devices  to  be  powered 
on,  increasing  current  injection. 

Finally,  the  annealing  factor  depends  on  the  energy  of  the  neutrons  that  are  irra¬ 
diating  the  silicon  device.  Srour  and  Curtis  [38]  found  that  the  ratio  of  the  annealing 
factor  for  14  MeV  neutrons  to  thermal  neutrons  at  10“^  s  was  about  1.25.  However, 
at  about  1  s,  this  difference  is  negligible. 

In  1981,  McMurray  and  Messenger  [35]  £t  the  device  annealing  factor  to  an  an¬ 
alytical  expression  that  was  derived  from  an  empirical  nomograph.  This  equation  is 
stated  by  Wrobel  and  Evans  [39]  as  Equation  6. 
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where 


t  =  time  to  device  turn  on  [s], 
r  =  rise  time  of  the  injection  signal  [s]. 

Hi  =  intrinsic  carrier  concentration, 
q  =  elementary  charge  constant  =  1.62  x  10“^®[C'], 
k  =  Boltzmann’s  constant,  8.62  x  ld~^[ev / K], 

T  =  temperature  of  operation  [K], 

Vheif)  =  base-emitter  voltage  [V]. 


(6) 


It  must  be  noted  that  this  work  was  done  using  a  high-flux,  single  pulse,  fast-burst 


reactor  based  neutron  source.  This  will  be  compared  to  the  results  of  such  analysis 
done  in  different  facility  types  to  see  how  the  annealing  factor  changes. 

2.3  Radiation  Environments  and  Hardness  Testing 

Semiconductors  are  used  in  a  wide  variety  of  applications  that  lead  to  their  use 
in  an  equally  wide  variety  of  environments.  Extremes  of  temperature,  humidity, 
pressure,  and  other  environmental  variables  are  not  uncommon  for  their  operation. 
Radiation  environments  are  no  exception.  Take,  for  example,  the  environment  on 
a  space-borne  satellite.  Charged  particles  from  the  Earth’s  Van  Allen  belt,  gamma 
photons  and  protons  from  the  sun,  and  exotic  cosmic  particles  from  space  are  a 
constant  radiation  source,  which  is  time  dependent. 

For  this  reason,  radiation  hardness  testing  must  be  accomplished  to  determine  the 
range  of  adverse  conditions  under  which  the  devices  will  continue  to  operate.  This 
testing  may  take  many  forms,  depending  on  the  device  and  environmental  expecta¬ 
tions.  No  matter  the  form,  however,  the  testing  is  necessary  to  completely  understand 
the  operational  limits  of  the  devices  in  question. 

Obviously,  if  operations  in  extreme  environments  are  being  tested,  it  is  not  always 
possible  to  test  in  the  exact  environment.  As  an  example,  the  South  Atlantic  Anomaly 
(SAA)  is  a  region  of  the  world  where  the  Earth’s  magnetic  held  and,  thus,  the  Van 
Allen  belts  dip  to  less  than  200  km  from  the  surface  [40].  For  hscal  and  logistical 
reasons,  it  would  not  be  viable  to  test  every  lot  of  BJTs  created  for  use  in  a  satellite 
by  sending  a  rocket  into  orbit  with  a  sample  set  in  its  cargo  bays  and  bound  for  the 
SAA.  In  addition,  the  devices  often  need  to  be  tested  for  longevity,  e.g.  spending 
many  years  passing  through  the  SAA.  Experimenters  could  not  wait  for  each  batch 
to  run  their  20  year  lifetime  before  sending  the  actual  satellite  into  space.  Therefore, 
the  test  environment  often  needs  to  be  accelerated  by  increasing  the  hux  to  reach  the 


29 


expected  fluence  faster.  More  feasible  test-beds  must  be  developed  so  that  multiple, 
repeatable  experiments  can  be  undertaken  to  test  for  the  environmental  hardness. 

Another  key  environment  for  military  and  civilian  uses  of  transistor  components  is 
that  of  a  nuclear  weapon.  This  environment  is  extremely  harsh  with  a  short  duration 
(but  high  flux)  pulse  of  neutrons,  gammas,  x-rays,  and  charged  particles.  Interna¬ 
tional  law  forbids  detonating  nuclear  weapons  for  verihcation,  thus  there  is  a  strong 
desire  to  hnd  similar  or  comperable  radiation  testing  environments  for  verhciation 
testing.  Additionally,  the  neutron  flux  and  spectrum  shape  are  unique  to  the  nuclear 
weapon  environment,  making  it  difficult  to  replicate. 

In  this  work,  only  radiation  hardness  is  investigated.  Therefore,  in  the  following 
sections,  the  testing  for  neutron  displacement  damage  studies  will  be  discussed  in 
detail,  and  a  few  of  the  facilities  that  are  used  for  neutron  radiation  testing  will  be 
described. 

Description  of  hardness  testing. 

Radiation  hardness  testing  is  the  act  of  evaluating  a  device  or  piece  of  equipment  to 
understand  the  influence  of  radiation  on  its  performance.  The  environment  may  single 
out  a  particular  radiation  or  particle  type,  or  it  may  be  a  test  bed  mixed  radiation 
environment  designed  to  more  closely  replicate  the  environment  of  interest,  and  thus, 
the  real-world  application  of  the  device.  These  types  of  testing  are  done  for  different 
reasons;  often  they  are  done  concurrently  to  fully  understand  device  performance.  By 
singling  out  a  particular  type  of  radiation,  the  experiment  can  provide  information 
related  to  a  particular  damage  mechanism.  Otherwise,  by  including  as  many  particle 
types  as  the  environment  of  interest,  the  experiment  may  create  coupling  effects  not 
measured  when  singling  out  particular  particles.  In  some  cases,  the  method  by  which 
the  particle  of  interest  is  created  at  the  facility  leads  to  the  creation  of  secondary 
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particles.  These,  and  many  other  considerations  must  be  made  when  selecting  a 
facility  used  for  radiation  hardness  testing. 

Radiation  Testing  Facilities. 

There  are  a  variety  of  radiation  testing  facilities  that  have  been  developed  to 
test  electronic  devices  and  systems  of  devices.  They  can  be  classihed  by  the  type  of 
particle(s),  radiation(s),  or  energy  of  radiation. 

For  charged  particles,  cyclotrons,  linear  accelerators,  and  their  variations  are  most 
common.  An  example  includes  the  Indiana  University  Cyclotron  Facility,  a  schematic 
of  which  is  found  in  Figure  11.  Charged  particles  are  relatively  easy  to  work  with, 
since  they  can  be  controlled  and  manipulated  using  electric  or  magnetic  helds. 


Figure  11.  A  schematic  of  the  lUCF  cyclotron  is  a  proton  accelerator  that  can  provide 
20  ns  pulses  of  200  MeV  protons  at  about  2  x  10^°  rad(Si)/s  [41]. 

For  a  gamma  only  environment,  facilities  typically  use  radioactive  materials  that 
decay  by  gamma  emission,  often  ®°Co.  By  choice  of  radioactive  material  and  ra¬ 
dioactivity,  the  experimenter  can  select  the  important  energy  and  dose  rate  desired. 
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Mixed-radiation  fields  that  include  gamma  rays  result  from  reactions  that  provide 
direct  or  indirect  emission  of  gamma  photons.  As  an  example,  in  a  nuclear  reactor, 
2-3  neutrons  are  produced  for  each  fission  of  .  In  addition  to  these  2-3  neutrons, 
energy  is  released  in  the  form  of  gamma  rays.  The  fission  products  are  generally  ra¬ 
dioactive  and  often  decay  by  gamma  emission  as  well.  In  nuclear  reactors,  the  gamma 
dose  rate  can  be  8-10  times  that  of  the  neutron  dose. 

There  are  also  many  facilities  that  are  capable  of  creating  exotic  particles,  such  as 
positrons  or  muons.  These  will  not  be  discussed  here,  but  they  are  prevalent  in  the 
academic  community.  They  are  not  of  much  interest  for  radiation  effects  experiments. 

It  is  very  difficult  for  a  neutron  irradiation  facility  to  provide  an  environment  free  of 
radiation  other  than  neutrons.  Therefore,  for  neutron  irradiation  testing,  mitigation 
techniques  must  be  used.  In  this  research,  two  different  facility  types  were  used  to 
induce  displacement  damage  in  silicon  devices:  the  thermal  neutron  research  reactor 
at  the  Ohio  State  Research  Reactor  (OSURR)  and  the  lU  Neutron  Radiation  Effects 
Facility  (NREF).  They  are  described  in  greater  detail  in  the  following  sections. 

Ohio  State  University  Research  Reactor. 

This  system  consists  of  a  light  water,  pool-type,  thermal  research  reactor.  It  is 
licensed  to  operate  up  to  500  kW  and,  at  maximum  steady-state  power,  the  average 
thermal  neutron  flux  is  approximately  5  x  10^^  n/cm^/s. 

There  are  a  number  of  separate  sub-facilities  within  the  reactor  core,  displayed 
in  Figure  12,  that  allow  for  different  spectra  and  flux  rates.  The  Central  Irradiation 
Facility  (GIF)  has  an  inner  diameter  of  1.3  inches  and  extends  from  the  top  of  the  pool 
down  to  a  central  grid  position  within  the  core  [42].  It  has  a  maximum  total  flux  of 
2.3  X  10^^  n/cm^/s  and  a  maximum  1-MeV  Equivalent  (Si)  flux  of  4.52  x  10^^  n/cm^/s 
with  a  hardness  oi  H  =  0.20.  The  calculation  of  1-MeV  Equivalent  (Si)  flux  and  the 
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Figure  12.  Reactor  Core  at  the  OSURR  shows  multiple  irradiation  facilities  (CIF,  AIF, 
PIF)  as  well  as  beam  ports  for  external  set-up  and  irradiation. 

determination  of  hardness  is  accomplished  using  the  method  described  in  Section  2.3. 
The  Auxiliary  Irradiation  Facility  (AIF)  is  larger,  at  2.5  inches  in  diameter,  but  the 
maximum  fluxes  are  lower  at  (f)tot  =  9.4  x  10^^  n/cm^/s  and  (pi-MeV  =  2.44  x  10^^ 
n/cm^/s  with  a  hardness  of  0.27.  An  advantage  of  this  is  that  the  gamma  flux  is  also 
lower.  The  Peripheral  Irradiation  Facility  (PIF)  is  also  2.5  inches  in  diameter  and  in 
a  different  position,  with  a  maximum  total  flux  of  5.3  x  10^^  n/cm^/s. 

Finally,  a  7  inch  tube  can  be  moved  anywhere  in  the  reactor  pool,  but  the  stan¬ 
dard  position  allows  a  maximum  total  flux  of  1.5  x  10^^  n/cm^/s  and  a  1-MeV  Equiv¬ 
alent  (Si)  flux  of  1.92  X  10^^  n/cm^/s  with  a  hardness  of  0.13  and  a  gamma  dose  of 
6.5  X  10^  rad(Si) /hr.  The  maximum  flux  and  hardness  decreases  because  the  7  inch 
tube  is  positioned  the  furthest  from  the  center  of  the  core  when  compared  with  the 
other  tubes.  The  additional  distance  leads  to  additional  neutron  thermalization  as 
neutrons  travel  through  the  cooling  water.  The  spectrum  has  been  calculated  using 
Monte  Carlo  simulations,  and  the  calculations  have  been  verihed  using  foil  activation 
analysis  and  spectrum  unfolding  using  the  Spectrum  Analysis  by  Neutron  Detectors 
(SAND-II)  code  [43].  The  spectrum  in  the  7  inch  tube  is  displayed  in  Figure 54. 
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Figure  13.  Neutron  differential  flux  spectrum  at  the  OSURR  7  inch  tube  at  full  power 
where  the  1-MeV  equivalent  (Si)  flux  at  this  power  level  is  1.92  x  10^^  n/cm^/s  and  the 
hardness  is  0.13. 


Indiana  University  Neutron  Radiation  Effects  Facility. 

At  the  lU  CEEM,  there  is  a  pulsed  neutron  source  that  produces  neutrons  via  the 
inelastic  scattering  of  protons  on  beryllium.  One  of  the  two  targets  is  used  for  neutron 
radiation  effects  in  electronics  studies,  and  it  is  called  the  NREF.  A  schematic  of  the 
facility  is  displayed  in  Figure  14. 

The  following  description  is  taken  from  the  Master’s  thesis  of  M.R.  Halstead  [44]: 

Protons  are  created  in  an  ion  source,  then  are  accelerated  and  directed  onto  a 
target  assembly  inside  a  moderator/reflector  stack  where  neutrons  are  created.  Each 
of  these  steps  will  now  be  explained  in  greater  detail. 

The  protons  are  created  by  stripping  an  electron  from  elemental  hydrogen  in  the 
ion  source,  displayed  in  the  lower  left  corner  of  Figure  14.  These  protons  are  collected 
and  focused  into  a  3  MeV  Radio  Freqnency  Qnadrupole  (RFQ)  accelerator  directly 
conpled  to  two  Drift  Tube  LINACs  (DTL),  which  accelerate  protons  to  their  hnal 
energy  of  13  MeV.  Proton  acceleration  is  facilitated  by  three  Litton  5773  klystron 
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Figure  14.  A  schematic  of  the  lU  NREF  beamline  at  the  LENS. 


RF  tubes  that  provide  425  MHz  microwave  power  to  the  RFQ  and  DTL  components. 
Power  delivered  to  the  LINAC  is  approximately  1  MW  with  a  0.9%  duty  factor  [45]. 
The  duty  factor  is  a  measure  of  the  time  a  pulsed  beam  is  running.  In  essence,  the 
duty  factor  can  be  described  by  DF  =  where  PW  is  the  pulse  width  and  tcyde  is 
the  total  time  for  one  pulse  cycle,  which  includes  one  proton  pulse  and  one  recovery 
period,  and  is  defined  as  tcyde  =  where  /  is  the  accelerator  frequency.  The  typical 
beam  operating  frequency  is  20  Hz,  making  tcyde  =  50  ms. 

After  acceleration,  the  protons  are  focused  using  octupole  magnets  and  directed 
down  one  of  two  beam  lines  using  a  selector  magnet.  The  existence  of  two  beam 
lines  gives  the  facility  increased  flexibility  in  experiment  set-up  and  operation.  The 
first  beam  line  is  called  the  Low-Energy  Neutron  Source  (LENS)  beam  line,  while  the 
second  is  called  the  Neutron  Radiation  Effects  Program  (NREP)  beam  line.  Neutron 
production  at  the  LENS  target  has  been  analyzed  extensively  [46],  however  the  mod¬ 
erator  conhguration  makes  the  LENS  target  impractical  for  use  in  radiation  effects 
studies. 
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The  NREP  beam  line  terminates  on  a  Target-Moderator-Reflector  (TMR)  stack, 
which  contains  a  beryllium  target  assembly,  space  for  a  moderator  containment  ves¬ 
sel, and  external  shielding.  The  shielding  primarily  consists  of  borated  polyethylene 
bricks  to  absorb  neutrons  and  lead  to  absorb  gamma  radiation,  with  the  goal  of 
limiting  radiation  exposure  outside  the  TMR. 

Spectral  Collapsing. 

Since  facilities  produce  neutrons  in  various  ways,  every  facility  is  unique  in  its 
own  neutron  flux  spectral  shape.  In  order  to  compare  results  of  experiments  between 
facilities,  a  concerted  effort  was  made  to  determine  a  standard  method  to  collapse 
any  spectral  shape  into  an  equivalent  number  of  monoenergetic  neutrons  by  means 
of  a  displacement  damage  kerma.  Messenger  [47]  was  the  first  to  suggest  a  curve  fit 
for  experimental  data  from  Smits  and  Stein  [48].  The  equation  was 

D{E)  =  AE[l-  ,  (7) 

where  D(E)  is  the  lifetime  damage  constant,  while  A  and  B  were  fitting  coefficients. 
From  there,  multiple  researchers  [49-53]  have  refined  the  curve. 

Currently,  the  American  Society  for  Testing  and  Materials  (ASTM)  publishes 
most  of  the  community-developed  standards  for  radiation  effects  testing,  to  include 
the  standard  by  which  a  neutron  spectrum  can  be  collapsed.  The  standard  num¬ 
ber  is  E722-09el,  and  the  description  is  Standard  Practice  for  Characterizing  Neu¬ 
tron  Fluence  Spectra  in  Terms  of  an  Equivalent  Monoenergetic  Neutron  Fluence  for 
Radiation- Hardness  Testing  of  Electronics.  It  uses  the  damage  kerma  folded  into  the 
incident  spectrum  to  determine  equivalence. 

The  equivalent  monoenergetic  fluence,  ^eq,Eref,mat  is  calculated  as  in  Equation 
Set  8. 
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/  ^E)Fo,mat{E)dE 

d^eq,E^f,f,mat  ^  (8) 

E  D,E^f,f , mat 

$(£')  =  incident  neutron  spectrum,  energy  dependent 
ED,mat{,E)  =  displacement  damage  function  for  material 
being  irradiated,  also  energy  dependent 
ED,Er^f,mat  =  reference  value  of  damage  function  for 
equivalent  energy,  Eref 

The  radiation  effects  community  typically  uses  1  MeV  as  the  reference  standard,  and 
the  accepted  reference  value  of  the  damage  function  at  this  energy  is  95  MeV-mb. 
Therefore,  the  1-MeV  Equivalent  (Si)  neutron  flux  can  be  calculated  as 

OO 

f  <S>(E)Fo,s,{E)dE 

=  "  gi[MeV-  mb] 

Additionally,  the  neutron  energy  spectrum  hardness  parameter  for  the  irradiated 
material,  Hmao  is  dehned  as  the  ratio  of  the  equivalent  monoenergetic  neutron  fluence 
to  the  true  total  fluence  [54],  or 


E-mat 


$ 


eq^Eref  ^mat 

f  ^(E)dE 


(10) 


0 

This  equivalence  principle  will  be  investigated  during  this  research,  and  the  results 
will  be  presented  later  in  this  document. 
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III.  First-Principles  Modeling  of  BJT  Peformance 


In  order  to  evaluate  the  relationship  between  short-pulse  irradiation  and  model 
BJT  gain  degradation,  a  physics-based  model  with  time-dependent  formation  and 
annealing  was  developed.  In  this  section,  the  models  used  to  describe  gain  degrada¬ 
tion  in  BJT  devices  will  be  presented,  the  basics  of  a  hrst-principles  model  will  be 
developed,  and  preliminary  validation  and  verihcation  will  be  described. 

3.1  Historical  Modeling  Efforts 

The  Non-Ionizing  Energy  Loss  (NIEL),  is  a  hgure  of  merit  that  describes  the  dam¬ 
age  potential  for  displacement  by  particle  radiation  in  solids.  NIEL  can  be  correlated 
between  particles  (protons,  neutrons,  etc.)  at  any  energy  via 

Na 

NIEL  =  {ae%  +  a,T,) ,  (11) 

where  Nav  is  Avagadro’s  number  and  A  is  the  atomic  weight.  In  Equation  11,  a  is 
the  cross  section  while  Ti  and  are  the  effective  average  recoil  energies  (corrected 
for  ionization  loss)  for  elastic,  e,  and  inelastic,  i,  scatterings. 

One  of  the  main  motives  for  use  of  NIEL  in  radiation  effects  experiments,  as 
Srour  et  ah  state  [1],  is  “that  the  signihcant  nuclear  weapons  effects  database  could 
be  mined  to  predict  a  device  response” .  The  copious  experiments  of  the  nuclear  test 
era  provided  a  treasure  trove  of  data  from  measurements  taken  across  a  broad  range 
of  nuclear  weapon  types  and  yields.  The  various  radiation  environments  explored 
allowed  researchers  to  empirically  correlate  device  responses  and  test  their  theories 
scaled  down  in  the  lab. 

In  these  experiments.  Equation  11  and  the  Messenger-Spratt  Equation  [2]  were 
applied  to  bipolar  technology  since  current  gain  scaled  linearly  with  fluence.  The 
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Messenger-Spratt  Equation  describes  semiconductor  material  parameter  response  to 
incoming  radiation.  In  this  case,  it  was  found  that  device  current  gain  after  irradia¬ 
tion,  can  be  modeled  as 


wri^KiEME) 


(12) 


where  /3o  is  the  initial  device  gain  before  irradiation,  K{E)  is  the  particle-  and  energy- 
dependent  damage  factor,  and  <h(E)  is  the  incident  particle  fluence  in  n/crri^. 

The  Messenger-Spratt  equation  was  proven  by  Messenger  [3]  in  1973  by  using  a 
series  expansion  of  the  base  transport  factor,  ctt,  per  the  work  done  by  Cover  [4], 
The  expansion  is 


ttT  =  1  — 


(13) 


where  W  is  the  depletion  width,  L  =  ^ Dt  is  the  diffusion  length,  and  7/  is  a  a 
coefficient  determined  by  successive  iteration  of  the  equations  for  excess  minority 
carrier  concentration  in  the  base  and  the  minority  current  density  (please  see  ??  for 
more  detail).  For  the  diffusion  length,  D  is  the  diffusion  constant  and  r  is  the  minority 
carrier  lifetime.  Messenger  made  the  assumption  that,  since  this  series  is  rapidly 
convergent,  it  can  be  truncated  after  the  Ui  term.  Thus,  from  Equation  13  and  a 
further  assumption  that  base  transport  is  the  dominant  factor  in  gain  degradation. 


1  _ 

P  ~  Dt  ’ 

Cover  also  showed  that,  if  n/r  is  described  by  a  single  time  constant,  u  = 
where  u  is  27r  times  the  common  emitter  gain  bandwidth  product.  From  Messenger 
and  Spratt’s  original  discussion  of  carrier  lifetime  degradation. 


1  _  1  <!> 
T  Ti^  K' 


(15) 
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where  K  is  the  carrier  lifetime  degradation  coefficient.  If  Equations  14  and  15  are 
combined  with  the  description  of  u,  the  Messenger- Spratt  equation  is 


/d/  /do  OOK 


(16) 


which  matches  with  Equation  12  if  it  is  noted  that  the  damage  constant,  K{E),  and 
{uK)~^  are  interchangeable  since  K  is  just  a  coefficient. 

There  have  been  numerous  studies  done  over  the  intervening  years  that  use  Equa¬ 
tion  12  to  predict  BJT  performance.  From  the  early  years  of  Frank  et  al.  [5]  and  Ram¬ 
sey  et  al.  [6]  to  the  more  recent  Williams  [7],  Bielejec  [8],  Li  [9],  and  Consolandi  [10], 
the  model  of  damage  response  in  a  silicon  BJT  has  remained  relatively  unchanged.  It 
is  an  adequate  measure  of  BJT  gain  degradation  in  certain  situations.  However,  it  is 
extremely  important  to  realize  that  there  are  numerous  assumptions  that  are  made. 

The  Messenger-Spratt  equation  assumes  that  the  hnal  gain,  /d/  is  measured  a 
long  time  after  irradiation.  In  essence,  the  damage  constant  K  only  deals  with  the 
permanent  defects  left  in  the  device  after  the  short-term  defects  have  annealed  away. 

The  most  important  assumption  used  in  the  Messenger-Spratt  equation  is  the  lack 
of  consideration  of  time.  Since  its  development  in  the  1950’s,  the  Messenger-Spratt 
equation  has  been  based  on  the  assumption  that  irradiation  times  for  displacement 
damage  experiments  are  very  short. 

The  SPR-III  was  an  often  used  and  well-characterized  irradiation  facility,  and 
it  had  two  modes  of  operation:  the  more  commonly  used  pulse  mode  and  the  less 
commonly  used  steady  state  mode.  The  steady  state  mode  was  less  commonly  used 
because  it  induced  a  large  build-up  of  hssion  fragment  inventory  in  the  core  [11].  For 
pulsed  operation,  bare  uranium  metal  plates  were  brought  together  and  the  heat  of 
hssion  caused  the  plates  to  expand  far  enough  apart  that  it  shut  down  the  reaction 
in  less  than  1  ms.  Neutron  huences,  in  terms  of  1-MeVequivalent  (Si)  neutrons,  in 
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the  “glory  hole”  of  the  core  were  originally  on  the  order  of  10^^ — 2"'  ,  .  Through  a 
number  of  upgrades  and  modihcations  to  the  system,  the  pulse  width  increased  from 
50  ps  for  SPR-II  [12],  to  approximately  100  ps  for  SPR-III  [13]. 

The  emphasis  on  using  a  pulsed  neutron  source  is  that,  if  neutrons  interact  over 
short  time  intervals,  the  defects  created  in  the  silicon  bulk  in  BJT  devices  do  not  have 
time  to  migrate  and  form  additional  complexes  before  the  pulse  is  complete.  Most  of 
the  point  defects,  vacancies  and  interstitials,  have  recombined  in  the  first  10  ns  [14]; 
however,  short  term  migration  and  rearrangement  of  the  permanent  defects  can  take 
10-2  to  10^  s. 

Another  point  to  note  about  the  application  of  the  Messenger- Spratt  equation  is 
that  it  contains  theory  that  is  based  solely  on  experimental  evidence.  Short  pulse, 
low-flux  experiments  would  underestimate  the  combined  effects  of  co-located  damage 
sites  and  long  interaction  times  would  overestimate.  The  expression  was  developed 
by  assuming  Shockley- Read-Hall  recombination  [15],  while  the  defects  introduced 
by  irradiation  were  considered  as  recombination  centers.  At  no  point  is  time  or 
defect  density  considered,  thus  it  is  intrinsically  time-independent.  It  is  an  adequate 
description  for  most  of  the  analysis  accomplished  over  the  last  60  years  because  it 
was  used  on  single-pulse  neutron  irradiation  facilities  and  reactors. 

3.2  Model  Development  Progression 

A  mathematical  model  was  sought  that  would  predict  the  degradation  of  BJT 
performance  parameters  with  neutron  irradiation.  This  model  had  a  number  of  steps 
in  its  progression,  which  will  now  be  described  in  greater  detail. 
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Step  One  —  General  Defect  Mode. 


The  Messenger-Spratt  relationship  for  gain  is  given  in  Equation  12.  Equation  12 
could  be  rewritten  as 

n 

(17) 


2=1 

•th 


1  1 

13  f  (3  o' 


where  Kd^i  and  Nd^  are  the  damage  constant  and  defect  population.  That  is 
to  say  that  the  damage  in  the  device  is  directly  proportional  to  the  defect  species 
concentration.  This  means  that  the  best  way  to  predict  the  gain  response  as  a  function 
of  a  time-dependent  flux  of  neutrons  such  as  the  one  found  at  the  lU  NREF  is  to 
properly  model  the  defects  as  they  are  created  and  anneal.  As  a  first  order  approach, 
a  Bateman-like  [16]  equation  can  be  developed  that  describes  the  time  rate  of  change 
of  the  number  of  defects  in  the  device,  Nd,  as 


— =  Created  —  Destroyed  =  ad(t> - Nd, 

at  Td 


(18) 


where  ad  is  the  cross  section  for  formation  of  a  defect,  0  is  the  defect-forming  particle 
flux,  and  Td  is  the  defect  lifetime.  This  is  a  relatively  simple  first-order,  linear,  ordinary 
differential  equation  that  has  an  exact  solution,  which  is 


(19) 


where  ^  is  the  constant  of  integration.  If  the  initial  population  is  Nd  =  0,  then 
Equation  19  becomes, 

Nd{t)  =  (1  -  .  (20) 

A  plot  of  this  result  is  displayed  in  Figure  15.  As  is  expected.  Equation  20  depicts 
growth  in  the  population  of  defects  that  approaches  an  asymptote  at  ad(f>Td  with  a 
characteristic  time  constant.  However,  if  the  flux  were  turned  off  (he.  0  — )■  0),  the 
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Figure  15.  Using  a  single  defect  species  results  in  the  growth  of  the  defect  population 
that  approaches  an  asymptote  with  a  time  constant  characteristic  of  the  relaxation  of 
defect  species  during  irradiation. 

model  predicts  an  exponential  drop  in  the  defect  population  to  0.  From  literature, 
this  cannot  be  the  case  as  some  displacement  damage  is  residual  in  the  material. 
Therefore,  the  model  must  be  developed  further,  and  the  next  step  is  to  include 
multiple  defect  species  interactions. 


Step  Two  —  Model  Including  Multiple  Defect  Species. 

If  this  concept  is  extended  to  two  defect  types,  vacancies  and  di-vacancies,  the 
physical  description  becomes  more  complicated.  Both  defects  have  a  hnite  proba¬ 
bility,  cr^  and  of  being  created  via  the  neutron  flux.  Additionally,  both  defect 
species  may  anneal  back  to  a  “perfect”  crystal  lattice  with  some  lifetime,  and 
However,  there  is  also  a  chance  that  a  vacancy  can  become  a  di-vacancy  by  capturing 
another  vacancy.  Additionally,  there  is  the  chance  that  a  di-vacancy  can  become  a 
single  vacancy  by  splitting  or  having  one  vacancy  recombine.  This  description  can  be 
translated  into  two  differential  equations  as  Equation  Set  21. 
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In  Equation  Set  21,  0  is  the  flux  of  incident  damaging  particles  and  Tv^w  and  Tw^v 
are  the  lifetimes  of  a  vacancy  capturing  another  vacancy  to  become  a  di-vacancy  and 
of  a  di-vacancy  splitting  into  two  vacancies,  respectively. 

As  this  step  of  the  process  is  still  hypothetical,  arbitrary  values  for  lifetime  and 
cross  section  were  selected  to  illustrate  the  model  behavior.  Figure  16  displays  an 
example  of  the  two  defect  model  showing  how  the  vacancy  population  rises  quicker 
than  the  di- vacancy  population,  but  it  reaches  an  equilibrium  faster  due  to  its  shorter 
lifetime.  In  this  example,  the  parameter  values  were  set  as  shown  below: 


Tv  T,,7;  Tj;. 


0.05s, 


0.5s, 


0  =  100cm“^s“^, 
av  =  O.lcm^, 
aw  =  O.Olcm^. 


This  system  of  equations  can  be  solved  analytically;  and,  using  the  parameters 
selected  above,  the  resultant  time-dependent  equations  are  described  as  Equation 
Set  22. 


Nvit)  =  0.56  +  0.12e-^‘  -  0.686“^* 
Nw{t)  =  0.61  +  0.58e-^‘  -  l.lOe”^* 


(22) 


These  first  steps  in  the  model  development  were  purely  hypothetical,  as  it  is 
known  that  displacement  damage  in  semiconductors  cannot  be  lumped  into  one  or 
two  categories  of  defects.  Even  a  model  using  two  species  does  not  incorporate  the 
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Figure  16.  For  this  arbitrary  selection  of  parameters,  the  vacancy  population  rises 
quickly,  but  reaches  a  secular  equilibrium  at  a  lower  level  than  the  di-vacancies  during 
irradiation. 


complexity  necessary  to  account  for  many  of  the  behaviors  of  device  performance 
under  irradiation. 


Step  Three  —  Incorporation  of  Precursors. 


Next,  the  model  was  further  developed  to  include  precursors.  In  any  crystalline 
matrix,  there  are  imperfections  and  bonding  sites  where  a  silicon  atom  or  atoms  are 
still  a  part  of  the  lattice  structure  but  are  not  properly  aligned.  This  can  essentially 
be  accounted  for  by  considering  these  sites  as  a  separate  species  with  a  cross  section 
for  conversion  into  vacancy  or  divacancy,  asv^v  and  asi'^w  respectively.  However,  if 
it  is  assumed  that  these  precursors  cannot  be  regenerated,  then  it  is  represented  by 
an  extension  of  Equation  21  as 


r 


dNy 

dt 


Ny  Nyy 

'^v  '^v — Vvv 


/  dNyy 
1  dt 


^vv  I 

'^vv  '^vv — >v 


dt 


“1“  Co, 


(23) 
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where  Cq  is  a  generation  mechanism  for  the  silicon  precursors.  Equation  Set  23  can  be 
rearranged  to  solve  for  each  of  the  three  species  and  collected  into  a  vector  equation 
as  in  Equation  Set  24. 


N'  {t)  =  AN  {t)  +  b 


'  A'y 

'  Ay 

N'  {t)= 

A'yy 

,Af(«)  = 

Ayy 

^  N's^'  y 

^  Nsi'  j 

/ 

Oj\  Q>2 

(pas 

i'^v 

A  = 

Oj2  Oj\ 

4^^  Si'  ^vv 

^  0  0 

0  Si' Si' ^vv ^ 

(  ^  \ 

b{t)= 

^vv4^ 

V  / 

(24) 


In  Equation  Set  24,  oi  =  ^  and  02  =  ^  The  general  solution  to 

Equation  Set  24  is, 

N  (t)  =  (25) 


In  Equation  25,  Ai,  A2,  and  A3  are  the  eigenvalues  found  by  det  (A  —  AI)  =  0;  ki, 
k2,  and  k^  are  constants,  and  .^1,  .^2,  and  .^3  are  eigenvectors  corresponding  to  the 
eigenvalues.  These  constants  are  determined  using  the  initial  conditions,  which  are 


Afv  (0)  =  Nvv  (0)  =  Nvv,o,  andNsi'  (0)  =  Nsi'^o- 

Equation  25  can  be  solved  for  three  eigenvalues  using  the  solution  for  the  roots  of 
a  cubic  equation.  This  indicates  that  the  analytic  solution  of  a  model  containing  just 
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three  defect  species  was  much  more  complicated  than  the  two-defect  model.  At  this 
point,  analytic  solutions  were  abandoned  in  favor  of  numerical  solution  methods  such 
as  Runge-Kutta  Ordinary  Differential  Equation  (ODE)  solvers  implemented  using  the 
Matlab  [17]  mathematical  software  package. 

Step  Four  —  Incorporating  Concentrations,  Cdef- 

In  order  to  better  incorporate  experimental  results  into  the  model,  the  equations 
were  converted  to  concentration  variables,  denoted  with  square  brackets,  [],  instead 
of  absolute  population  variables.  This  approach  is  described  in  the  following  sub¬ 
sections. 

Change  to  Concentration  Variables. 

The  following  was  done  after  considering  the  work  of  Verner,  Gerasimenko,  and 
Corbett  [18]  in  which  they  consider  the  reaction  kinetics  of  defects  in  semiconductor 
processing.  Thus,  there  are  formation  and  dissociation  reactions,  driven  by  captnre 
and  dissociation  rates.  For  example,  a  vacancy  combining  with  another  vacancy  to 
form  a  divacancy  can  be  written  as. 


IV]  +  [V] 


ivv.]- 


The  capture  rate,  Ocap,  and  dissociation  rate,  aaiss,  are  given  by 


(26) 


^cap  DyV 

^diss  ^dissj  5 


(27) 


is  the  diffusion  coefficent  for  vacancies  (and  the  assnmption  is  made  that  the 
diffusion  rate  of  a  higher  order  vacancy  complex  is  very  small  in  comparison,  to  make 
it  negligible  and  set  equal  to  0),  reap  is  the  capture  radius  of  one  species  by  another. 
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i^diss  is  the  vibrational  frequency  of  dissociation  for  a  particular  species,  Ediss,  is  the 
energy  associated  with  that  dissociation,  ks  is  Boltmann’s  constant,  and  T  is  the 
temperature. 

Additionally,  annihilation  occurs  when  a  vacancy  [V]  combines  with  an  intersti¬ 
tial  [I]  as  [V^]  4-  [/]  — )■  Lattice.  In  this  case,  the  annihilation  rate  coefficient  would  be 
«ann  =  47r  +  Di)  Tcapy+i-  Siuce  the  interstitial  and  vacancy  are  capable  of  diffus¬ 

ing  with  comparable  speeds  throughout  the  semiconductor,  both  diffusivities  should 
be  considered  in  creating  the  rate  coefficient  of  annihilation.  This  may  not  be  the 
case  for  point  defects  interacting  with  larger  clusters  or  large  impurity  atoms  whose 
diffusivities  are  much  smaller  and,  thus,  negligible  in  comparison  with  those  of  the 
point  defects. 

Impurity  interactions  also  occur.  For  example,  a  vacancy  may  interact  with  an 

impurity,  [X],  and  this  can  be  denoted  as  [V]  -|-  [X]  i - >  [V  ■  X]  with  relevant  capture 

and  dissociation  coefficients.  The  impurity  is  often  considered  “stationary”  in  regards 
to  the  vacancy,  so  it  has  no  diffusion  coefficient.  The  same  development  can  be  applied 
to  all  possible  defect  species,  so  the  problem  has  a  large  number  of  unknowns  and 
equations.  Therefore,  the  problem  must  be  focused  on  those  that  affect  the  device 
operation  in  order  for  the  solution  to  be  tractable. 

For  this  research,  many  combinations  of  defect  species  were  considered.  Initially, 
the  list  of  species  was  limited  to  vacancies,  divacancies,  interstitials,  vacancy-impurity 
complexes,  divacancy-impurity  complexes,  and  impurities.  There  are  nine  interactions 
that  account  for  creation  or  destruction  of  these  six  species.  They  are  tabulated  in 
Table  3. 
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Table  3.  Description  of  nine  interactions  used  for  six  defect  species  model. 


Formula 

Creation 

Annealing 

Gv  — ^  [V] 

Generation  of  Vacancies 

N/A 

Gi  [-^] 

Generation  of 
Interstitials 

N/A 

[F]  +  [F]  ^  [FF] 

Gapture  of  vacancy  by 
another  vacancy  to 
become  divacancy 

Dissociation  of 
divacancy  into  two 
vacancies 

[V]  +  [/]  -^Lattice 

N/A 

Annihilation  of  vacancy 
with  interstitial  to 
become  lattice 

[F]  +  [X]  ^  [FX] 

05 

Gapture  of  vacancy  by 
impurity 

Dissociation  of 
vacancy-impurity 
complex 

[F]  +  [FX]  ^  [FFX] 

0:7 

Gapture  of  vacancy  by 
vacancy-impurity 
complex 

Dissociation  of 
divacancy-impurity 
complex 

[I]  +  [FX]  [X] 

Gapture  of  interstitial 
by  vacancy-impurity 
complex 

N/A 

[I]  +  [Fl/X]  ^  [FX] 

Gapture  of  interstitial 
by  divacancy-impurity 
complex 

N/A 

[FF]  +  [X]  ^  [FFX] 

ail 

Gapture  of  divacancy  by 
impurity 

Dissociation  of 
divacancy-impurity 
complex 
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Each  of  these  interactions  can  be  converted  into  differential  equations  and  collected 
into  a  system  of  six  ODEs.  As  an  example,  [V]  +  [E]  [VV]  models  a  reaction 

a2 

going  to  the  right  as  a  destruction  of  2  vacancies.  Therefore,  they  can  be  written  as 


dC^ 

dt 


-2aiC^C^. 


(28) 


This  reaction  also  models  the  creation  of  one  divacancy,  so. 


da, 


dt 


2aiC^a. 


(29) 


After  considering  each  reaction  and  reversibility /irreversibility  of  the  reaction,  the 
system  of  differential  equations  becomes: 


dCy 

dt 


dC^y 

dt 

da 

dt 

dCyx 

dt 

dC,,'i,x. 


dt 

dC^ 

dt 


Gy  2(y.\CyCy  T  20i2C yy  Oi^CyC^  ^  ^G yC x  T  CX^Gyx  ^qG yG yx  , 

T  CXlGyyx, 

C^lGyGy  Oi2Gyy  (^l^GyyGx  T  C^lQGyyx, 

Gi  —  a^GyGi  —  agGiGyx  —  agGiGyyx,  (30) 

Oi^GyGx  Cii^Gyx  ^qG yG yx  T  Oi'j'Gyyx  T  Oi^GiGyyx, 

^CiGyGyx  Oi'lGyyx  Oi^G j^Gyyx  ^]_oGyyx  T  ^  1 1  , 

CX^GyGx  ^sGyx  T  ^%GiGyx  ^1\G yyG x  T  Ciil()GyyX‘ 


In  this  system  of  equations,  the  concentration  notation,  [],  was  replaced  with  the  letter 
G  and  a  subscript.  For  example,  the  concentration  of  vacancies  noted  previously  as 
[V]  was  changed  to  Gy.  Additionally,  while  impurities  cannot  be  created,  complexes 
that  include  clusters  may  dissociate  into  the  constituent  parts,  leaving  the  impurity 
atoms  behind.  For  the  sake  of  the  model,  this  is  equivalent  to  “creating”  them. 
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Generation  Rate  of  Defects. 


In  the  system  of  equations,  Gy  and  G,  are  the  generation  rates  of  the  vacancies  and 
interstitials,  respectively.  These  are  functions  of  the  neutron  flux  and  lattice  proper¬ 
ties.  The  model  provides  a  method  for  generation  of  Frenkel  pairs  by  incorporating 
the  displacement  damage  kerma  along  with  the  expected  neutron  flux  spectrum.  The 
displacement  damage  kerma  in  silicon  is  tabulated  in  ASTM  Standard  E722-09,  the 
Standard  Practice  for  Characterizing  Neutron  Energy  Fluence  Spectra  in  Terms  of  an 
Eguivalent  Monoenergetic  Neutron  Fluence  for  Radiation- Hardness  Testing  of  Elec¬ 
tronics  [19],  and  is  displayed  in  Figure  17.  In  this  standard,  the  user  is  directed  to 
ICRU  reports  13  [20]  for  more  information  on  the  kerma. 


Figure  17.  The  displacement  damage  kerma  for  a  silicon  semiconductor  as  tabulated 
in  the  ASTM  Standard  E722  [19]. 


In  the  model  developed  in  this  research,  it  is  assumed  that  Frenkel  pairs  can  be 
created  when  the  interacting  particle  deposits  energy  surpassing  a  threshold  value, 
Efp.  Therefore,  the  generation  rate  of  Frenkel  pairs  is  defined  as 


mE)F(E)dE 

7-1  ? 


(31) 
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where  F{E)  is  the  displacement  kerma  for  neutrons  in  silicon.  This  results  in  a 
constant  generation  rate  in  the  bulk  silicon  that  is  directly  dependent  upon  the  flux. 

Step  Five  —  Development  of  Reaction  Rate  Parameters. 

The  system  of  non-linear  differential  equations  developed  in  Section  3.2  can  be 
solved  if  the  reaction  coefficients  are  known.  Each  reaction  has  a  reaction  rate  co¬ 
efficient  that  dictates  the  process  progression  and  can  be  listed  in  two  categories: 
capture  and  dissociation  rates.  Every  a  in  Equation  30  designates  a  rate  coefficient. 
Rates  are  most  often  often  empirically  derived  due  to  their  many  dependencies,  most 
on  material  quality.  Therefore,  they  have  large  uncertainties.  In  the  following  sub¬ 
sections,  the  diffusion  coefficients  of  vacancies  and  interstitials,  capture  radii,  and 
the  dissociation  frequency  are  investigated  in  more  detail.  These  parameters  fit  into 
model  as  a  part  of  many  rate  coefficients,  a. 

Vacancy  and  Interstitial  Diffusion  Coefficient. 

As  an  example  of  the  uncertainty  in  knowledge  of  the  rate  parameters,  Wang  et 
al.  [21]  stated  that  the  diffusion  coefficient  of  vacancies  in  Si  at  room  temperature 
is  between  3.4  x  10“®  cm^/s  and  14.0  x  10“®  cm^/s.  However,  in  1998,  Hallen  et  al. 
[22]  measured  the  vacancy  diffusivity  in  Si  at  room  temperature  as  4.2  x  10“®  cm^/s. 
Voronkov  and  Falster  [23]  have  measured  a  value  at  room  temperature  of  3.0  x  10“^° 
cm^/s  and  indicate  it  has  an  Arrhenius  behavior  with  a  pre-factor  of  2.3  x  10“^  cm^/s 
and  an  associated  energy  barrier  of  0.35  eV.  This  leads  to  a  room  temperature  value 
of  1.85  X  10“^°  cm^/s.  In  essence,  there  is  no  single  definitive  value  for  the  diffusion 
coefficient  of  vacancies  in  silicon  at  room  temperature. 

The  diffusion  coefficient  for  interstitials  is  similar.  Hallen  et  al.  [22]  measured  a 
value  of  3.2  x  10“^^  cm^/s  at  room  temperature,  whereas  Libertino  and  Coffa  [24,25], 
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measured  values  of  1.5  x  10  cm^/s  and  2.0  x  10  cm^/s.  Consistent  in  the 
literature  is  that  vacancies  are  substantially  more  mobile  than  interstitials. 

Capture  Radii. 

Less  information  exists  regarding  captnre  radii  for  the  various  species.  Cowern 
et  al.  [26]  report  effective  capture  radii  for  boron  interstitials  and  snbstitntionals  on 
the  order  of  0.25-1  nm.  This  is  reasonable,  as  the  typical  spacing  between  two  silicon 
atoms  is  abont  0.2  nm. 


Dissociation  Frequency. 


The  “vibrational  frequency”  of  dissociation  describes  the  rate  at  which  larger  de¬ 
fect  species,  snch  as  the  vacancy-impnrity  complex,  obtain  enough  vibrational  energy 
to  break  part  into  their  constitnent  parts.  In  their  seminal  1965  paper,  Watkins  and 
Corbett  [27]  state  a  valne  of  approximately  10^^  s“^  for  something  they  call  the  “fre- 
qnency  factor”  for  silicon  defect  dissociation.  With  this  exception,  there  is  no  other 
information  on  the  rate  at  which  defect  species  may  dissociate. 

Fnrnhashi  and  Taniguchi  [28]  state  that  the  activation  energy  of  dissociation  is 
defined  as  the  snm  of  the  binding  energy  of  the  complex  with  the  smaller  migration 
energies  of  the  separated  defects.  Using  values  for  a  divacancy-oxygen  complex,  one 
gets 


V2O  ^  U2  +  Oi  2.98eV  =  1.3eV  +  1.68eV, 

(32) 

no  ^  12  +  VOi  2.24eV  =  0.33eV  +  1.91eV, 

where  O*  represents  an  oxygen  interstitial.  Both  reactions  in  Eqnation  Set  32  can 
occnr,  bnt  the  second  is  more  likely  because  of  its  smaller  activation  energy.  In  this 
same  paper,  they  provide  barrier  energies  for  mnltiple  species  as  Uf,  {Vacancy)  =  0.33 
eV,  Eb  (Divacancy)  =  1.3  eV,  Eb  {VO)  =  2.0  eV,  Eb  (O*)  =  2.53  eV. 
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Step  Six  —  Gain,  (3,  Calculation. 


From  Section  2.1,  device  gain  in  a  semiconductor  is  defined  as  in  Equations  3  and 
4.  Each  parameter  in  Equations  is  device-dependent,  in  that  they  rely  upon  the 
characteristics  of  a  particular  device  that  has  been  fabricated;  and  defect-dependent, 
in  that  they  rely  on  the  amount  and  type  of  defects  in  the  material. 

Semiconductor  Parameters. 

Minority  carrier  concentrations  are  determined  from  basic  principles.  For  doped, 
or  extrinsic,  semiconductor  material,  a  region  is  n-  or  p-type  material  if  the  majority 
carrier  is  electrons  or  holes,  respectively.  This  is  done  by  adding  electron  donor 
or  acceptor  impurities.  In  an  n^pn  BJT  such  as  the  2n2222A,  the  base  is  p-type, 
while  the  emitter  and  collector  are  n-type.  Physical  characterization  of  the  2n2222A 
devices  was  not  possible  for  this  work.  However,  Keiter,  Russo,  and  Hembree  [13] 
used  spreading  resistance  measurements  as  well  as  Secondary  Ion  Mass  Spectroscopy 
(SIMS)  to  characterize  the  2n2222A  silicon  BJTs  that  they  were  testing,  and  this 
profile  is  adopted  as  a  “typical”  device  doping  profile  for  the  sake  of  analysis.  A 
replotting  of  their  data  is  displayed  in  Figure  18. 

From  a  depth  of  0-1.5  /rm,  the  phosphorous  donor  impurity  atoms  dominate  the 
boron  acceptor  impurity  atoms,  thus  the  region  is  n-type  and  represents  the  emitter 
region.  On  the  other  hand,  from  1. 5-3.0  /rm,  the  opposite  is  true,  making  the  region 
p-type  and  the  base.  Since  both  impurity  types  are  present  in  each  region,  it  is 
compensation  doped,  meaning  some  impurity  sites  are  canceled  out  by  those  of  the 
opposite  type.  Assuming  full  ionization  of  donors  and  acceptors,  the  majority  and 
minority  carrier  concentrations  in  the  base  are 
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Figure  18.  The  dopant  concentrations  (5+  and  P  )  in  the  base-emitter  region  of  a 
2n2222A  silicon  BJT  were  determined  from  spreading  resistance  measurements  and 
SIMS  data  [13], 


Pp  = 

Up  = 


Na-Nd  + 


nj 


Pp 


while,  in  the  n-type  emitter, 


nn 

Pn 


1 

2 

2 


Nd-Na  +  \j (Na  -  Nnf  +  4nf 


Ur. 


(33) 


(34) 


where  Na  and  Nr,  are  the  acceptor  and  donor  dopant  concentrations,  respectively,  and 
rii  is  the  intrinsic  carrier  concentration,  which  is  9.65  x  10®  cm“^  for  room  temperature 

silicon. 

From  Figure  18,  the  phosphorous  donor  concentration  varies  from  approximately 
10^®  cm“^  to  10^®  cm“®  across  its  depth.  However,  using  the  abrupt  junction  approx¬ 
imation  to  simplify  the  analysis,  a  constant  value  of  5  x  10^®  cm“®  is  assumed.  The 
same  approximation  is  done  for  the  base,  resulting  in  an  average  value  of  5  x  10^^  cm“® 
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in  the  region.  Since  Figure  18  does  not  explicitly  describe  the  collector  concentration, 
it  is  assumed  that  they  are  10^^  cm“^  P  donors  and  10^^  cm“^  B  acceptors.  Therefore, 
using  Equations  33  and  34,  the  emitter,  base,  and  collector  carrier  concentrations  are 

Uns  =  4.6  X  10^®  [crrT^]  =  20  , 

PpB  =  ^  ^  10^®  [cm~^]  ,  Upg  =  10^  [cm~^]  ,  (35) 

Unc  =  9  X  10^"^  ,pnc  =  10^  . 

Of  particular  note  from  these  results  is  that  the  donor  concentration  drives  the  emitter 
electron  concentration,  leaving  few  holes  in  this  region,  which  is  a  very  desirable 
condition  because  it  increases  the  emitter  efficiency. 

Reexamining  Equation  3,  the  two  important  concentration  values  for  calculating 
gain  are  Png  and  or  the  minority  carrier  concentrations  for  both  emitter  and 
base.  An  assumption  already  made  but  restated  here  is  that  the  neutron-induced 
defect  concentration  will  be  on  the  order  of  the  dopant  concentrations  of  the  emitter, 
and  thus  the  hole  concentration  in  the  emitter  will  not  be  affected  by  the  neutron 
radiation  environment.  From  the  results  displayed  in  Equation  35, 

Pngo  =  Png  it)  =  20  [cm~^]  ,  (36) 

for  all  times,  t.  On  the  other  hand,  the  electron  concentration  in  the  base  will  be 
strongly  dependent  on  the  neutron  flux,  since  the  number  of  defects  in  a  typical 
irradiation  cycle  will  approach  and  possibly  exceed  the  dopant  concentration. 

The  diffusivity  of  the  minority  carriers  in  the  emitter  and  base,  Dp^  and  Dng, 
depend  on  the  carrier  trap  density  in  these  regions.  In  the  emitter,  it  is  assumed 
that  the  radiation-induced  defect  concentration  never  approaches  that  of  the  donor 
concentration,  A^Db-  Therefore,  the  donor  concentration  can  be  used  to  calculate  the 
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diffusivity  and  mobility  of  holes  in  the  emitter  by  interpolating  from  the  empirical 
data  in  Sze  [29],  displayed  in  Figure  19.  In  the  base,  however,  the  concentration  of 
neutron-induced  trap  sites  caused  by  displacement  damage  can  be  on  the  same  level 
or  higher  than  the  number  of  acceptors.  Therefore,  the  sum  of  the  two  is  necessary 
when  calculating  the  interpolated  value  of  diffusivity. 

The  diffusion  length  is  the  average  length  that  the  particular  carrier  will  travel  in 
the  device  region,  which  is  a  function  of  the  square  root  of  the  product  of  diffusivity 
and  carrier  lifetime.  Extending  this,  the  carrier  lifetime  for  some  carrier  x  is  calculated 
using  the  mobility  of  that  carrier,  as  =  {Hx  ■  ^n)  /q  where  is  the  rest  mass 
of  the  electron  and  q  is  the  fundamental  charge  constant,  set  at  1.602  x  10“^®  C.  This 
particular  parameter  will  be  affected,  albeit  indirectly,  by  neutron-induced  trap  sites 
due  to  the  change  in  mobility  and  diffusivity. 

The  change  in  the  physical  dimensions,  i.e.  Wb  and  We,  are  assumed  to  be  neg¬ 
ligible.  While  the  introduction  of  additional  trap  sites  may  alter  the  regions  of  the 
BJT,  it  will  be  small  compared  to  the  actual  dimensions  of  these  regions.  Therefore, 
they  are  set  at  a  constant  value  and  are  unchanged  during  simulation. 


Figure  19.  Diffusivity  (left)  and  mobility  (right)  as  a  function  of  impurity  concentration 
in  the  semiconductor  bulk  is  flat  till  10^®  cm“^;  it  can  be  affected  by  radiation  induced 
defects  adding  to  the  impurity  concentration. 
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Results  of  Initial  Model  Implementation. 


At  this  stage  of  development,  the  model  contained  6  coupled,  non-linear  ODEs  for 
the  interactions  of  vacancies,  divacancies,  interstitials,  impurity-vacancy  complexes, 
impurity-divacancy  complexes,  and  impurities.  A  test  of  the  model  using  hypothetical 
parameters  was  accomplished  to  determine  the  validity  of  the  results,  and  this  is 
plotted  in  Figure  21. 


Figure  20.  The  neutron  flux  for  the  hypothetical  situation  in  Figure  21  corresponds  to 
a  1-MeV  equivalent  (Si)  flux  of  5  x  10^°  n/cm^/s  using  the  ASTM  standard  method  of 
calculation. 

The  spectrum  that  was  used  for  generation  of  defects  was  the  time-averaged  neu¬ 
tron  spectrum  from  the  lU  NREF  beamline.  This  was  determined  using  MCNP-X 
and  verified  using  multiple  foil  activation  analysis.  It  is  plotted  in  Figure  20.  The 
initial  concentration  of  all  species  except  impurities  was  set  to  10^°  cm“^  to  simulate 
a  non-zero  but  negligible  amount  in  relation  to  the  number  density  of  silicon  atoms. 
The  impurity  concentration  was  the  one  exception,  as  it  was  set  to  5  x  10^®  at/cm^ 
to  correspond  to  the  typical  dopant  concentration  in  the  2n2222A  base  region.  The 
vacancy,  interstitial,  and  divacancy  diffusivities  were  set  to  3  x  10“^°  cm^/s,  2  x  10“^^ 
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cm^/s,  and  1  x  10“^®  cm^/s,  respectively.  These  values  were  selected  because  they 
fall  within  the  typical  range  found  in  literature.  The  vibrational  frequency  of  disso¬ 
ciation  was  also  set  to  a  value  found  in  literature  of  10^^  s“^  for  all  defect  species 
capable  of  dissociation  (divacancies,  vacancy-impurities,  and  divacancy-impurities). 
The  capture  radii  were  set  to  1  A  for  all  species. 


Figure  21.  The  first  simulation  using  the  model  with  six  defect  species  and  a  time- 
averaged  1-MeV  equivalent  neutron  fiux  of  5  x  10^°  n/cm~^/s  and  an  initial  impurity 
concentration  of  5  x  10^®  at/cm^. 

Figure  21  presents  the  result  of  species  formation  during  a  simulated  constant 
neutron  flux.  Most  species  initially  rise  in  concentration,  but  there  is  a  decline  of 
vacancies  after  peaking  at  about  <l>  =  3  x  10^°  cm“^.  As  in  the  previous  model 
development  steps,  most  of  the  defect  species  rise  asymptotically  as  they  approach 
the  secular  equilibrium  between  their  respective  creation  and  destruction  mechanisms. 
However,  vacancies  and  impurities  are  the  exception.  Vacancies  rise  asymptotically, 
but  begin  to  decline  after  a  time.  This  is  due  to  an  increase  in  the  concentrations 
of  defect  species  to  which  they  can  combine,  such  as  other  vacancies,  interstitials, 
or  other  higher-order  complexes.  The  impurity  concentration  declines  because  there 
are  many  avenues  of  removal  with  very  few  “creation”  mechanisms  for  the  impurity. 
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and  those  few  creation  paths  have  low  probabilities  of  occurrence  with  the  current 
parameter  values  used  in  the  model. 

Utilizing  the  results  from  this  simulation  and  applying  Equations  with  the  pa¬ 
rameters  already  described,  the  device  gain  can  be  calculated.  The  result  of  this  is 
plotted  in  Figure  22.  The  device  gain  response  should  be  linear  with  respect  to  the 
fluence  if  the  Messenger-Spratt  equation  is  applicable.  In  examining  Figure  22,  the 
relation  is  initially  linear  but  becomes  sub-linear  at  approximately  1  x  10^'^  cm“^  and 
peaks  around  3  x  10^°  cm“^.  This  indicates  that  the  modeled  device  anneals  while 
still  being  irradiated  due  to  the  increase  in  vacancies  and  recombining  events.  It  is 
because  of  this  behavior  that  the  model  was  further  amended. 


Figure  22.  Results  of  the  calculation  of  gain  change  as  a  function  of  fluence  for  the 
hypothetical  situation  uncover  an  issue  with  the  model  as  it  shows  that  the  change  in 
gain  would  begin  to  recover  during  irradiation,  which  is  not  observed  in  literature. 


Step  Seven  —  Adding  Additional  Species. 

In  the  previous  step,  the  high  diffusivity  of  vacancies  coupled  with  the  additional 
destruction  paths  caused  the  vacancy  concentration  to  decrease  during  irradiation. 
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This  behavior  lead  to  A  (1//3)  getting  smaller  while  the  device  was  being  irradiated. 
This  phenomenon  has  not  been  observed  in  literature. 

Chakravarthi  and  Dunham  discuss  the  fact  that  a  “high  supersaturation  of  va¬ 
cancies. ..  can  lead  to  their  aggregation  into  vacancy  clusters”  [30].  One  aspect  of 
neutron  displacement  damage  that  was  not  addressed  in  the  initial  model  was  the 
formation  of  vacancy  clusters.  These  clusters  form  when  one  or  more  vacancies  are 
in  close  proximity  to  other  vacancies  without  filling  adjacent  lattice  sites.  By  occu¬ 
pying  close  but  not  adjacent  sites,  they  retain  their  original  properties  and  do  not 
become  a  divacancy,  which  occupies  a  different  location  within  the  silicon  bandgap. 
However,  because  the  silicon  lattice  is  perturbed,  vacancy  mobility  is  decreased,  and 
the  vacancies  do  not  anneal  as  quickly. 

Therefore,  two  new  species  were  created  in  the  model  representing  called  2-  and  3- 
vacancy  dusters.  These  clusters  were  added  to  the  model  in  two  ways.  First,  clusters 
were  able  to  grow  or  dissolve  with  the  collection  or  expulsion  of  a  vacancy.  Second, 
the  clusters  could  interact  with  interstitials.  These  two  processes  are  described  by 
Equation  Set  37. 


K_1  +  V^Vn, 
Vn  +  I  ^  Vn-l- 


(37) 


The  rate  of  increase  of  a  species  is  given  by  aincrease  =  where  Dy  is  the 

vacancy  diffusion  coefficient  and  dZ^ture  is  the  capture  radius  of  the  n-1  cluster.  On 
the  other  hand,  the  rate  of  decrease  for  a  cluster  is  given  by  Equation  38, 


1  „  f  1 


(38) 


where  Ei,  is  the  barrier  energy  of  breaking  apart  the  cluster,  is  the  Boltzmann 
constant,  and  T  is  the  temperature.  The  pre-factor  of  (5  x  10^^)  ^  was  determined 
from  atomistic  calculations  performed  by  Bongiornio  et  ah  [31].  One  could  add  as 
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many  of  these  clusters  as  is  feasible  for  time  and  computational  resources.  However, 
for  this  analysis,  only  the  2-  and  3-vacancy  clusters  were  considered.  Taking  the  new 
processes  into  account  yields  two  new  equations. 


dt 

dC‘i„ 

dt 


ai2CjjCy  —  ai3C2v  —  cyiiCiC2v  +  +  ai^CiC-^^  —  oinC2v) 

(^isC^C2v  —  +  cx17C2v  —  OimCiC^v- 


(39) 


However,  this  does  not  fully  address  the  formation  of  defect  clusters,  since  the 
model  only  allows  them  to  be  formed  by  vacancies  being  captured  by  other  vacancies. 
However,  this  method  was  investigated  to  analyze  the  effects  of  clusters  on  the  model. 

The  2-  and  3-vacancy  cluster  calculations  lead  to  8  total  species  being  addressed. 
In  order  to  verify  the  completeness  of  the  system  of  differential  equations,  a  Q'^'^species 
was  added:  the  annihilation  species.  This  additional  species  accounts  for  the  con¬ 
centration  of  annihilations  that  occur  between  vacancies  and  interstitials  in  the  bulk. 
Since  the  total  number  of  vacancies  that  are  created  in  the  material  can  be  directly 
tracked  as  a  function  of  time,  it  can  be  compared  against  the  number  of  vacancies 
that  currently  reside  in  the  material  at  any  time  to  ensure  conservation  of  defects  in 
the  system  of  equations.  Since  this  species  is  only  included  for  accounting  purposes, 
it  can  only  be  created.  Therefore, 


dC„ 


dt 


—  agCiC^x  +  CXgCiCyxx  +  0'i4CiC2v  +  CXigCiC^^. 


(40) 
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Now,  if  Equations  39  and  40  are  added  to  Equation  30,  the  model  becomes 


dC^ 

dt 


dC^y 

dt 

da 

dt 


dCyX 

dt 

dCyyx 

dt 

dCx 


dt 

dC2v 

dt 

dC^y 

dt 

dCann 

dt 


Gy  20ilGyGy  “t“  2Q:2Gyy  Oi^CyC^  ^  ^G yC x  ^^GyX^ 

—  a^GyGyx  +  CXjGyyx  “  2Q;i2C'^C't,  +  2Q;13C'2«  +  Q;i4C'iC'2«, 

—  ai^Gy  —  ai8GyG2v  +  aigG^y, 

C^lGyGy  Oi2Gyy  OinGyyGx  C^loGyyx-, 

Gi  CX^GyGi  CX^GiGyx  ^S^G {G yyx  ^lAGiG2y  ^l^Gy^ 

—  ai^GiG^v  +  Oil7G2y, 

Cy^GyGx  ^bGyX  ^qG yG yx  Oi^Gyyx  “t“  Oi^G j^Gyyx  CX^G iG yx 

CX^GyGyx  CX^Gyyx  CXgGiGyyx  Oi^oGyyx  O^nGyyGx^ 

(X4GyGx  4“  Oi^Gyx  4“  ^sGiGyx  CXwGyyGx  4“  Ot-loGyyx^ 
oiuGyGy  —  ai8G2y  —  ai4GiG2y  4-  cti^Gy  4-  aieGiGsy  —  ai7G2y, 
cxi8GyG2y  —  4-  q;17C'2^  —  ai^GiG^y, 

a^GyGi  +  agGiGyx  +  agGiGyyx  4-  Q;i4C'iC'2v  4-  aigGiG^y. 


(41) 


With  the  model  in  this  form,  the  conservation  of  vacancies  can  be  verified  if 


Gy^O  4"  G(t)  ■  t  —  Gy(t)  4"  2Gyy(t)  4"  Gyxif)  4"  2Gyyx(t)  4"  2G2y(t)  4“  3C*3„(t),  (42) 


which  shows  that  the  sum  of  the  initial  concentration  and  the  total  generation  of 
vacancies  at  any  time,  t,  is  equal  to  the  sum  of  all  the  concentrations  of  species  with 
vacancies,  subject  to  the  number  of  vacancies  within  each  one.  Stated  in  another 
fashion,  if 


E ClCaccounted  f  or  E dCcreated 

E  dCcreated 


<  tolerance. 


then  conservation  criteria  have  been  satisfactorily  met.  Figure  23  is  a  plot  of  this 


Figure  23.  The  calculation  of  residuals  between  vacancies  “created”  and  vacancies 
“accounted  for”  is  below  the  tolerance  set  for  noise  in  the  calculation,  meaning  that 
the  value  is  effectively  zero  at  all  x-axis  values. 

calculation,  showing  that  the  residual  stays  under  8  x  10“^^  at  all  times,  which  is 
below  the  tolerance  set  for  machine  precision  in  the  numerical  calculation  of  the 
residual  itself. 

Having  ensured  conservation  criteria,  there  is  one  last  couplet  of  species  that 
was  modeled.  Originally,  the  carrier  concentrations  within  the  base  were  assumed 
constant.  This  is  not  realistic,  as  the  radiation-induced  defects  very  strongly  affect 
the  number  of  charge  carriers.  Therefore,  two  additional  species  were  modeled:  the 
holes  (majority  carrier)  and  electrons  (minority  carrier)  in  the  base. 

Typically,  divacancy  defect  species  he  at  the  approximate  center  of  the  bandgap 
in  silicon  semiconductors.  Therefore,  thermal  generation  of  carriers  occurs  from  these 
trap  sites  adding  electrons  and  holes.  Additionally,  divacancies  act  as  recombination 
centers  for  the  minority  carrier,  removing  electrons  due  to  recombination.  These 
phenomena  are  implemented  in  the  following  equations, 
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(43) 


dCh+ 

dt 

dCe- 

dt 


+  (y.2oC.u 


—  Oi2lGvvGe-  +  a22Gv 


where  G^  is  the  generation  rate  of  electron-hole  pairs  due  to  deposition  of  ionization 
energy  from  gamma  interactions  in  the  silicon  semiconductor.  In  this  equation, 


^20  =  VthCTp+Uiexp  {{Et  -  Ei)  /  (/cbT)}  , 

a2i  =  Vth(Jn- ,  (44) 

^22  =  Vth(Jn-niexp  {{Ei  -  Et)  /  (/cbT)}  . 

where  Vth  is  the  thermal  velocity  of  the  electron,  Et  is  the  mid-gap  energy,  and  Et  is 
the  intrinsic  Fermi  level  within  the  bandgap.  The  thermal  velocity  is  calculated  as 
Vth  =  \/^kBT /rrin  =  2.29  x  lO^cm/s.  For  silicon,  Et  =  1.1/2  =  0.55eV.  The  intrinsic 
Fermi  level  can  be  calculated  as 

=  ^  +  (45) 

where  =  2.66  x  and  =  2.86  x  are  the  density  of  states  in 

the  valence  and  conduction  bands,  respectively,  for  room  temperature  silicon. 

Finally,  Equation  43  can  be  combined  with  Equation  41  to  reach  the  model  state 
with  eleven  tracked  species.  This  is  given  in  Equation  Set  46. 
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dC^ 

dt 


dCvv 

dt 

dCi 

dt 


dCyx 

dt 

rir' 

^^vvx 

dt 

dCr 


dt 

dC2v 

dt 

dCsv 

dt 

rir' 

dt 

dCh+ 

dt 

dCe- 

dt 


Gy  ~i“  ^^GyCi  CX^CyCx  Oi^Gyx 

—  a^CyCyx  +  CXjCyyx  “  2ai2CyCy  +  2aisC2v  +  0'l4CiC2y 

—  ai^Cy  —  ai^CyC2y  +  aigC^y 
Oi^CyCy  Oi2Gyy  Oi\\CyyCx  Oi^gGyyx 

Gi  CX'^GyGi  (y.%G iG yx  ^qG ^G yyx  OC\4GiG2v  “1“  Oil^Gy 

—  aigGiG^y  +  anG2v 

CX^GyGx  ^^Gyx  CXqG yG yx  H“  ^jGyyX  “t“  OigG {G yyx  ^^G {G yx 

^aGyGyx  ^jGyyx  ^qG iG yyx  0^4QGyyx  Oi^GyyGx 

Oi^GyG X  ^^Gyx  C^sGlGyx  Oi\\GyyGx  Oi^gGyyx 

duGyGy  —  0'i3G2y  —  cxi4GiG2y  +  oiibGy  +  oii^GiG^y  —  cxnG2v 
Oii%GyG2y  —  aigG^y  +  al7G2y  —  ai^G^G^y 
a^GyGi  +  agGiGyx  +  agGiGyyx  +  ai4GiG2y  +  cxigGiG^y 
G^  +  a2oGyy 

G^  —  0'2lGyyGe-  +  0'22Gyy 


(46) 


Step  Eight  —  Temporal  Profile  Capability. 

The  final  step  in  the  model  development  involved  adding  the  ability  to  simulate  a 
continuously  pulsing  radiation  environment.  While  the  model  could  address  constant 
fluence  sources  such  as  a  nuclear  reactor,  accelerator-driven  neutron  sources  such  as 
the  Indiana  University  (lU)  Neutron  Radiation  Effects  Facility  (NREF)  operate  in 
pulsed  modes.  During  each  cycle,  there  is  a  period  of  beam  “on”  followed  by  a  period 
of  beam  “off” .  Therefore,  the  model  must  be  able  to  account  for  this. 
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An  approach  was  taken  by  breaking  the  calculation  steps  into  individual  cycles  of 
operative  and  inoperative  mode.  During  each  irradiation  cycle,  the  generation  rates 
are  set  to  positive  values  and  calculation  progresses  for  the  period  of  time  correlating 
to  the  pulse  width  of  the  proton  source.  Then,  when  switching  to  the  “off”  cycle, 
the  initial  conditions  of  the  simulation  are  set  to  the  hnal  conditions  of  the  previous 
cycle,  the  generation  rates  are  set  to  zero,  and  the  calculation  again  progresses  for 
the  period  of  time  that  correlates  to  the  remainder  of  a  single  pulse  period.  As  an 
example,  if  the  beam  frequency  is  20  Hz  and  the  pulse  width  is  400  ps,  this  means 
that  the  operative  mode  would  be  400  /is  while  the  “off”  state  would  last  49.6  ms. 

3.3  Initial  Model  Results 

An  initial  set  of  model  simulations  was  accomplished  to  assess  the  performance 
and  computational  intensity  of  the  model.  The  parameter  set  for  this  analysis  are 
tabulated  in  Tabled. 

The  simulation  was  a  10  s  irradiation  by  a  neutron  source  with  energy  and  flux  sim¬ 
ilar  to  the  lU  NREF  accelerator-driven  source  followed  by  a  10  s  annealing  period.  The 
spectrum  was  as  shown  in  Figure  20,  and  it  has  a  time-averaged  1-MeV  equivalent  (Si) 
neutron  flux  of  approximately  5  x  10^°  n/cm^/s. 

The  simulation  was  then  conducted  for  irradiation  and  annealing  times  from  10 
through  1000  s,  in  order  to  monitor  model  computation  times  and  ensure  that  the 
simulation  was  stable  and  repeatable.  Figure  24  displays  the  results  of  these  simula¬ 
tions. 

As  expected,  the  required  computational  time  was  roughly  linear  with  respect  to 
the  simulated  irradiation  time.  The  data  can  be  £t  to  a  plane  by  tcomp  =  O.Tfj^r  + 
0.128tanneai  +  13.36.  If  this  result  can  be  extrapolated  to  longer  irradiation  and  an¬ 
nealing  intervals,  a  20  minute  irradiation  followed  by  a  20  minute  annealing  interval 
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Table  4.  Model  parameters  selected  for  test  run  of  final  model  design.  The  lattice 
constant  elo  for  silicon  at  room  temperature  is  5.4307  x  10“®  cm. 


Parameter 

Species 

Value 

Units 

D 

Vacancy 

Interstitial 

Divacancy 

10-10 

10-15 

cm^/s 

^barrier 

Divacancy 

Vac-Imp 

Divac-Imp  to  Vac/ Vac-Imp 
Divac-Imp  to  Divac/Imp 

1.6 

1.7 

2.24 

2.98 

eV 

^diss 

Divacancy 

Vac-Imp 

Divac-Imp 

o  o  o 

CO  CO  CO 

1/s 

f  cap 

V  and  I 

V  and  V 

V  and  Imp 

V  and  Vac-Imp 

I  and  Vac-Imp 

I  and  Divac-Imp 

Divac  and  Imp 

2oo 

2oo 

2cio 

lOtto 

lOtto 

lOtto 

lOtto 

cm 

Init.Cond. 

Vacancies 

Divacancies 

Interstitials 

Vac-Imp 

Divac-Imp 

Impurities 

2-  Vac  Clusters 

3-  Vac  Clusters 
Annihilations 

Electrons 

Holes 

1013 

1013 

1013 

lOlo 

lOlo 

5  X  1013 

lOlo 

lOlo 

lOlo 

lOi 

Cimp 

cm“3 

^fp 

N/A 

40 

eV 

Ibeam 

N/A 

20 

mA 

fbeam 

N/A 

20 

Hz 

P^^beam 

N/A 

400 

fJ,S 
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Figure  24.  The  test  of  computational  time  versus  irradiation  and  annealing  time  showed 
a  nearly  linear  relationship  from  0  s  to  1200  s. 

would  take  almost  17  minutes,  and  an  hour  long  irradiation  followed  by  a  3  hour 
long  annealing  interval  would  require  just  over  an  hour  of  computational  time.  Thus, 
simulating  short  irradiation  periods  is  feasible  but  longer  irradiations  become  time¬ 
intensive  and  may  require  code  modihcations  to  accelerate  convergence. 

Also,  for  longer  irradiation  times,  the  longer  it  takes  to  fully  anneal  defects,  which 
increases  overall  computation  time.  This  was  also  expected  due  to  the  way  that  an¬ 
nealing  time  was  addressed  in  the  code.  Since  there  was  no  pulsing  during  annealing, 
the  entire  annealing  interval  could  be  approached  as  a  single  solution  step.  This  is 
in  contrast  to  the  irradiation  time,  where  each  interval  of  irradiation  and  annealing 
must  be  calculated  separately.  In  the  end,  the  model  was  found  to  be  stable  and  the 
computational  times  reasonable  for  the  current  application. 

After  testing  for  stability  and  repeatability,  one  of  the  simulation  sets  was  ana¬ 
lyzed  for  important  features.  The  10  s  irradiation  and  annealing  interval  simulation 
is  plotted  in  Figure  25.  The  hrst  feature  to  note  is  that  the  vacancy  concentration, 
which  starts  at  10^^  cm“^  rises  quickly  to  10^®  cm“^,  and  there  is  a  very  distinct 


74 


0  2  4  6  8  10  12  14  16  18  20 

time  [s] 


Figure  25.  The  plot  of  species  concentrations  after  simulating  a  10  s  irradiation  period 
followed  by  a  10  s  annealing  period  displays  a  noticeable  vacancy  defect  concentration 
increase  and  subsequent  annealing  during  each  pulse  cycle. 

growth  and  annealing  of  the  vacancies  during  each  pulse  cycle.  This  leads  to  a  dis¬ 
tinctive  sawtooth  pattern  in  the  A(l//d)  versus  fluence  graph.  Additionally,  when 
the  irradiation  is  stopped  at  10  s,  the  vacancy  concentration  drops  quickly  as  they 
are  consumed  by  annealing.  The  other  defect  species  rise  asymptotically,  as  in  the 
previous  versions  of  the  model;  and  after  irradiation,  they  remain  constant  over  the 
time  interval  simulated. 

In  Figures  26  and  27,  gain  is  calculated  and  plotted  for  the  irradiation  and  an¬ 
nealing  phases  of  the  model.  In  Figure  26,  the  change  in  gain  with  neutron  fluence  is 
plotted  during  the  irradiation  interval.  The  slope  of  this  curve  is  2.5  x  10“^^  cm^/n 
and  is  a  direct  measurement  of  the  damage  constant  for  these  simulation  parame¬ 
ters,  due  to  the  Messenger-Spratt  equation.  This  value  of  the  damage  constant  is 
about  two  orders  of  magnitude  larger  than  the  1.5  x  10“^®  stated  in  literature  [32]  for 
2N2222  Si  BJTs  for  1-MeV  equivalent  neutrons  in  silicon.  This  may  be  indicative  of 
parameter  values  that  are  incorrect.  It  may  also  indicate  that  there  are  other  physical 
processes  at  work  that  are  not  captured  in  the  model. 
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Figure  26.  The  change  in  gain  as  a  function  of  fluence  during  irradiation  calculated 
using  the  model  developed  in  this  research.  It  simulates  a  10  s  irradiation  period 
and  displays  a  generally  linear  response  with  a  sawtooth-like  pattern  corresponding  to 
irradiation  and  annealing  during  a  pulse  cycle. 


When  considering  the  annealing  behavior,  the  annealing  factor  is  an  additional 
measure  of  the  transient  annealing  response  in  devices.  The  annealing  factor  is  de- 


hned  [33]  as 


AF  it)  = 


1/13  (t)  -  l/l3o 
1/I3f  -  \//3o 


The  initial  gain  in  the  device  is  /do,  which  acts  as  a  normalizing  factor  that  takes 


variability  in  device  parameters  into  account,  /3f  is  the  hnal  gain  in  the  device  at  long 
time,  and  (3  (t)  is  the  gain  at  a  particular  time,  t. 

The  hrst  trend  to  notice  in  Figures  26  and  27  is  that  A  (l//d)  approaches  its  hnal 
value,  i.e.  that  value  where  l/ft  — )■  1//3/  and  annealing  is  complete,  in  under  1  s. 
This  is  not  supported  by  experiments  or  literature.  However,  the  model  is  stable  and 


results  are  repeatable.  With  these  initial  results,  experimental  results  were  needed. 


76 


Annealing  Factor 


0.1245 

0.1244 

0.1244 

0.1243 

0.1243 

0.1242 

0.1242 

0.1241 


Time  After  Irradiation  [s] 


1.004 

1.0035 

1.003 

1 .0025 

1.002 

1.0015 

1.001 

1.0005 

1 


10'^  lO"'  lO" 
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Figure  27.  (Top)  The  simulated  change  in  gain  during  annealing.  (Bottom)  The  sim¬ 
ulated  annealing  factor  after  irradiation.  Annealing  in  the  model  occurs  quite  quickly 
and  approaches  the  final  value  within  less  than  1  s. 
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IV.  Experimental  Verification 


In  the  first  step  of  this  work,  a  time-dependent  system  of  coupled,  non-linear  ODEs 
was  solved  using  reaction  kinetics  to  determine  defect  species’  concentrations  as  they 
formed,  interacted,  and  were  annihilated.  The  solution  of  this  system  of  equations 
provides  insight  into  the  concentration  and  type  of  defects  that  affect  device  gain 
through  their  affect  on  the  base  dopant  and  minority  carrier  parameters. 

While  the  development  of  this  model  is  an  important  step  in  the  process  of  inves¬ 
tigating  temporal  and  spectral  dependence  in  neutron  displacement  damage  studies, 
the  next  step  is  to  use  experimental  measurements  to  validate  the  findings.  This 
data  will  help  to  predict  the  device  gain,  using  a  controlled  laboratory  environment. 
Additionally,  the  question  arises  as  to  whether  the  flux  temporal  and  spectral  prohles 
have  an  effect  on  the  device  performance.  This  can  be  investigated  by  comparing 
results  of  measurements  from  two  different  facilities:  the  accelerator-driven  neutron 
source  at  lU  NREF  and  a  steady-state  neutron  source  at  the  Ohio  State  University 
Research  Reactor  (OSURR).  The  facility  at  lU  NREF  has  a  much  harder  spectrum 
approximately  8  neutrons  for  every  gamma  ray,  while  the  thermal  spectrum  at  the 
OSURR  has  approximately  10  gamma  photons  emitter  per  neutron.  These  two  facil¬ 
ities  will  provide  a  good  contrast  in  both  time  and  neutron  energy  spectrum  for  this 
comparison.  It  is  expected  that  the  neutron  source  temporal  prohle  will  affect  the 
rate  at  which  damage  is  occurring,  and  that  the  neutron  energy  spectrum  shape  will 
affect  the  defect  species  being  created,  and  thus  the  damage  manifested  in  the  device. 

In  the  following  sections,  the  experimental  work  that  was  completed  during  this 
research  will  be  presented  in  detail,  to  include  the  procedures,  the  equipment  used, 
and  the  methodology  employed.  The  results,  which  show  that  device  performance 
degradation  during  irradiation  is  dehnitely  different  for  a  varying  flux  prohle,  will  be 
presented  and  analyzed. 
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4.1  lU  NREF  Experiments 


In  this  set  of  experiments,  the  lU  NREF  neutron  source,  which  is  a  continuously- 
pulsed  neutron  source,  was  used  to  damage  silicon  BJTs  and  the  device  gain  was 
recorded  to  measure  the  amount  of  displacement  damage  that  occurred  in  the  device. 
Three  types  of  experiments  were  accomplished:  an  interval  test,  a  multiple-fluence 
test,  and  an  ASTM  standard-based  test.  Gain  was  measured  in  two  ways,  as  described 
in  Tables.  First,  the  collector  and  base  were  sourced  with  a  constant  voltage,  the 
emitter  was  sourced  with  a  constant  current,  and  the  current  of  the  collector  and 
base  was  measured.  Second,  the  base  voltage  and  emitter  current  were  sourced,  but 
the  voltage  on  the  collector  was  swept  from  -1  to  -|-10  V,  while  the  base  and  collector 
currents  were  again  measured. 

Table  5.  BJT  operating  parameters  for  two  measurement  methods. 


Leg 

#1 

#2 

Collector  [V] 

+10 

Swept  from  -1  to  +10 

Base  [V] 

0 

0 

Emitter  [pA] 

-220 

-220 

The  results  of  these  experiments  show  that  the  device  gain  degradation  follows  a 
Messenger-Spratt  linear  relationship.  However  the  damage  constant  is  greater  than 
the  value  specihed  in  the  ASTM  E1855-10,  the  Standard  Test  Method  for  Use  of 
2N2222A  Silicon  Bipolar  Transistors  as  Neutron  Spectrum  Sensors  and  Displacement 
Damage  Monitors  [1] .  In  the  following  subsections,  the  facilities  at  lU  NREF  will  be 
explained  in  greater  detail,  the  equipment  used  for  each  of  the  experiment  sets  will 
be  described,  and  the  results  and  analysis  of  each  experiment  set  will  be  presented. 

Facility  Description. 

At  the  lU  Center  for  Exploration  of  Energy  and  Matter  (CEEM),  there  is  a  pulsed 
neutron  source  that  produces  neutrons  via  the  inelastic  scattering  of  protons  on  beryl- 
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Hum.  Protons  are  created  in  an  ion  sonrce,  then  are  accelerated  and  directed  onto 
a  target  assembly  inside  a  moderator/reflector  stack  where  nentrons  are  created.  A 
more  detailed  description  of  the  facility  can  be  fonnd  in  other  sonrces  [2-5],  and  a 
schematic  of  the  facility  is  displayed  in  Fignre28. 


Beam  line 
Selector  Magnet 


NREF 

Target 


Proton 

Source 


LIN  AC 


Figure  28.  The  lU  LINAC  Facility  schematic  shows  the  proton  source  and  LINAC  at 
bottom  left;  the  beam  is  split  to  targets  in  lower  right  and  upper  middle  [6]. 


On  a  basic  level,  the  lU  NREF  facility  is  made  np  of  a  linear  accelerator  that 
bombards  a  beryllinm  target  with  protons,  prodncing  nentrons  primarily  via  the 
(p,n)  inelastic  reaction.  The  linear  accelerator  is  capable  of  a  maximnm  energy  of  13 
MeV  nsing  one  radio-freqnency  qnadrnpole  and  two  drift  tnbe  Linear  Accelerators 
(LINAC).  The  LINAC  is  powered  by  klystrons  and  mnst  be  pnlsed  in  a  repetitions 
manner,  leading  to  a  distinct  dnty  factor.  This  dnty  factor  is  variable,  bnt  the  NREF 
is  typically  operated  with  a  dnty  factor  of  aronnd  0.8%.  The  freqnency  of  operation  is 
20  Hz  and  the  pnlse  width  of  the  proton  packets  is  between  100-1000  /rs.  An  example 
of  a  proton  pnlse  flnx  prohle  is  displayed  in  Fignre29. 

After  acceleration  throngh  the  RFQ  and  LINACS,  the  protons  can  be  directed 
down  one  of  two  beam  lines  nsing  a  selector  magnet.  The  existence  of  two  beam 
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Figure  29.  The  lU  NREF  beam  is  pulsed,  and  the  duty  cycle  is  typically  around  0.8% 
with  the  current  configuration.  Ideally,  the  beam  pulse  is  a  square  wave;  however,  each 
pulse  has  a  rise  and  fall  transient. 

lines  gives  the  facility  increased  flexibility  in  experiment  set-up  and  operation.  The 
hrst  beam  line  is  called  the  Low-Energy  Neutron  Source  (LENS)  beam  line,  while  the 
second  is  called  the  Neutron  Radiation  Effects  Facility  (NREF)  beam  line.  Neutron 
production  at  the  LENS  target  has  been  analyzed  extensively  [7],  however  the  mod¬ 
erator  conhguration  makes  the  LENS  target  impractical  for  use  in  radiation  effects 
studies  of  this  type. 

In  previous  work  by  this  author  and  others  [2-4,7-10],  the  neutron  energy  spectra 
for  the  LIN  AC’s  two  beamlines  (the  LENS  and  the  NREF)  has  been  characterized. 
The  spectrum  of  the  LENS  is  signihcantly  moderated,  i.e.  slowed  down,  so  that  the 
neutrons  may  interact  with  and  probe  the  structure  of  materials.  Since  one  of  our 
main  goals  in  doing  neutron  displacement  damage  studies  is  to  simulate  (as  closely 
as  possible)  the  environment  in  which  the  devices  may  have  to  operate,  this  facility 
is  not  of  great  interest.  Additionally,  the  extremely  low  energy  of  the  neutrons  of  the 
LENS  will  be  so  small  that  very  little  damage  will  occur.  It  is  for  these  reasons  that 
the  NREF  is  the  beamline  of  greater  interest  in  this  research. 


84 


Figure  30.  Using  Monte  Carlo  particle  simulation  software  and  an  input  deck  provided 
by  lU  NREF  personnel  [11],  the  neutron  spectrum  was  calculated  on  a  per  proton  basis 
and  at  different  distances  from  the  beryllium  target;  this  can  then  be  scaled  based  on 
the  proton  source  specifications  for  a  particular  run. 

The  NREF  has  been  modified  recently  to  accommodate  the  radiation  effects  on 
electronics  community.  The  spectrum  was  hardened,  meaning  that  the  average  neu¬ 
tron  energy  is  skewed  to  higher  energy;  and  the  maximum  neutron  energy  is  approxi¬ 
mately  11  MeV.  In  Figure  30,  the  neutron  energy  spectrum  is  plotted  as  a  function  of 
energy  and  per  a  single,  incident  proton.  In  this  way,  the  spectrum  can  be  scaled  for 
any  proton  beam  current,  which  is  dependent  on  the  beam’s  operating  parameters. 
Usually,  the  beam  is  operated  at  a  current  of  20  mA,  and  this  scaled  spectrum  is 
plotted  in  Figure  31. 

NREF  Experiment  1:  Exploration  Experiments. 

The  first  set  of  experiments  at  the  lU  NREF  was  an  investigation  of  the  feasibility 
of  the  experimental  method  devised  to  collect  BJT  gain  degradation  data.  To  this 
end,  a  number  of  procedures  were  tested  for  use  in  the  second  set  of  experiments. 
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Figure  31.  When  operating  the  proton  source  at  a  current  of  20  mA,  the  resultant 
neutron  spectrum  is  scaled  to  a  1-MeV  equivalent  (Si)  neutron  flux  of  >  6x10^°  n/cm^/s. 

Test  Equipment. 

The  equipment  used  for  this  experiment  included  three  Keithley  237  Source  Mea¬ 
surement  Units  (SMU),  a  switchbox,  and  a  test  fixture.  The  latter  two  pieces  of 
equipment  were  built  specihcally  for  this  experiment.  Figure  32  displays  a  picture  of 
the  equipment  set-up  for  the  irradiation  portion  of  the  experiment. 


Figure  32.  A  close-up  of  the  DuT  fixture  (left),  switchbox  (middle),  and  SMUs  (right). 
In  the  text  fixture  and  switchbox,  the  aluminum  box  acted  as  an  EMI  shield. 


Using  the  standard  DuT  mount,  the  1-MeV  neutron  flux  at  the  lU  NREF  is 
2  X  10^°  n/cm^/s.  A  modification  was  made  to  increase  the  flux  by  using  an  acrylic 
jack  to  reposition  the  DuT  approximately  4-6  cm  from  the  beryllium  target.  This 
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modifcation  increased  the  flux  by  as  much  as  a  factor  of  2.5.  The  jack  was  supported 
by  a  6  inch  tall  block  of  acrylic.  This  positioned  the  BJTs  along  the  neutron  source 
center  line,  maximizing  the  flux  incident  on  the  device  by  minimizing  streaming  loss. 
A  picture  of  this  set-up  is  displayed  in  Figure  33. 


Figure  33.  The  test  set-up  inside  the  TMR  included  a  block  of  acrylic  to  position  the 
DuT  fixture  at  the  centerline  of  the  beam.  The  acrylic  jack  sat  atop  this  block.  The 
aluminum  test  rig  is  positioned  to  the  left  of  the  jack  and  is  facility  the  beryllium  target 
(which  is  out  of  sight  below  and  to  the  left  of  the  set-up). 


Test  Procedures. 

Two  types  of  experiments  were  accomplished:  the  Multiple  Sequential  Interval 
(MSI)  Test  and  the  Single  Long  Interval  (SLI)  Test.  These  two  experiments  were 
conducted  on  the  same  day  with  the  same  beam  parameters  in  order  to  minimize 
variations  in  the  neutron  flux  spectrum. 

The  gamma-induced  photocurrent  was  mitigated  by  using  a  second  shunting  BJT 
[12].  Figure  34  displays  a  schematic  of  the  experimental  circuit  in  the  common  base 
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configuration.  The  shunt  BJT  was  selected  to  most  closely  match  the  DuT  gain.  The 
collector  of  the  shunting  BJT  was  connected  to  the  base  in  the  test  BJT  to  ensure  that 
the  base-collector  junction  remained  reverse-biased  during  the  irradiation.  Because 
the  gamma  flux  from  the  NREF  is  relatively  low,  especially  when  compared  with  a 
nuclear  reactor,  its  effect  was  minimal. 


Test 

BJT 


Figure  34.  A  schematic  of  the  shunt  and  test  BJTs.  The  BJT  being  analyzed  is  in 
Common  Ground  configuration. 


In  the  MSI  experiment,  the  beam  was  set  to  the  parameters  found  in  Table  6.  For 
the  initial  interval,  the  device  was  irradiated  until  the  gain  was  reduced  to  75%  of 
the  range  between  its  starting  value  and  30.  Then,  the  beam  was  shut  off  for  30  s  of 
annealing.  This  was  repeated  three  more  times,  and  the  gain  was  reduced  by  25%  of 
the  range  each  time  until  the  gain  was  30  after  the  fourth  irradiation  interval.  Then, 
the  device  was  allowed  to  anneal.  For  example,  the  BJT  pre-irradiation  gain  was 
147.4.  The  beam  was  shut  down  when  the  device  gain  was  118,  88.5,  59,  and  stopped 
when  it  reached  30. 

Table  6.  Beam  parameters  on  the  day  of  experiment. 


Beam  Parameter 

Value 

Current  [mA] 

20 

Max.  Proton  Energy  [MeV] 

13 

Frequency  [Hz] 

20 

Pulse  Width  [ps] 

400 

The  SLI  experiment  followed  the  MSI  experiment.  The  beam  was  turned  on  for 
a  specihed  interval  of  time  at  the  same  flux,  thereby  corresponding  to  a  desired  total 


fluence.  The  BJT  being  irradiated  was  constantly  biased  with  +10  V  on  the  collector, 
0  V  on  the  base,  and  -220  /iA  on  the  emitter.  After  the  desired  time  was  reached,  the 
beam  was  turned  off  while  measurements  of  the  powered  devices  continued.  After  an 
interval  of  measuring  the  device  annealing,  power  from  the  SMU  was  turned  off  and 
the  devices  were  removed  from  the  TMR  vault. 

Neutron  Dosimetry. 

The  neutron  source  energy  spectrum  has  been  quantihed  via  Monte  Carlo  sim¬ 
ulations  using  MCNP-X  Version  2.7.0  [13]  along  with  verihcation  using  sulfur  foil 
activation  analysis. 

Table  7.  Activation  analysis  results  from  sulfur  foils. 


Irr  Time  (min.) 

4^7neasured 
\  "  1 

4>mcnp 

\  "  1 
Lcm^-sJ 

Ratio 

Avg. 

30 

3.29  X  10^^ 

3.78  X  10^° 

3.55  X  10^^^ 

3.55  X  10^° 

0.927 

1.065 

0.996 

60 

3.26  X  10^° 
3.47  X  10^° 

3.55  X  10^° 
3.55  X  10^° 

0.918 

0.977 

0.948 

90 

3.59  X  lO^^^ 
3.27  X  10^° 

3.55  X  10^^^ 

3.55  X  10^° 

1.011 

0.921 

0.966 

Table  7  displays  the  results  of  neutron  activation  analysis  completed  for  the  30, 
60,  and  90  min  irradiations.  Sulfur  provided  good  agreement  with  simulation  (within 
5%,  which  is  a  common  factor  of  systemic  uncertainty  for  neutron  activation  analysis) 
so  that  the  simulated  neutron  energy  spectrum  can  be  accepted.  Thus,  the  resulting 
1-MeV  equivalent  (Si)  neutron  flux  for  the  new  position  dictated  by  the  acryllic  jack 
was  found  to  be  4.42(±0.22)  x  10^°  n/cm^/s,  or  about  a  factor  of  2  larger  than  the 
reported  value  of  flux  for  the  standard  position  without  the  jack. 


Results  and  Analysis. 


The  MSI  experiment  raw  results  are  displayed  in  Figure  35.  The  SLI  experiment 
raw  results  are  presented  in  Figure  36. 


Figure  35.  BJT  gain  as  a  function  of  time  for  the  MSI  test.  The  intervals  of  irradia¬ 
tion  and  intervals  of  annealing  were  fit  to  the  Messenger-Spratt  equation  and  a  power 
function,  respectively. 


The  gain  for  each  interval  of  irradiation  was  £t  to  the  Messenger-Spratt  equation 
using  the  equation, 


Pit) 


(47) 


where  D  is  the  fitting  coefficient,  which  encompasses  both  the  damage  constant  K  [E) 
and  the  flux  0  [E).  The  first  three  intervals  of  annealing  were  30  s  long.  Not  seen  in 
Figure  35  is  the  final  interval  of  annealing,  which  was  also  30  s.  Annealing  was  found 
to  follow  a  power  law  function  as  in  Equation  48, 


/3  (t)  =  At^  +  B 


(48) 


90 


(a) 


(b) 


(c) 


(d) 


Figure  36.  Gain  as  a  function  of  time  for  continuous  irradiations  of  (a)  20  min,  (b) 
30  min,  (c)  60  min,  and  (d)  90  min.  Note  that  in  (a)  there  is  a  significant  annealing 
time  (approximately  60  min),  while  in  the  other  three  cases  (b-d)  there  is  only  a  5  min 
annealing  measured  following  irradiation. 
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where  A,  B,  and  C  are  fitting  coefficients.  By  inspection,  B  is  the  device  gain  at  the 
start  of  the  annealing  interval. 

After  nonlinear  regression  was  accomplished  and  the  coefficients  were  determined, 
D  provided  an  estimate  of  the  damage  constant,  K,  and  how  it  compared  with  the 
accepted  value  in  literature,  Ki^ueV  ~  1-5  x  10“^^  cm^/n.  This  information  is 
tabulated  in  Table  8. 

Table  8.  Fitting  coefficients  from  gain  degradation  intervals  with  extracted  damage 
constant,  K. 


7.55x10-^ 

19.6% 

1.72x10-15 

5.84x10-® 

6.95x10-^ 

18.3% 

1.58x10-15 

1.09x10-5 

7.09x10-5 

6.6% 

1.61x10-15 

4.57x10-® 

6.70x10-5 

2.4% 

1.52x10-15 

4.65x10-® 

7.31x10-5 

0.04% 

1.65x10-15 

1.65x10-^ 

8.30x10-5 

0.05% 

1.88x10-15 

5.09x10-^ 

7.87x10-5 

0.02% 

1.78x10-15 

1.56x10-^ 

6.81x10-5 

0.02% 

1.54x10-15 

1.60x10-^ 

7.32x10-5 

47.03% 

1.66x10-15 

1.27x10-^ 

The  same  nonlinear  regression  fitting  was  accomplished  for  annealing,  however 
it  was  fit  to  a  power  function  as  presented  in  Equation  48.  The  results  of  this  are 
displayed  in  Table  9. 


Table  9.  Fitting  coefficients  from  annealing  intervals. 


Interval 

A 

Err. 

B 

Err. 

C 

Err. 

1-R2 

First 

2.9  X  10-1 

6.7% 

116.2 

0.022% 

0.45 

3.3% 

4.0  X  10-® 

Second 

1.9  X  10-2 

5.8% 

88.21 

0.018% 

0.50 

2.7% 

3.6  X  10-5 

Third 

6.2  X  10-2 

4.1% 

60.07 

0.008% 

0.69 

1.5% 

1.6  X  10-^ 

Final 

3.0  X  10-2 

4.7% 

30.22 

0.008% 

0.62 

1.9% 

8.5  X  10-9 

20  min 

2.2  X  10-5 

5.5% 

10.33 

0.003% 

0.79 

1.8% 

3.4  X  10-9 

60  Min 

2.6  X  10-^ 

12.0% 

3.405 

0.002% 

0.84 

3.8% 

2.3  X  10-9 

90  Min 

2.8  X  10-^ 

10.3% 

2.652 

0.002% 

0.71 

3.8% 

1.2  X  10-9 

The  fitting  coefficient  B  is  observed  to  be  the  B  JT  gain  at  the  start  of  the  annealing 
interval.  The  A  and  C  coefficients  follow  a  general  downward  trend  as  irradiation  time 
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Figure  37.  Examples  of  the  curve  fitting  for  both  degradation  (top)  and  annealing 
(bottom).  Note  that  in  the  top  chart,  gain  was  normalized  to  the  precharacterization 
gain,  while  on  the  bottom  it  has  not  been  normalized. 
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Time  [s]  Time  [s] 

Figure  38.  The  A  and  C  coefficients  are  plotted  as  a  function  of  time.  The  connecting 
dotted  lines  in  the  figures  are  just  for  presentation  and  analysis. 

increases.  If  it  is  possible  to  determine  the  functionality  of  A  and  C  with  fluence  and 
annealing,  it  may  lead  to  a  better  understanding  of  how  to  incorporate  relaxation  into 
the  model.  However,  this  still  has  not  been  resolved  and  requires  further  investigation. 
The  coefficients  are  plotted  as  a  function  of  time  in  Figure  38. 

In  Figure  39,  a  comparison  of  the  change  in  gain  relative  to  the  beginning  of  the 
annealing  interval  is  presented  as  a  function  of  time  from  the  start  of  the  annealing 
interval  for  the  MSI  experiment.  In  this  format,  one  can  see  that  annealing  in  the 
BJT  devices  occurs  faster  when  the  device  has  been  damaged  less.  This  is  verihcation 
of  a  prediction  by  the  mathematical  model  developed  in  Chapter  III. 

After  removing  the  devices  from  the  lU  NREF  vault,  three  of  the  irradiated  BJTs 
were  placed  in  a  laboratory  oven  and  baked  from  at  150°C.  Afterward,  they  were 
removed  from  the  oven  and  allowed  to  cool  to  room  temperature.  During  this  entire 
process,  the  gain  was  measured  in-situ  using  an  HP4145A  Semiconductor  Parameter 
Analyzer  [14].  The  results,  after  being  normalized,  are  plotted  in  Figure 40. 

Comparison  with  traditional  single-pulse  system. 

One  of  the  main  interests  in  this  research  was  to  compare  the  NREF  neutron 
source  with  other  neutron  sources.  One  of  the  best  ways  to  do  this  is  to  obtain 
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Figure  39.  Comparison  of  the  annealing  of  BJT  gain  for  multiple  sequential  interval 
(top)  and  single  long  interval  (bottom)  tests.  The  y-£Lxis  is  the  change  in  gain  relative 
to  the  start  of  annealing  interval. 
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Figure  40.  (Top)  Gain  as  a  function  of  temperature  during  a  post-irradiation  thermal 
anneal.  Gain  was  measured  in-situ  while  being  baked  and  normalized  to  the  gain  at 
150°C.  (Bottom)  Gain  as  a  function  of  temperature  as  the  devices  were  allowed  to 
come  back  to  room  temperature  after  the  thermal  anneal.  Gain  in  both  plots  has  been 
normalized  to  the  maximum  value  (at  150°C). 


gain  degradation  data  at  multiple  1-MeV  equivalent  (Si)  neutron  fluences  and  make 
a  direct  comparison  of  the  change  in  gain,  A  (1//3). 

Figure  41  displays  a  plot  that  compares  data  as  collected  from  four  widely-used 
and  well  characterized  single  pulse  neutron  sources  [12].  These  four  systems  are  the 
Sandia  Pulsed  Reactor  (SPR-III),  White  Sands  Missile  Ranges  Fast  Burst  Reactor 
(WSMR  FBR),  the  Annular  Core  Research  Reactor  (ACRR),  and  the  Los  Alamos 
Neutron  Science  Center  (LANSCe).  These  four  systems  are  traditional,  single  pulse, 
high  intensity  neutron  sources  that  provide  neutrons  with  a  small  pulse  width.  Data 
from  these  systems  match  well  with  the  typically  accepted  value  for  the  damage 
constant  for  2N2222A  BJTs  in  the  Messenger-Spratt  Equation,  from  1  to  2x10“^^ 
cm^/n. 


_2 

Fluence  [n  cm’  ] 


X  10 


Figure  41.  A  comparison  of  gain  degradation  as  a  function  of  fluence  for  multiple  sys¬ 
tems.  ’’Legacy  data”  includes  data  from  SPR-III,  WSMR  FBR,  ACRR,  and  LANSCe. 
The  dotted  line  represents  the  A  (1/d)  using  the  accepted  value  of  the  damage  constant, 

K. 


However,  when  the  lU  NREF  data  are  plotted  along  with  the  legacy  data,  the 
device  degradation  at  higher  fluences  is  different  due  to  annealing.  In  all  cases,  the 
NREE  data  fall  below  the  expected  values;  thus,  the  device  final  gain  is  higher  than 
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would  be  expected  to  correspond  to  the  neutron  fluence  levels.  Unfortunately,  each  of 
these  lU  NREF  data  points  in  Figure  41  is  a  single  measurement,  and  it  could  just  be 
experimental  or  systemic  error.  Therefore,  a  second  set  of  experiments  was  planned  at 
lU  NRFF  with  more  rigor  applied  and  more  devices  irradiated  to  get  better  statistics. 

NREF  Experiment  2:  Expanded  Dataset. 

The  purposes  of  the  second  set  of  experiments  at  the  lU  NREF  beamline  were 
twofold.  First,  the  second  set  of  experiments  employed  more  rigorous  dosimetry,  to 
include  additional  neutron  activation  foils  as  well  as  thermoluminescent  dosimeters 
(TLD)  as  prescribed  by  ASTM  standard  methods.  Second,  better  statistics  were 
obtained  by  irradiating  more  devices,  thereby  expanding  the  dataset  for  analyiss. 

As  was  done  during  the  first  experimental  series  at  lU  NREF,  BJTs  were  irradiated 
by  neutrons  and  the  performance  degradation  was  measured  in-situ.  In  the  following 
sections,  the  equipment  used,  procedures  followed,  and  results  obtained  are  presented. 

Test  Equipment. 

During  irradiation,  the  BJTs  being  tested  were  under  constant  electrical  bias.  To 
accomplish  this,  three  Keithley  237  SMUs  were  set  up  inside  the  vault  but  outside 
the  TMR.  Each  BJT  leg  was  connected  to  an  SMU  via  approximately  10  to  15  ft  of 
BNC  cabling.  The  SMU  output  was  fed  to  a  laptop  that  was  running  a  custom-built 
LabView  virtual  instrument,  via  General  Purpose  Interface  Bus  (GPIB)  cables  and  a 
National  Instruments  Ethernet-to-GPIB  converter. 

As  with  the  first  experiment,  a  custom-built  acrylic  jack  was  placed  on  the  TMR 
positioning  platform  inside  the  open-air  cavity  to  properly  position  the  DuT  for  max¬ 
imum  neutron  flux.  Additionally,  the  BJTs  being  irradiated  were  contained  within 
a  custom-built  aluminum  enclosure  for  EMI  protection.  This  enclosure  used  SMA 


through-bulkhead  connectors  to  connect  to  the  measurement  equipment  positioned 
outside  the  TMR,  but  inside  the  vault.  Also  similar  to  the  first  experiment,  a  sec¬ 
ond  BJT  was  connected  to  the  base  leg  of  the  original  device  in  order  to  shunt  the 
gamma-induced  photocurrent. 

Test  Procedures. 

For  all  irradiations,  the  BJT  for  which  gain  was  being  measured  was  under  con¬ 
stant  bias  with  J-IOV  on  the  collector,  -220  A  on  the  emitter,  and  the  base  was  set  to 
SMU  ground.  The  collector  and  base  SMUs  were  set  to  source  voltage  and  measure 
current,  while  the  emitter  SMU  was  set  to  source  current  and  measure  voltage.  This 
conhguration  ensured  that  the  collector-base  junction  was  reverse-biased,  while  the 
emitter-base  junction  was  forward-biased,  causing  the  BJT  to  be  in  forward  active 
mode  at  all  times.  It  is  in  this  mode  that  BJTs  are  designed  to  provide  the  greatest 
current  gain;  this  also  means  that  gain  sensitivity  to  changes  in  the  defect  species 
populations  is  maximized. 

For  four  experiment  runs,  the  lU  NREF  beamline  was  switched  on  with  a  nominal 
beam  current  of  20  mA,  a  pulse  width  of  400  s,  and  a  beam  repetition  rate  of  20  Hz. 
The  SMU  output  from  all  three  legs  was  collected  using  a  custom-built  LabView 
virtual  instrument.  After  20  min,  the  beamline  was  switched  off  but  the  SMU  output 
continued  recording  for  5-10  min  afterward.  The  purpose  of  this  experiment  set  was  to 
ensure  a  consistent,  repeatable  experiment  that  complies  with  community  standards, 
specihcally  the  American  Society  for  Testing  and  Materials  (ASTM)  standards  such  as 
E265-07(2013)  Standard  Test  Method  for  Measuring  Reaetion  Rates  and  Fast-Neutron 
Fluenees  by  Radioaetivation  of  Sulfur-32  [15]  and  E1855-10  Standard  Test  Method 
for  Use  of  2N2222A  Silieon  Bipolar  Transistors  as  Neutron  Speetrum  Sensors  and 
Displaeement  Damage  Monitors  [1]. 
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In  a  second  set  of  irradiations,  devices  were  exposed  to  diffferent  total  neutron 
fluences,  starting  with  40  and  increasing  to  60,  90,  and  120  min  total  irradiation 
periods.  BJTs  were  irradiated  and  the  gain  degradation  measured  in-situ  before, 
during,  and  after  irradiation.  Because  the  beam  capacitor  faulted  down  during  the 
hrst  40  min  irradiation,  an  additional  irradiation  period  of  4  min  was  added  to  the 
experiment.  The  neutron  flux  from  the  hrst  set  of  experiments  at  the  lU  NREF 
was  4.4  X  10^°  n/cm^/s  when  the  DuT  was  approximately  4-6  cm  from  the  beryllium 
target,  which  corresponds  to  huences  of  0.96,  1.44,  2.16,  and  2.88  x  10^^  n/cm^/s. 

The  lU  NREF  beamline  operational  parameters  for  this  experiment  was  a  nominal 
beam  current  of  20  mA,  a  pulse  width  of  400  s,  and  a  beam  repetition  rate  of  20  Hz. 

The  SMU  output  from  all  three  legs  was  collected  using  the  same  custom-built 
LabView  program.  After  the  specihed  time,  the  beamline  was  switched  off  but  the 
SMU  output  continued  recording  for  5-10  min  afterward.  The  purpose  of  this  ex¬ 
periment  set  was  to  verify  the  data  points  on  Figure  41  and  hll  in  additional  huence 
values  to  better  understand  the  lU  NREF  damage  coefficient. 

Neutron  Dosimetry. 

The  hrst  set  of  4  BJTs  was  irradiated  for  20  min  each.  One  S  pellet,  one  Fe  foil, 
and  one  TLD  were  positioned  within  6  cm  of  the  BJT  being  irradiated  to  determine 
the  neutron  and  gamma  doses  incident  on  the  BJTs.  Before  irradiation,  the  mass  of 
the  S  pellets  and  Fe  foils  used  in  neutron  activation  analysis  were  measured  using  a 
standard  digital  laboratory  scale. 

This  information  is  tabulated  in  column  3  of  Table  10  and  Table  11.  The  average 
S  mass  was  0.590  ±  0.01  g  while  the  average  Fe  mass  was  0.303  ±  0.003  g.  This 
information  was  then  used  to  calculate  the  atomic  number  density  of  the  relevant 
isotopes  for  the  sulfur  pellets  and  ®^Fe  and  ^®Fe  for  the  iron  foils).  Number 
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density,  iV^,,  is  where  is  the  density  of  the  isotope,  Nav  is  Avagadros  number 
at  6.022  X  10^^  at/mol,  and  is  the  atomic  mass  of  the  isotope.  The  atomic 
mass  of  S  and  Fe  are  assumed  to  be  32.065  ±  0.005  g/mol  and  55.845  ±  0.002  g/mol, 
respectively. 


Table  10.  Iron  foil  mass  measurement. 


# 

Mass  (g) 

+/- 

N_Fe-54 

+/- 

Npe-se 

+/- 

1 

0.2998 

1.89x1020 

8.60x10^0 

2.97x1021 

1.26x101® 

2 

0.3062 

0.0001 

1.93x1020 

8.65x10^0 

3.03x1021 

1.26x101® 

3 

0.3036 

1.91x1020 

8.63x10^0 

3.00x1021 

1.26x101® 

4 

0.3002 

1.89x1020 

8.60x10^0 

2.97x1021 

1.26x101® 

Table  11.  Sulfur  pellet  mass  measurement. 


# 

Mass  (g) 

+/- 

1 

CO 

to 

+/- 

1 

0.5913 

1.05x1022 

3.98x101® 

2 

0.599 

0.0001 

1.07x102^ 

4.01x101® 

3 

0.5972 

1.07x102^ 

4.01x101® 

4 

0.574 

1.02x102^ 

3.92x101® 

The  Fe  foil  is  used  to  best  measure  the  epithermal  region  of  a  neutron  spectrum, 
that  is  the  region  between  thermal  (<0.25  eV)  and  fast  (>1  MeV)  neutron  energies. 
There  are  two  main  reactions  of  interest  in  Fe  foils  in  the  spectrum  at  the  lU  NREF, 
which  are  described  in  the  following  list  along  with  some  of  the  pertinant  information 
about  these  reactions: 

1.  ^'^Fe{n,pY^Mn  -  There  is  a  low  percentage  of  atoms  in  each  iron  sample, 
and  the  834.8  keV  gamma  ray  can  be  difficult  to  differentiate  from  gamma 
emitted  by  ®®Mn. 

2.  ^^Fe{n,p)^^Mn  -  This  isotope  of  manganese  has  two  dominant  gamma  rays  at 
846.8  keV  and  1810.7  keV.  The  second  (photopeak)  is  easily  identihed  and  used 
in  calculation  of  the  radioactivity. 
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Table  12.  Gamma  peak  information  for  the  iron  foils. 


# 

Epea/c 

a  EpCak 

C  Frac. 

a  C 

Peak  Area 

PA  ± 

1 

834.838 

0.003 

0.9998 

0.001 

1.49x103 

38.71 

846.7638 

0.0019 

0.989 

0.3 

3.03x103 

54.32 

1811 

0.004 

0.287 

0.8 

3.57x10^ 

22.66 

2 

834.838 

0.003 

0.9998 

0.001 

1.41x103 

38.57 

846.7638 

0.0019 

0.989 

0.3 

3.75x103 

60.48 

1811 

0.004 

0.287 

0.8 

4.15x10^ 

24.08 

3 

834.838 

0.003 

0.9998 

0.001 

1.37x103 

37.78 

846.7638 

0.0019 

0.989 

0.3 

3.65x103 

59.8 

1811 

0.004 

0.287 

0.8 

4.18x103 

23.04 

4 

834.838 

0.003 

0.9998 

0.001 

9.85x10^ 

32.37 

846.7638 

0.0019 

0.989 

0.3 

2.82x103 

52.68 

1811 

0.004 

0.287 

0.8 

3.28x10^ 

21.79 

The  use  of  Fe  foils  as  a  neutron  activation  analysis  dosimeter  is  delineated  by  the 
ASTM  Standard  E263-13,  Standard  Test  Method  for  Measuring  Fast-Neutron  Reac¬ 
tion  Rates  by  Radioactivation  of  Iron  [16].  To  calculate  and  correct  for  the  activity 
at  inhnite  irradiation  time,  the  total  counts  under  the  peak  is  used  in  Equation  49: 


A(C'  —  B)  exp 

e  [1  -  exp  {-AtJ]  [1  -  exp  ’ 


(49) 


where  7  is  the  isotopic  decay  constant;  C  and  B  are  the  measured  counts  under 
the  peak  and  background  counts  in  the  detector,  respectively;  ti,  tw,  and  tc  are  the 
irradiation,  waiting,  and  counting  times,  respectively;  and  e  is  the  detector  efficiency 
factor.  Figure  42  displays  a  visual  description  of  the  respective  times  for  Equation  49. 
The  tabulated  results  of  this  analysis  are  displayed  in  Table  13. 


t  t  t  t 

Irr.  Irr.  Count  Count 

Start  Stop  Start  Stop 


Figure  42.  Timeline  representing  irradiation,  waiting,  and  counting  times  for  neutron 
activation  analysis. 
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Table  13.  Calculation  of  measured  1-MeV  Equivalent  (Si)  neutron  flux  using  iron  foils. 


# 

Epea/c 

[keV] 

C 

e 

Aoo 

[dps] 

^isotope 

[cm~^] 

Aoo 

^ isotope 

^eff 

[cm-2s-i] 

1 

835 

1.85x10^ 

1.67x10® 

1.89x102® 

8.85x10-1® 

6.00x101® 

847 

3.89x10^ 

1.01x10® 

2.97x1021 

3.40x10-1® 

2.31x10® 

1811 

4.03x10^ 

1.05x10® 

2.97x1021 

3.52x10-1® 

2.39x10® 

2 

835 

1.75x10^ 

1.51x10® 

1.93x102® 

7.81x10-1® 

5.29x101® 

847 

4.82x10^ 

9.26x10® 

3.03x1021 

3.06x10-1® 

2.07x10® 

1811 

4.69x10^ 

9.01x10® 

3.03x1021 

2.97x10-1® 

2.02x10® 

3 

835 

1.70x10^ 

1.54x10® 

1.91x102® 

8.03x10-1® 

5.45x101® 

847 

4.69x10^ 

9.84x10® 

3.00x1021 

3.28x10-1® 

2.22x10® 

1811 

4.72x10^ 

9.91x10® 

3.00x1021 

3.30x10-1® 

2.24x10® 

4 

835 

1.22x10^ 

1.10x10® 

1.89x102® 

5.82x10-1® 

3.95x101® 

847 

3.62x10^ 

6.60x10® 

2.97x1021 

2.22x10-1® 

1.51x10® 

1811 

3.70x10^ 

6.75x10® 

2.97x1021 

2.27x10-1® 

1.54x10® 

To  determine  the  effective  ffux,  ^e//,  the  activity  calculated  in  Equation 49  was 
used  via  Equation  50. 

^oo  =  ^isotope'^isotope^’ef  f  ■  (50) 

In  Equation  50,  Nisotope  is  the  number  density  of  the  parent  isotope  and  cr^^otope  i®  the 
spectrum-averaged  microscopic  cross  section  of  the  reaction  of  interest  in  the  parent 
isotope  that  creates  the  radioactive  daughter.  The  spectrum-averaged  cross  section 
is  defined  as 

/  u  {E')  0  {E')  dE' 

^isotope  j^[E')dE'  ’  ^  ’ 

E 

where  a  {E')  is  the  energy-dependent  cross-section  and  0  {E')  is  the  energy-dependent 
neutron  spectrum  determined  using  MCNP.  From  MCNP  calculations,  the  expected 
spectrum-averaged  cross  section  for  Fe  is  1.47  x  10  cm^  or  0.147  b. 

The  S  pellet  was  used  to  measure  the  fast  region  of  a  neutron  spectrum.  There  is 
one  reaction  of  interest  in  S  pellets,  which  is  the  (n,p)  inelastic  collision  resulting  in 
radioactive  ^^P.  This  phosphorus  isotope  decays  by  beta  emission  with  a  maximum 
energy  of  1709  keV.  At  lU,  the  induced  activity  in  the  S  pellets  after  irradiation  was 
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Figure  43.  Beta  organic  scintillation  counter  at  lU  (left  arrow)  for  measurement  of 
beta  decay  products  from  phosphorus  daughter;  S  pellet  was  placed  in  an  aluminum 
holding  fixture  (right  arrow)  14  cm  from  the  detector. 


measured  using  a  scintillation  detector  (displayed  in  the  left  portion  of  Figure  43). 
Before  measuring  the  irradiated  S,  a  series  of  four  5- min  background  measurements 
was  made  to  determine  an  average  background  count,  which  is  tabulated  in  Table  14. 
The  average  5-min  background  measurement  was  2478  (±50)  counts. 


Table  14.  Sulfur  pellet  background  measurement  for  5  minute  counting  intervals;  note 
that  Cps  =  counts  per  second. 


# 

Counts 

Avg  Counts 

Counts  ± 

Cps 

Avg  Cps 

Avg  Cps  ± 

b.l 

2716 

9.05 

b.2 

2381 

2477.5 

49.77 

7.94 

8.26 

0.55 

b.3 

2457 

8.19 

b.4 

2356 

7.85 

Additionally,  a  104  year  half-life,  standardized  ^^Si  source  was  used  to  calibrate 
detector  efficiency.  This  source  was  produced  by  Eckert  &  Ziegler  and  is  NIST- 
traceable.  The  standard  source  was  placed  in  the  apparatus  and  counted  for  5  min. 
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Total  beta  counting  efficiency  was  then  calculated  as, 


^beta  >S*>S*  ■ 

C std^src  B 


=  ss- 


i-A., 


(52) 


where  S' S'  is  the  self-shielding  factor;  Sdet  is  the  detector  efficiency;  Cstd_src  and  B  are 
the  standard  source  and  background  counts,  respectively;  is  a  geometric  correction 
factor  accounting  for  solid  angle  of  the  standard  source  versus  the  irradiated  S  pellet 
samples;  and  Agource  is  the  activity  of  the  standard  source  in  units  of  decays/min. 

The  self-shielding  factor  corrects  for  betas  emitted  from  within  the  S  pellet  that 
must  travel  through  the  S  before  being  emitted  from  the  pellet.  A  percentage  of 
the  betas  are  captured  by  the  pellet  as  they  travel,  reducing  the  measured  value  of 
radioactivity.  By  including  this  self-shielding  factor,  which  is  always  <  1,  we  can 
correct  for  this  reduction.  From  literature  and  by  using  MCNP,  the  best  value  of  SS 
for  sulfur  is  0.88,  which  was  subsequently  used  in  this  dosimetry  analysis. 

The  geometric  factor  ^  corrects  for  the  difference  between  the  calculation  of  effi¬ 
ciency  in  the  standard  source  and  the  calculation  of  efficiency  in  the  S  pellet  source.  In 
the  calculation  of  geometric  detector  efficiency,  the  correction  factor  is  ^det  = 

■^source 

if  the  standard  source  has  the  same  physical  dimensions  as  the  S  pellet.  This  is  equiv¬ 
alent  to  saying  that  .^  =  1.  In  the  case  of  this  experiment,  this  is  not  true,  so  the 
geometric  factor  ^  must  be  added  in  the  denominator.  To  calculate  solid  angle,  hi, 
for  a  detector  with  radius  a  at  a  distance  d  from  a  source  of  radius  s,  Equation  Set  53 
was  used. 
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Figure  44.  Schematic  showing  variables  involved  in  calculating  solid  angle  for  a  source 
of  beta  particles  with  radius,  s,  a  distance,  d,  away  from  a  detector  with  radius,  a.  The 
geometric  factor,  is  dependent  on  the  solid  angle. 


The  geometric  factor  is  thus  the  ratio  of  solid  angle  for  the  standard  source  ^std^src 
to  the  solid  angle  for  the  irradiated  S  pellet  sample  sample  or  ^  .  The 

detector  face  radius  is  a  =  3  cm,  the  distance  from  the  beta  source  is  d  =  14  cm, 
the  standard  source  radius  is  Sstd^src  =  1-02  cm,  and  the  S  pellet  sample  radius  is 
Ssam  =  0.64  cm.  Therefore,  the  geometric  factor  is  .^  =  0.998.  Therefore,  since  the 
average  background  counts  are  2478,  the  calculated  detector  efficiency  is  found  as 


. 2.204  X  10 

0.99773  ■  63180  J  63037 


(54) 
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This  leads  to  a  total  system  efficiency  of 


ettta  =  0.88  ■  2.204  x  10“^  =  1.940  x  10“^ 


(55) 


Each  irradiated  S  pellet  was  measured  for  four  5-min  counting  intervals.  The  ac¬ 
tivity  was  corrected  to  inhnite  irradiation  time  using  Equation  49.  The  results  are  tab¬ 
ulated  in  Table  15.  The  effective  flux  was  determined  as  (pavg  =  ^oo  [Nisotope<^ isotope]  ^ ■ 
The  spectrum-averaged  cross  section  for  S  is  1.23  x  10“^^  cm^  or  0.123  b. 


Table  15.  Calculation  of  measured  1-MeV  Equivalent  (Si)  neutron  flux  using  sulfur 
pellets. 


C-B 

Aoo 

4^00 

^eff 

ID# 

[dps] 

[cm~^] 

^isotope 

[cm-2s-i] 

1.1 

5.041x10^ 

2.625x10^ 

1.054x10^^ 

2.490x10-1® 

2.023x101^ 

1.2 

4.973  xlO^ 

2.590x10^ 

1.054x10^2 

2.456x10-1® 

1.996x10® 

1.3 

4.976  xlO^ 

2.593x10^ 

1.054x1022 

2.459x10-1® 

1.998x10® 

1.4 

4.889  xlO^ 

2.548  xlO^ 

1.054x1022 

2.416x10-1® 

1.963x10® 

2.1 

6.476  xlO^ 

3.201x10® 

1.068x1022 

2.998x10-1® 

2.436x10® 

2.2 

6.419x10^ 

3.174x10® 

1.068x1022 

2.972x10-1® 

2.415x10® 

2.3 

6.495x10^ 

3.213x10® 

1.068x1022 

3.008x10-1® 

2.444x10® 

2.4 

6.321x10^ 

3.127x10® 

1.068x1022 

2.928x10-1® 

2.379x10® 

3.1 

5.836  xlO^ 

3.026x10® 

1.065x1022 

2.842x10-1® 

2.309x10® 

3.2 

5.852  xlO^ 

3.035x10® 

1.065x1022 

2.851x10-1® 

2.316x10® 

3.3 

5.753  xlO^ 

2.985x10® 

1.065x1022 

2.803x10-1® 

2.277x10® 

3.4 

5.908x10^ 

3.065x10® 

1.065x1022 

2.879x10-1® 

2.339x10® 

4.1 

5.005x10^ 

2.577x10® 

1.023x1022 

2.518x10-1® 

2.046x10® 

4.2 

5.092x10^ 

2.623x10® 

1.023x1022 

2.563x10-1® 

2.082x10® 

4.3 

5.044  xlO^ 

2.599x10® 

1.023x1022 

2.539x10-1® 

2.063x10® 

4.4 

5.040  xlO^ 

2.597x10® 

1.023x1022 

2.538x10-1® 

2.062x10® 

The  results  displayed  in  Tables  13  and  15  show  a  wide  spread  of  results  for  the 
effective  flux.  However,  one  must  remember  that  this  is  not  an  absolute  measure  of 
the  neutron  flux,  but  a  relative  measure  of  the  flux  that  the  Fe  isotope  ean  respond 
to  based  on  their  cross  section. 

A  more  intuitive  measure  of  the  dosimetry  methodology  would  be  to  compare 
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the  measured  activity  at  infinite  irradiation  as  calculated  in  Equation  49  with  the 
expected  activity  at  infinite  irradiation  as  calculated  using  the  Monte  Carlo  derived 
spectrum.  This  is  done  using  Equation  56  by, 


^oo  ^isotope 


(66) 


In  Equation  56,  Nisotope  is  the  number  density  of  parent  isotopes  in  the  dosimeter 
in  cm“^  and  (TlsS.ope  (-^)  ^^e  energy-dependent  cross  section  for  the  reaction  that 

creates  the  radioactive  daughter  product.  Using  the  spectrum  in  Figure  31  and  cross 
sections  for  the  three  reactions  of  interest  as  described  in  Section  4.1,  a  comparison 
can  be  made  and  the  expected  neutron  flux  can  be  verified.  The  three  cross  sections 
are  plotted  in  Figure  45.  Note  that  the  units  of  the  y-axis  in  all  three  plots  are  mb 
or  10“^^  cm^,  meaning  that  the  cross  sections  remain  very  small,  even  at  their  peak 
value. 


£ 

b 


10“  I - ^ ^ ^ — . . ^ - 1  10“ 


Energy  [MeV]  Energy  [MeV]  Energy  [MeV] 


Figure  45.  Cross  sections  for  neutron  capture  by  (left)  (middle)  ^^Fe,  and  (right) 
®®Fe.  Note  that  the  units  of  the  y-axis  are  in  mb. 


In  Figure  46,  the  measured  activities  at  infinite  irradiation  for  all  4  dosimetric 
irradiations  are  plotted  along  with  the  calculated  activity  for  various  levels  of  1-MeV 
equivalent  (Si)  neutron  fiuence.  The  dashed  line  signifies  the  expected  value  for 
neutron  fiuence  for  the  lU  NREF  spectrum,  which  is  4.96  x  10^*^  n/cm^/s.  The  Fe 
foils  agree  well,  within  experimental  error. 

The  flux,  as  determined  using  S  pellets,  does  not  fit  the  other  measurements.  Un- 
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fortunately,  the  beta  detection  system  used  for  this  measurement  included  a  number 
of  uncertainties  for  which  approximations  were  made.  First,  there  was  no  accounting 
for  absorption  in  the  detector  window,  meaning  the  result  would  be  skewed  to  under 
calculate  the  true  flux.  Second,  when  calculating  the  geometric  factor  used  to  correct 
for  the  difference  between  the  standard  calibrating  source  and  the  unknown  S  pellet 
sample,  the  radius  of  the  calibrating  source  disk  was  used.  Most  likely,  however,  the 
beta  source  itself  was  a  very  small  radius  sample  positioned  somewhere  on  the  disk. 
Because  of  these  assumptions,  the  efficiency  was  calculated  to  be  larger  than  it  ac¬ 
tually  was,  and  therefore  the  S  pellet  measurements  again  under  calculated  the  true 
flux.  Using  only  the  Fe  foils,  the  1-MeV  equivalent  (Si)  neutron  flux  is  found  to  be 
4.96  X  10^°  n/cm^/s  for  these  experiments. 

For  each  of  the  four  20-min  dosimetry-focused  irradiations,  a  TLD  crystal  was 
included  in  the  dosimetry  package  placed  next  to  the  BJTs.  This  TLD  was  then  kept 
in  a  light-tight  package  while  it  was  sent  to  the  Naval  Surface  Warfare  Center  (Crane 
Detachment)  for  gamma  flux  determination.  The  average  gamma  dose  was  found 
to  be  (5.3  ±  1.3)  x  10^  rad(Si).  From  Monte  Carlo  calculations,  the  TLD  response 
expected  based  on  the  calculated  gamma  spectrum  is  4.75  x  10^  rad(Si).  One  reason 
the  measured  TLD  response  would  be  higher  than  the  calculated  value  may  be  due 
to  neutron-induced  electron  excitation  that  is  not  taken  into  account  in  the  simulated 
value,  so  having  a  higher  average  gamma  fluence  is  not  unexpected. 

Results  and  Analysis. 

For  the  four  dosimetric  irradiations  of  20  min  each,  B  JT  current  gain  was  measured 
as  a  function  of  time  using  the  Keithley  SMUs.  Using  the  calculated  neutron  flux, 
a  plot  of  the  A(l//3)  is  displayed  in  Figured?.  The  1-MeV  equivalent  (Si)  fluence 
is  calculated  using  ASTM  Standard  E722-09el,  Standard  Practice  for  Characterizing 
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[dps] 


Figure  46.  A  comparison  of  the  measured  and  expected  activity  at  infinite  irradiation 
time  for  the  lU  NREF  neutron  spectrum  provides  good  agreement  for  the  Fe  foils  for 
4.96  X  10^°  n/cm^/s  while  the  S  pellet  does  not  agree. 


110 


Neutron  Fluence  Spectra  in  Terms  of  an  Equivalent  Monoenergetic  Neutron  Fluence 
for  Radiation- Hardness  Testing  of  Electronics.  The  Messenger-Spratt  equation  states 
that  the  performance  degradation  of  a  semiconductor  transistor  when  subjected  to 
displacement  damage  such  as  that  from  neutrons  is  governed  by  a  linear  response  as 
described  by  Equation  12. 


O  [n/cm^] 


Figure  47.  The  A  (1//3)  with  fluence  for  the  4  dosimetric  shots  with  error  corresponding 
to  beam  current  uncertainty  in  the  x-direction  and  equipment  error  in  the  y-direction. 
The  average  slope  was  3.06  (±0.55)  x  10“^®  cnF/n. 

The  response  of  the  BJT  under  irradiation  is  linear  within  experimental  error  and 
uncertainty  with  an  average  slope  of  3.06(±0.55)  x  10“^^  cm^/n.  This  is  greater  than 
the  slope  in  the  ASTM  Standard  E1855  [1],  however  these  devices  have  not  been 
baked  at  80° C  for  the  standard-prescribed  2  hour  time  interval. 

The  second  set  of  irradiations  in  this  experimental  series  investiaged  effects  of  time 
and  variations  to  total  neutron  fluence.  Figure  48  displays  the  A{l//3)  with  fluence 
for  three  of  the  hve  variable  fluence  irradiation  intervals,  which  have  an  average 
slope  of  3.43  x  10“^®  cm^/n.  Unfortunately,  the  beam  current  was  not  constant  from 
irradiation  to  irradiation  so  that  the  total  fluence-to-time  varied,  causing  the  slope 
to  fluctuate.  Additionally,  the  beam  current  varied  in  the  middle  of  the  5th  and 
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Figure  48.  The  A  (1//3)  with  fluence  for  three  of  the  five  variable  fiuence  shots  with  error 
corresponding  to  beam  current  uncertainty  in  the  x-direction  and  equipment  error  in 
the  y-direction.  Shots  5  and  8  were  removed  because  the  beam  current  fiuctuated 
during  the  shot. 

8th  irradiations,  so  they  were  removed  from  Figure  48  and  presented  separately  in 
Figure  50.  Irradiation  6F  was  for  four  min  of  total  irradiation  time,  as  the  proton 
beam  capacitor  bank  tripped  the  safety  breaker,  automatically  shutting  the  beamline 
down.  Irradiations  7F  and  9F  were  for  40  and  60  min,  respectively. 

To  correct  for  the  variation  in  the  beam  proton  flux  from  shot-to-shot,  a  scaling 
factor  was  sought  using  the  beam  monitor  log  displayed  in  Figure  49.  The  beam 
current  for  each  irradiation  was  integrated  over  the  time  interval  of  that  irradiation 
to  get  a  total  proton  current.  Since  the  y-axis  in  Figure  49(a)  is  relative,  absolute 
calculation  of  the  total  current  is  not  possible.  However,  comparison  from  one  ir¬ 
radiation  to  the  next  is  possible  by  calculating  the  ratio.  Shot  IT  was  selected  as 
the  normalizing  irradiation,  and  the  results  of  this  calculation  are  also  displayed  in 
Figure  49(b).  Table  16  displays  the  time  intervals  as  well  as  the  expected  and  actual 
beam  current  ratios,  dehned  as  R  =  In/Iir,  for  each  shot.  It  is  these  values  that  are 
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used  to  scale  the  expected  flux. 


Table  16.  Fitting  coefficients  from  gain  degradation  intervals  with  extracted  flux,  (/>• 


Shot 

tirr  [min] 

Exp.  Ratio 

Act.  Ratio 

Diff 

I/t 

IT 

20 

1 

1 

~ 

2.343 

2T 

20 

1 

0.98 

-2% 

2.304 

3T 

20 

1 

0.95 

-5% 

2.301 

4T 

20 

1 

1.04 

4% 

2.38 

5F 

120 

6 

6.73 

12% 

2.664 

6F 

4 

0.2 

0.17 

-15% 

1.883 

7F 

40 

2 

2.25 

13% 

2.451 

8F 

90 

4.5 

5.15 

14% 

2.687 

9F 

60 

3 

3.19 

6% 

2.514 

From  Figure 49(b),  most  irradiation  intervals  have  a  relatively  constant  beam  cur¬ 
rent  during  the  entire  interval.  However,  in  irradiations  5F  and  8F,  the  proton  source 
jumped  between  two  average  current  levels  during  irradiation  of  the  BJTs.  These  two 
irradiations  are  displayed  in  Figure  50.  They  are  broken  down  into  intervals,  three 
for  5F  and  four  for  8F,  corresponding  to  the  periods  in  which  the  beam  was  in  a  low 
or  high  current  state.  A  linear  regression  was  performed  on  the  three  intervals  in 
irradiation  5F  and  the  four  intervals  in  irradiation  8F;  the  overall  damage  constant 
was  2.94(±0.35)  and  2.810(±  0.39)  xl0“^^  cm^/n,  respectively.  In  chosing  the  change 
in  the  “average  flux”,  the  ratios  determined  from  Figure 49  and  tabulated  in  Table  16 
were  used.  This  analysis  assumes  that  neutron  flux  is  approximately  proportional  to 
beam  current  incident  on  the  beryllium  target.  A  scaling  factor,  SF,  was  determined 
by  comparing  the  average  values  for  each  interval  against  the  normalizing  average 
value  selected  as  the  first  dosimetric  irradiation,  IT. 

The  results  for  the  damage  constant  calculated  using  gain  data  for  each  irradiation 
are  tabulated  in  Table  17.  Overall,  the  average  value  of  K  is  3.17(±0.40)  x  10“^^ 
cm^/n.  This  is  higher  than  the  value  quoted  by  the  ASTM  standard  by  a  factor  of 
approximately  2,  however  these  devices  have  not  been  annealed  in  the  80° C  oven. 
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Figure  49.  (a)  lU  NREF  beam  current  as  a  function  of  time  for  the  experiment,  (b) 
The  proton  beam  current  was  not  constant  throughout  the  experimental  series,  and  an 
example  of  proof  of  this  is  that  the  beam  current  ratio  of  5F,  which  should  be  6.0,  is 
12%  higher. 
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1.5 


Figure  50.  Shots  5F  (top)  and  8F  (bottom)  of  the  A(l//3)  versus  fluence.  In  these 
two  shots,  the  proton  beam  current  jumped  between  low  and  high  values  during  the 
irradiation.  The  blue  dotted  vertical  lines  correspond  to  the  time/fluence  at  which  the 
beam  current  changed  during  the  irradiation. 
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Unfortunately,  these  devices  could  not  be  recovered  to  perform  the  prescribed  bake. 


Table  17.  The  calculated  damage  constant  and  linear  regression  parameter  for  neutron 
irradiation  at  the  lU  NREF  experimental  set  2. 


Shot 

K  (xlO-^5) 
[cm^/n] 

IT 

4.186 

2T 

2.991 

3T 

2.739 

4T 

2.329 

5F 

3.086 

6F 

3.992 

7F 

2.640 

8F 

2.950 

9F 

3.646 

AVG 

3.173 

Comparison  with  traditional  single-pulse  systems. 

In  Figure  75,  a  plot  of  the  change  in  inverse  gain  as  a  function  of  fluence  is  plotted. 
However,  this  time,  the  results  differ  from  the  hrst  experimental  set  at  lU  NREF.  Most 
of  the  data  sets  he  above  the  line  traced  by  the  legacy  data  set,  and  the  variation  is 
most  noticeable  at  higher  fluences. 

4.2  OSURR  Experimental  Series 

After  accomplishing  the  experiments  at  lU  NREF,  the  research  questions  were 
investigated  further  by  utilizing  a  different  neutron  source;  the  Ohio  State  University 
Research  Reactor  facility.  Silicon  BJTs  were  irradiated  with  neutrons  and  the  current 
gain  measured  in  the  same  fashion.  The  gain  was  then  compared  with  experiments 
at  other  facilities. 

The  OSURR  facility  was  used  for  this  experiment  because  the  model  presented  in 
Chapter  III  predicted  a  different  response  to  radiation  than  in  the  lU  NREF.  Three 
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Figure  51.  The  A(l//3)  versus  fluence  for  the  experimental  data  collected  from  the 
second  experimental  series  at  lU  NREF  as  compared  to  data  from  various  single-pulse, 
short-duration  irradiation  facilities  such  as  SPR-III  and  WSMR  FBR. 

irradiation  conditions  were  compared  via  simulation: 

1.  Typical  lU  NREF  conditions  with  beam  parameters  of  20  Hz  and  400  ps 

2.  The  lU  NREF  operating  running  with  beam  frequency  “de-tuned”  to  15  Hz, 
yet  still  having  the  same  time-averaged  flux  as  (1) 

3.  An  OSURR-type  neutron  spectrum  with  the  same  1-MeV  equivalent  neutron 
flux  as  (1)  and  (2). 

In  the  third  simulation,  the  reactor  simulation  included  operation  at  a  percentage  of 
the  maximum  power  to  ensure  the  1-MeV  flux  was  the  same  as  the  time-averaged 
neutron  flux  at  HI  NREF.  It  is  expected  that  using  this  type  of  facility  would  provide 
a  device  response  with  a  damage  constant  lower  than  that  found  at  the  lU  NREF. 

The  most  important  qualitative  feature  of  the  results  of  these  simulations  is  that 
the  NREF  damage  constant  is  different  from  the  OSURR.  At  lower  fluence,  this  is 
difficult  to  differentiate.  However,  as  the  simulation  reaches  higher  total  neutron 
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Figure  52.  Simulation  of  three  different  irradiation  conditions  at  the  lU  NREF  and 
OSURR. 

fluence,  the  three  simulations  begin  to  separate.  The  “detuned”  lU  NREF  conditions 
result  in  the  largest  damage  constant.  This  is  because,  to  have  the  same  1-MeV 
equivalent  (Si)  neutron  flux  with  a  lower  beam  frequency,  the  instantaneous  neutron 
flux  must  be  higher.  With  the  higher  instantaneous  flux,  more  damage  is  created  with 
the  same  amount  of  total  irradiation  time.  However,  higher  flux  increases  the  amount 
of  time  in  which  no  damage  is  occurring  (he.  time  when  the  beam  is  off).  Similarly, 
the  OSURR  has  the  smallest  damage  constant.  All  three  simulations  still  follow  the 
linear  response  with  fluence  as  predicted  by  the  Messenger-Spratt  equation. 

Facility  Description. 

The  OSURR  consists  of  a  light  water,  pool-type  research  reactor.  It  is  licensed 
to  operate  up  to  500  kW  and,  at  maximum  steady-state  power,  the  average  thermal 
neutron  flux  is  approximately  5  x  10^^  n/cm^/s. 

There  are  several  sub-facilities  within  the  reactor  core,  displayed  in  Figure  53, 
that  allow  researchers  to  select  different  neutron  spectrum  shapes  and  flux  rates.  The 
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Central  Irradiation  Facility  (GIF)  has  an  inner  diameter  of  1.3  inches  and  extends 
from  the  top  of  the  pool  down  to  a  central  grid  position  within  the  core.  It  has  a 
maximum  total  flux  of  2.3  x  10^^  n/cm^/s  and  a  maximum  1-MeV  equivalent  (Si)  flux 
of  4.52  X  10^^  n/cm^/s  with  a  hardness  oi  H  =  0.20.  The  Auxiliary  Irradiation  Facility 
(AIF)  is  2.5  inches  in  diameter,  but  the  maximum  fluxes  are  lower  at  0tot  =  9.4  x  10^^ 
n/cm^/s  and  0i_Mey  =  2.44  x  10^^  n/cm^/s,  with  a  hardness  oi  H  =  0.27.  An 
advantage  of  this  is  that  the  gamma  flux  is  also  lower  because  the  additional  water 
between  the  core  and  the  tube  attenuates  the  low  energy  gamma  rays.  The  Peripheral 
Irradiation  Facility  (PIF)  is  also  2.5  inches  in  diameter  and  in  a  different  position, 
with  a  maximum  total  flux  of  5.3  x  10^^  n/cm^/s.  Finally,  the  thermal  column  has 
a  thermalized  neutron  held  with  a  maximum  total  hux  of  6.8  x  10^^  n/cm^/s.  Flux 
values,  except  for  the  calculated  1-MeV  equivalent  (Si)  huxes,  come  from  the  OSURR 
web  site  [17]. 


Figure  53.  Image  looking  down  in  the  OSURR  reactor  pool  at  the  core  during  opera¬ 
tion.  This  shows  the  multiple  tubes  and  beam-ports  available  for  use.  The  7  inch  tube 
used  for  irradiations  was  positioned  at  the  red  X.  [Photo  obtained  from  OSURR  Staff.] 

All  of  these  facilities  have  advantages  and  disadvantages,  however  the  one  that 
was  used  was  the  7  inch  tube.  Its  position  near  the  core  is  displayed  as  the  red  X  in 
Figure  53.  The  7  inch  tube  can  be  moved  within  the  reactor  pool,  but  the  standard 
position  allows  a  maximum  total  flux  of  1.5  x  10^^  n/cm^/s  and  a  1-MeV  equivalent 
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(Si)  flux  of  1.92  X  10^^  n/cm^/s  with  a  hardness  of  0.13  and  a  gamma  dose  of  6.5  x  10"^ 
rad(Si)/hr.  The  maximum  flux  and  hardness  are  lower  than  the  other  tubes  because 
the  7  inch  tube  is  positioned  the  furthest  from  the  center  of  the  core  when  compared 
with  the  other  major  tubes. 

The  greater  distance  from  the  center  of  the  core  leads  to  additional  neutron  ther- 
malization  as  neutrons  travel  through  the  light  water  coolant.  The  spectrum  has  been 
calculated  using  Monte  Carlo  simulations,  and  the  calculations  have  been  verihed  us¬ 
ing  foil  activation  analysis  and  spectrum  unfolding  using  the  Spectrum  Analysis  by 
Neutron  Detectors  (SAND-II)  code  [18].  The  spectrum  in  the  7  inch  tube  is  displayed 
in  Figure  54. 


Energy  [MeV] 


Figure  54.  Neutron  differential  flux  spectrum  at  the  OSURR  7  inch  tube  at  full  power 
where  the  1-MeV  equivalent  (Si)  flux  at  this  power  level  is  1.92  x  10^^  n/cm^/s  and  the 
hardness  is  0.13. 

OSURR  Test  Equipment. 

During  irradiation,  the  two  BJTs  being  measured  were  under  constant  electrical 
bias.  To  accomplish  this,  hve  Keithley  237  SMUs  and  one  Agilent  B2901A  SMU 
were  set  up  at  the  side  of  the  reactor  pool.  Each  BJT  leg  was  connected  to  an  SMU 
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via  approximately  10  to  15  feet  of  BNC  cabling.  The  SMUs  data  was  output  to  a 
laptop,  which  was  running  a  custom-built  LabView  virtual  instrument,  via  General 
Purpose  Interface  Bus  (GPIB)  cables  and  a  National  Instruments  Ethernet-to-GPIB 
converter.  The  setup  is  displayed  in  Figure  55. 


Figure  55.  SMU  set  up  next  to  the  OSURR  reactor  pool. 

At  the  OSURR,  there  is  an  aluminum  frame  with  an  eye  loop  at  the  top  which 
can  be  hooked  to  a  lead  plug  for  the  top  of  the  7  inch  tube.  A  simple  cardboard 
frame  with  BNG  cable  throughputs  was  attached  to  this  aluminum  frame.  A  printed 
circuit  board  (PGB)  was  soldered  directly  to  the  BNG  throughputs,  and  power  sockets 
were  soldered  to  the  PGB.  BJTs  could  be  socketed  to  the  PGB  and  removed  easily 
for  repeatability.  The  dosimetry  packets  were  adhered  to  the  cardboard  frame.  An 
image  of  the  test  rig  set-up  is  displayed  in  Figure  56. 

In  addition  to  the  BJT  being  tested,  a  second  BJT  was  connected  to  the  base 
leg  of  the  device  in  order  to  shunt  the  gamma-induced  photocurrent,  as  displayed  in 
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Figure  56.  The  test  rig  built  for  the  OSURR  experiments  was  simple  and  made  from 
materials  that  would  minimally  activate  so  that  the  experiment  could  be  accomplished 
as  many  times  as  possible  while  requiring  the  minimum  downtime  for  radioactive  decay 
between  shots. 

Figure  34.  The  gamma  flux  at  the  OSURR  is  signiflcant,  and  the  electrons  promoted 
into  the  conduction  band  by  gamma  ionization  in  the  base  can  mask  current  gain 
degradation  during  the  measurement.  Therefore,  this  shunting  procedure  must  be 
done  in  order  to  properly  measure  BJT  gain  degradation. 

OSURR  Test  Procedure. 

Three  irradiations  were  accomplished:  two  at  1%  of  maximum  power,  correspond¬ 
ing  to  1.92  X  10®  n/cm^/s,  and  one  at  23.4%  of  maximum  power,  corresponding  to 
4.96  X  10^®  n/cm^/s.  For  each  of  the  three  irradiations,  two  BJTs  were  placed  in  the 
printed  circuit  board  rig  and  lowered  into  the  core  in  the  7  inch  tube.  There  were  also 
dosimetry  packets  to  validate  the  neutron  flux  for  one  of  the  1%  and  for  the  23.4% 
power  irradiations.  These  packets  consisted  of  2  S  pellets,  2  Fe  foils,  and  2  An  foils. 
For  the  1%  power  dosimetry  packet,  there  were  also  14  TLDs.  The  dosimetry  packet 
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for  the  23.4%  irradiation  did  not  contain  TLDs  because  they  would  have  become  too 
radioactive  and  would  present  a  safety  hazard  when  removed  from  the  test  fixture. 

Initially,  the  core  was  subcritical  and  data  recording  commenced.  The  SMUs  were 
set  to  source  +10  V  on  the  collector,  -220  +A  on  the  emitter,  and  the  base  was  set  to 
SMU  ground.  The  collector  and  base  SMUs  were  set  to  source  voltage  and  measure 
current,  while  the  emitter  SMU  was  set  to  source  current  and  measure  voltage.  This 
conhguration  ensured  that  the  collector-base  junction  was  reverse-biased,  while  the 
emitter-base  junction  was  forward-biased,  causing  the  BJT  to  be  in  forward  active 
mode  at  all  times.  It  is  in  this  mode  that  BJTs  are  designed  to  provide  the  greatest 
current  gain;  this  also  means  that  it  is  most  sensitive  to  changes  in  defect  species 
concentration  during  irradiation. 

During  reactor  start-up,  measurement  continued.  Then,  when  the  reactor  was 
at  the  desired  power  level,  the  twenty  min  timer  began.  The  SMU  output  from  all 
three  legs  was  collected  and,  after  20  min,  the  reactor  was  shut  down.  Output  was 
recorded  until  the  devices  were  removed  from  the  pool.  After  a  period  of  time  to 
allow  short-lived  radioactive  isotopes  to  decay  away,  the  test  rig  was  pulled  from  the 
tube  and  the  devices  and  dosimeters  were  removed  from  their  positions. 

Neutron  Dosimetry. 

The  hrst  step  in  dosimetry  was  to  measure  the  masses  of  the  pellets  and  foils. 
This  data  is  tabulated  in  Table  18.  The  two  neutron  capture  reactions  for  ^^Fe  and 
^®Fe  were  used  for  dosimetry  purposes.  Additionally,  the  gamma  emitter  ^®Fe  was 
detected  with  the  OSURR  data  due  to  neutron  capture  by  ^®Fe.  However,  the  ^®Fe 
isotope  has  an  isotopic  abundance  of  0.2%  in  natural  iron.  Therefore,  its  gamma 
emission  line  is  well-pronounced. 

Equation  49  was  used  to  calculate  the  saturation  activity.  For  the  Fe  dosimeters. 
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Table  18.  Mass  measurement  for  the  OSU  dosimeters. 


Mat’l 

ID# 

Mass  (g) 

+/- 

N 

32S,54Fe,i97Au 

N 

56Fe 

N 

58Fe 

S 

SL40 

0.5923 

0.0001 

1.056x10^2 

SL48 

0.5915 

0.0001 

1.055x102^ 

SLl 

0.6006 

0.0001 

1.071x102^ 

SL7 

0.5954 

0.0001 

1.062x102^ 

Fe 

1 

0.3075 

0.0001 

1.938x102^^ 

3.043x1021 

9.351x101^ 

2 

0.3096 

0.0001 

1.952xl02‘J 

3.063x1021 

9.415xl0i« 

3 

0.3060 

0.0001 

1.929x102^ 

3.028x1021 

9.305xl0i« 

4 

0.3043 

0.0001 

1.918x102^ 

3.011x1021 

9.254xl0i« 

An 

G 

0.5993 

0.0001 

1.533x1021 

H 

0.5194 

0.0001 

1.537x1021 

I 

0.1227 

0.0001 

3.751  X1021J 

J 

0.1208 

0.0001 

3.693  X1021J 

the  counting  time  was  86,400  s,  or  1  full  day.  For  the  S  pellets,  the  counting  time 
was  300  s.  The  results  of  the  measurement  are  tabulated  in  Tables  19, 20,  and  21.  In 
these  tables,  the  acronym  NPA  stands  for  Net  Peak  Area  and  is  the  total  number  of 
counts  under  the  curve  for  a  gamma  spectrum  after  background  adjustment  or  the 
total  number  of  counts  for  the  beta  detection  system  in  the  time  interval. 

Utilizing  this  information,  the  measured  values  of  A^o  can  be  compared  against 
what  is  expected  based  on  Monte  Carlo  particle  simulation.  As  before,  the  activity  at 
inhnite  irradiation  is  based  on  the  Monte  Carlo  calculated  spectrum  and  the  isotopic 
capture  cross  section  using  the  equation 

-^max 

Aoo  =  n,  j  a:^P{E)<P{E)dE,  (57) 

0 

where  is  the  number  of  isotope  x  in  the  sample  of  interest,  is  the  energy- 
dependent  microscopic  neutron  capture  cross  section  for  the  isotope  x,  and  0  {E)  is 
the  energy-dependent  differential  nentron  flux  spectrum. 

Figures  58  and  59  display  plots  of  the  measured  dosimetric  data  as  black  dots  on 
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Figure  57.  Cross  sections  for  neutron  capture  by  ®®Fe;  note  that  the  units  of  the  y-arxis 
are  in  b  this  time. 

the  blue  curves  of  A^o  as  a  function  of  1-MeV  equivalent  (Si)  flux  for  the  hve  isotopes  of 
interest  in  this  analysis.  The  vertical  dotted  line  corresponds  to  the  1-MeV  equivalent 
(Si)  flux  for  the  particular  irradiation,  1.92  x  10®  n/cm^/s  for  the  1%  irradiation  and 
4.96  X  10^®  n/cm^/s  for  the  23.4%  irradiation. 

At  1%  power,  both  dosimetry  packets  agree  well  with  the  expected  1-MeV  equiv¬ 
alent  (Si)  flux  of  1.92  X  10®  n/cm^/s  .  However,  only  dosimetry  packet  has  good 
agreement  for  the  23.4%  power  irradiation.  Dosimetry  packet  ^1  is  consistently  below 
the  expected  value  of  1-MeV  equivalent  (Si)  neutron  flux.  For  the  23.4%  irradiation, 
the  dosimetry  packets  were  positioned  differently  than  for  the  1%  irradiation.  While 
the  packets  were  taped  to  cardboard  next  to  the  BJTs  in  the  1%  irradiation,  the 
packets  were  taped  to  the  aluminum  rig  for  the  23.4%  irradiation.  The  rig  can  rotate 
freely,  and  this  is  most  likely  the  cause  of  the  discrepancy.  If  the  rig  rotated  such  that 
hrst  dosimetry  packet  were  further  away  from  the  core,  the  flux  would  consequently 
be  lower  due  to  1/r^  streaming.  With  this  understanding  and  recognizing  that  the 
BJTs  were  positioned  between  the  two  dosimetry  packets  for  the  23.4%  irradiation, 
the  predicted  values  were  accepted  for  these  two  irradiations.  The  values  of  1-MeV 
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Figure  58.  Comparison  of  measured  dosimetric  data  (black  dots)  for  packet  1  (left)  and 
packet  2  (right)  with  curves  of  saturation  activity  as  a  function  of  1-MeV  equivalent 
flux  with  the  expected  value  of  1.92  x  10®  n/cm®/s  for  1%  power  in  the  7  inch  tube 
at  OSURR  shows  that  the  pellet/foil  set  agrees  well  with  the  Monte  Carlo  simulation 
results  and  OSURR  reactor  estimates,  within  experimental  error. 


Figure  59.  Comparison  of  measured  dosimetric  data  (black  dots)  for  packet  1  (left)  and 
packet  2  (right)  with  curves  of  saturation  activity  as  a  function  of  1-MeV  equivalent 
flux  with  the  expected  value  of  4.96  x  10^®  n/cm®/s  for  23.4%  power  in  the  7  inch  tube 
at  OSURR  shows  that  the  pellet/foil  set  agrees  well  with  the  Monte  Carlo  simulation 
results  and  OSURR  reactor  estimates,  within  experimental  error. 
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Table  19.  Iron  foil  dosimetry  analysis  for  the  OSU  dosimeters. 


ID# 

Rxn. 

Prod. 

Epeak 

[keV] 

^wait 

Det. 

e 

(C-B) 

€ 

Aoo 

[dps] 

Aoo 

Unc.  [%] 

1 

Mn-54 

835 

4371 

2.00x10-3 

3.35x10^ 

3.63x10^ 

52.32% 

Mn-56 

847 

2.28x10-3 

6.79x10^ 

1.05x10^ 

5.85% 

Fe-59 

1099.251 

1.40x10-3 

1.38x10^^ 

1.07x103 

8.66% 

Fe-59 

1291.596 

1.25x10-3 

1.34x10^ 

1.04x103 

9.62% 

Mn-56 

1811 

1.00x10-3 

7.06x10^^ 

8.63x103 

7.37% 

2 

Mn-54 

835 

4350 

2.00x10-3 

5.70x10^ 

6.17x10^ 

22.13% 

Mn-56 

847 

2.28x10-3 

4.66x10^^ 

7.19x103 

6.02% 

Fe-59 

1099.251 

1.40x10-3 

9.00x103 

6.98x10^ 

9.36% 

Fe-59 

1291.596 

1.25x10-3 

1.06x10^ 

8.21x10^ 

9.78% 

Mn-56 

1811 

1.00x10-3 

4.96x10^ 

6.05x103 

7.82% 

3 

Mn-54 

835 

4371 

2.00x10-3 

3.35x10^ 

3.63x10^ 

52.32% 

Mn-56 

847 

2.28x10-3 

6.79x10^ 

1.05x10^ 

5.85% 

Fe-59 

1099.251 

1.40x10-3 

1.38x103 

1.07x103 

8.66% 

Fe-59 

1291.596 

1.25x10-3 

1.34x103 

1.04x103 

9.62% 

Mn-56 

1811 

1.00x10-3 

7.06x103 

8.63x103 

7.37% 

4 

Mn-54 

835 

4350 

2.00x10-3 

5.70x10^ 

6.17x10^ 

22.13% 

Mn-56 

847 

2.28x10-3 

4.66x103 

7.19x103 

6.02% 

Fe-59 

1099.251 

1.40x10-3 

9.00x103 

6.98x10^ 

9.36% 

Fe-59 

1291.596 

1.25x10-3 

1.06x103 

8.21x10^ 

9.78% 

Mn-56 

1811 

1.00x10-3 

4.96x103 

6.05x103 

7.82% 

equivalent  (Si)  flux  were  thus  used  for  the  analysis  to  be  completed  in  the  next  section. 


BJT  Gain  Degradation  Experiment. 

For  this  experiment,  there  were  two  methods  of  BJT  current  gain  data  collection. 
For  one  of  the  1%  reactor  power  irradiations,  gain  was  recorded  as  was  done  in  the 
lU  NREF  experiment  set,  i.e.  +10  V  of  DC  voltage  on  the  collector,  -220  /iA  of  DC 
current  on  the  emitter,  the  base  grounded  to  SMU  ground,  and  the  base  and  collector 
currents  measured  to  calculate  gain.  For  the  23.4%  irradiation,  the  emitter  and  base 
legs  were  biased  at  the  same  levels.  However,  the  collector  voltage  was  swept  from 
-1  V  to  20  V  with  40  data  points  between  -1  and  +1  and  20  data  points  between 
+1  and  +20.  The  reason  for  dividing  up  the  range  in  this  manner  is  that  the  BJT 
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Table  20.  Sulfur  pellet  dosimetry  analysis  for  the  OSU  dosimeters. 


ID# 

^wait 

NPA 

(C-B) 

£ 

Aoo 

[dps] 

Aoo 

Unc.  [%] 

S140.1 

68087 

7.83x10^ 

1.03x10® 

1.53x10® 

5.39% 

S140.2 

68405 

7.80x10^ 

1.02x10® 

1.52x10® 

5.39% 

S140.3 

68720 

7.85x10^ 

1.03x10® 

1.53x10® 

5.39% 

S140.4 

69660 

7.78x10^ 

1.02x10® 

1.52x10® 

5.39% 

S140.5 

71185 

7.72x10^ 

1.01x10® 

1.51x10® 

5.40% 

S148.1 

71605 

7.38x10^ 

9.70x10^ 

1.44x10® 

5.40% 

S148.2 

71920 

7.40x10^ 

9.72x102 

1.44x10® 

5.40% 

S148.3 

72952 

7.40x10^ 

9.73x102 

1.45x10® 

5.40% 

S148.4 

74105 

7.40x10^ 

9.74x102 

1.45x10® 

5.40% 

S148.5 

74417 

7.40x10^ 

9.75x102 

1.45x10® 

5.40% 

Sll.l 

72040 

7.89x10^ 

4.48x10® 

2.31x10^ 

5.15% 

S11.2 

72370 

7.88x10^ 

4.48x10® 

2.31x10^ 

5.15% 

S11.3 

72775 

7.87x10^ 

4.47x10® 

2.31x10^ 

5.15% 

S11.4 

73105 

7.88x10^ 

4.48x10® 

2.31x10^ 

5.15% 

S11.5 

73435 

7.87x10^ 

4.47x10® 

2.31x10^ 

5.15% 

S11.6 

73705 

7.88x10^ 

4.48x10® 

2.31x10^ 

5.15% 

S11.7 

74110 

7.86x10^ 

4.46x10® 

2.30x10^ 

5.15% 

S17.1 

74460 

1.39x10^ 

7.92x10® 

4.09x10^ 

5.12% 

S17.2 

74800 

1.39x10^ 

7.92x10® 

4.09x10^ 

5.12% 

S17.3 

75300 

1.39x10® 

7.92x10® 

4.09x10^ 

5.12% 

S17.4 

75640 

1.39x10® 

7.91x10® 

4.08x10^ 

5.12% 

S17.5 

75979 

1.39x10® 

7.91x10® 

4.09x10^ 

5.12% 

S17.6 

76330 

1.39x10® 

7.93x10® 

4.09x10^ 

5.12% 

S17.7 

76670 

1.39x10® 

7.90x10® 

4.08x10^ 

5.12% 

switches  operation  modes,  essentially  turning  on,  between  -1  and  0  V.  Therefore,  in 
order  to  best  measure  the  transistion  from  negative  current  to  positive  current  in  the 
collector,  the  hner  grid  of  points  was  selected.  Then,  from  +1  to  +20  V,  there  is  little 
change  so  a  coarser  grid  was  selected  to  minimize  the  time  required  for  each  sweep. 


1%  Power  /  DC  Measurements. 


Using  the  SMUs,  BJT  current  gain  was  measured  as  a  function  of  time  for  both 
DuTs  before,  during,  and  after  the  20  min  irradiation  at  1%  of  maximum  reactor 
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Table  21.  Gold  foil  dosimetry  analysis  for  the  OSU  dosimeters. 


ID# 

Epea/c 

^wait 

NPA 

{C-B) 

£ 

Aoo 

[dps] 

Aoo 

Unc.  [%] 

AuG 

411.802 

66650 

3.20x10® 

4.18x10^ 

4.77x10^ 

5.18% 

678.8 

66650 

1.68x10^ 

4.17x10^ 

4.75x10^ 

7.50% 

1087.4 

66650 

2.36x10® 

4.12x10^ 

4.70x10^ 

11.52% 

AuH 

411.802 

67565 

3.13x10® 

4.09x10^ 

4.68x10^ 

5.18% 

678.8 

67565 

1.68x10^ 

4.17x10^ 

4.77x10^ 

5.29% 

1087.4 

67565 

1.92x10® 

3.35x10^ 

3.84x10*^ 

12.35% 

Aul 

411.802 

758872 

6.96x10® 

1.08x10^ 

9.66x10^ 

5.12% 

678.8 

758872 

3.60x10® 

1.08x10^ 

9.66x10^ 

6.75% 

1087.4 

758872 

4.42x10^ 

1.08x10^ 

9.65x10^ 

10.01% 

AuJ 

411.802 

758872 

5.15x10® 

7.23x10® 

6.48x10^ 

5.14% 

678.8 

758872 

2.73x10® 

7.22x10® 

6.47x10^ 

6.98% 

1087.4 

758872 

3.58x10^ 

7.22x10® 

6.46x10*^ 

10.44% 

power.  The  raw  data  is  displayed  in  Figure  60.  The  hrst  25  min  show  reactor  start 
up.  Following  this,  the  reactor  is  maintained  at  the  designated  power  for  20  min,  at 
which  time  annealing  is  monitored  when  the  reactor  is  turned  off. 

Figure  61  displays  the  change  in  inverse  gain  for  the  20  min  irradiation  as  a  function 
of  1-MeV  equivalent  (Si)  ffuence  as  calculated  from  ASTM  Standard  E722-09el  [19]. 
In  this  figure,  there  are  two  distinct  curves:  one  with  and  one  without  gamma  displace¬ 
ment  damage  correction.  It  is  possible  that  the  damage  induced  in  a  semiconductor 
device  can  be  inflated  by  sufficiently  high  gamma  doses.  Therefore,  an  additional 
term  representing  damage  due  to  gamma  recoil  was  included  in  the  Messenger-Spratt 
equation  as  K<h  (t)  =  A-d —  where  is  the  damage  factor  and  is  the 

gamma  dose  in  Gy.  The  value  of  the  damage  factor  stated  in  ASTM  standard  E722  is 
1.5  X  10“^  Gy“^.  This  correction  factor  was  not  applied  for  the  lU  NREF  gain  degra¬ 
dation  analysis  because,  contrary  to  thermal  nuclear  reactors,  the  gamma  flux  at  the 
lU  NREF  is  quite  low.  Therefore,  the  correction  factor  was  negligible  in  comparison 
to  experimental  error. 

In  Figure  61,  the  time-averaged  flux  was  assumed  to  be  the  value  determined  by 
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Raw  Data 


Figure  60.  Raw  data  from  irradiation  of  2  BJTs  at  the  OSURR  at  1%  power  for  20 
minutes. 

dosimetry.  Therefore,  the  slope  of  the  line,  which  corresponds  directly  to  the  damage 
factor  K,  is  calculated  to  be  3.77  x  10“^^  cm^/n.  It  is  interesting  to  note  in  Figure  61 
how  little  of  an  effect  that  gammas  have  on  gain  degradation. 

23.4%  Power  /  Sweep  Measurement. 

Again,  BJT  current  gain  was  measured  as  a  function  of  time  for  2  BJTs  before, 
during,  and  after  the  20  min  irradiation  at  1%  of  maximum  reactor  power.  However, 
this  time,  the  collector  voltage  was  swept  from  -1  V  to  +20  V  once  approximately 
every  1.5  s. 

Figure  62  displays  a  plot  of  all  the  sweeps  from  the  23.4%  power  sweep  measure¬ 
ment.  The  important  features  of  this  graph  are  the  turn-on  region  around  -0.5  V 
and  the  instability  in  the  sweep  at  voltages  greater  than  5  V.  The  turn-on  region 
is  important  because  several  features  about  it  will  be  used  in  later  analysis.  Those 
features  are  voltage  range  from  10  to  90%  of  operating  voltage  and  voltage  at  zero 
current.  As  more  defects  are  injected  into  the  device,  these  parameters  should  get 
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Al/p  (xlO-^) 


Figure  61.  The  A  (1//3)  with  fluence  for  the  20  minute  irradiation;  gain  change  has  been 
corrected  for  gamma  dose  using  the  ASTM  standard  value  of  gamma  displacement 
damage  factor. 
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larger  as  it  will  require  a  greater  potential  to  flip  the  collector-base  region  to  reverse 
bias  and  into  an  “on”  state. 
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Figure  62.  The  collector  current  as  a  function  of  collector  voltage  for  the  26%  power 
shot. 


From  these  sweeps,  the  device  gain  can  be  sampled  at  a  voltage  of  interest.  While 
-|-10  V  was  the  voltage  sourced  to  the  collector  in  the  lU  NREF  experiments,  this 
voltage  could  not  be  used  in  this  analysis.  There  was  an  unresolved  issue  with  auto¬ 
ranging  in  the  SMU  that  was  not  noticed  until  post-experiment  analysis.  This  caused 
destabilization  of  the  SMU  measurement  capability  at  high  voltages  and  caused  the 
measurement  of  collector  current  at  voltages  greater  than  -|-5  V  to  be  unsuitable  for 
analysis.  Therefore,  +0.9  V  was  selected  for  this  analysis,  as  the  collector  current 
changes  very  little  before  reaching  instability  at  +5  V.  Figure  63  is  a  plot  of  the 
Messenger-Spratt  response  for  the  higher  reactor  power  irradiation.  In  this  hgure, 
the  two  irradiated  BJTs  have  differing  responses.  It  is  believed  that  BJT2  shorted  to 
the  power  socket  connections.  Additionally,  BJTl  responded  in  a  non-linear  fashion. 
However,  this  is  likely  due  to  the  reactor  operator  performing  adjustments  in  order  to 
stabilize  the  reactor.  Therefore,  time  is  not  a  linear  function  of  fluence,  which  affects 
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the  A  (1//3)  plot. 


Figure  63.  The  Messenger-Spratt  response  during  irradiation  by  4.96  X  cm? / s 

1-MeV  equivalent  (Si)  neutrons. 

Neutron  dosimetry  proved  that  the  average  1-MeV  equivalent  (Si)  flux  was  4.96  x 
10^°  n/cm^/s  for  this  irradiation;  and  if  the  final  value  of  A  (1//?)  is  0.1254,  then  the 
damage  constant  is  2.11  x  10“^^  cm^/n.  However,  this  value  has  not  been  gamma 
displacement  damage  corrected.  The  gamma  dose  was  217  Gy  leading  to  a  correction 
factor  of  3.26  x  10“^,  changing  the  damage  factor  to  2.05  x  10“^®  cm^/n.  This  value 
of  the  damage  constant  is  much  more  consistent  with  that  found  in  literature,  and  it 
is  different  from  that  found  during  the  lU  NREF  experiments. 

ASTM  Bake. 

The  devices  from  OSURR  versions  1  and  2  were  placed  in  an  oven  and  heated  to 
80° C  for  2  hours  as  prescribed  by  the  ASTM  E1855  standard.  The  device  gain  was 
measured  for  three  devices  after  this  period  of  temperature-assisted  annealing:  (1) 
OSURR  version  1  at  low  power,  (2)  OSURR  version  2  at  low  power,  and  (3)  OSURR 
version  2  at  high  power.  Then,  the  annealing  factor  was  calculated  for  these  three 
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devices.  The  results  are  displayed  in  Figure  64.  The  slope  of  each  line  was  determined 
using  linear  regression,  and  the  average  value  was  — 0.026(±5.97%). 
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Figure  64.  The  annealing  factor  of  the  BJTs  irradiated  at  the  OSURR  as  a  function  of 
time  after  irradiation  is  completed. 


The  A  (1//9)  was  calculated  after  baking  and  the  results  are  plotted  in  Figure  65. 
The  damage  constant,  K,  is  calculated  using  the  initial  gain  and  hnal  gain  in  Equa¬ 
tion  58. 

A=(A1//3)<F  (58) 

The  result  of  applying  Equation  58  is  that  K  =  2.36  x  10“^®  cm^/n  for  the  OSURR 
Version  1  BJT  at  low  power,  K  =  1.72  x  10“^®  cm^/n  for  the  OSURR  Version  2  BJT  at 
low  power,  and  K  =  1.60  x  10“^®  cm^/n  for  the  OSURR  Version  2  BJT  at  high  power. 
It  is  interesting  to  note  that  the  two  Version  2  BJTs  have  a  lower  damage  constant 
than  the  Version  1  BJT,  but  this  is  also  manifested  in  the  in-situ  measurements.  As 
already  mentioned,  this  measurement  suffers  from  not  having  enough  data  points. 
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Time  [s] 

Figure  65.  The  A(l//3)  as  a  function  of  the  point  in  time  where  the  measurement  was 
taken;  note  that  the  a;-axis  is  an  arbitrary  scale. 

Comparison  with  traditional  single-pulse  systems. 

Now  that  data  has  been  collected  from  both  the  lU  NREF  in  pulsed  mode  and 
the  OSURR  at  constant  power,  a  comparison  can  be  made  between  these  two  types 
of  facilities  and  the  legacy,  single-pulse  systems,  as  shown  in  Figure  66. 

The  1//9  data  collected  at  the  OSURR  falls  below  the  data  from  lU  NRFF  at  the 
end  of  irradiation.  Since  there  is  no  accounting  for  real-time  reactor  flux  changes  due 
to  operations,  it  is  best  to  ignore  the  curvature  and  focus  on  the  end  value,  which 
falls  below  the  lU  NRFF  measurements  but  consistent  with  legacy  data.  However, 
as  happened  before  with  the  lU  NRFF  data,  this  is  only  one  set  of  BJTs,  and  only 
one  of  the  two  irradiated  BJTs  provided  an  interpretable  result.  More  experimenting 
is  needed  to  truly  verify  the  damage  constant  as  found  in  Figure  66. 
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Figure  66.  The  devices  respond  differently  for  different  neutron  spectra;  however,  it 
is  yet  to  be  proven  whether  the  ASTM  elevated  temperature  annealing  interval  will 
have  a  significant  enough  impact  to  bring  all  devices  back  to  the  literature  value  for 
the  damage  constant. 
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V.  Comparison  of  Model  to  Experiment 


The  experimental  evidence  was  compared  to  the  mathematical  model  developed  in 
Chapter  III.  Initially,  a  sensitivity  analysis  and  parametric  study  were  accomplished 
to  determine  the  set  of  parameters  required  to  £t  the  experimental  results.  In  this 
chapter,  these  two  studies  are  described,  the  results  are  analyzed,  and  the  resultant 
parameter  set  used  for  comparison  against  the  experimental  data  collected  from  the 
lU  NREF  and  OSURR. 

5.1  Sensitivity  Analysis 

Since  it  was  determined  in  Section  3.2  that  the  values  for  model  parameters  have 
wide  variation  in  literature,  a  sensitivity  analysis  was  warranted  to  determine  which 
were  most  important.  This  sensitivity  study  was  accomplished  when  the  model  con¬ 
tained  six  equations,  so  the  model  used  was  as  shown  in  Equation  Set  30.  In  Equa¬ 
tion  30,  the  six  species  are  vacancies  (V),  divacancies  (VV),  interstitials  (I),  vacancy- 
impurity  complexes  (VX),  divacancy- impurity  complexes  (VVX),  and  impurities  (X). 

The  system  of  six  coupled  differential  equations  has  a  total  of  17  parameters  that 
were  analyzed  for  sensitivity.  These  seventeen  parameters  were  selected  because  of 
their  use  in  the  determination  of  the  rate  coefficients  related  to  Equation  30.  They 
are  tabulated  in  Table  22.  In  this  table,  there  is  the  expectation,  or  most  common, 
value  of  each  parameter  derived  from  the  literature.  Likewise,  the  upper  and  lower 
boundaries  of  the  sensitivity  analysis  were  selected  to  encompass  the  range  of  values 
found  in  literature. 

If  ten  values  were  selected  from  within  each  parameter’s  range,  the  problem  would 
require  10^^  permutations  of  the  parameter  subsets  to  be  solved,  stored,  and  analyzed. 
This  is  not  computationally  feasible  from  a  resource  requirement  or  an  analysis 
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Table  22.  Parameters  (P),  Expectation  values  (<P>),  and  Ranges  (R)  for  17  parameters 
varied  in  sensitivity  study.  Units  of  D  are  cm^/s,  units  of  E  are  eV,  units  of  v  are  s“^ 
and  units  of  Radius  are  A. 


±_ 

P 

<p> 

R 

# 

P 

<p> 

R 

1 

3  X  10-1° 

3  X  10-13  -  3  X  10-° 

10 

^VVX 

1013 

1012  _  1014 

2 

2  X  10-11^ 

2  X  10-1°  -  2  X  10-1^ 

11 

1 

0.1-5 

3 

DvV 

1  X  10-11^ 

1  X  10-1°  -  1  X  10-1^ 

12 

1 

0.1-5 

4 

E 

^VV 

1.6 

1.0  -2.5 

13 

Rv,x 

1 

0.1-5 

5 

E 

^vx 

1.7 

1.0  -2.5 

14 

Rv,vx 

1 

0.1-5 

6 

E 

■^v—vx 

2.24 

o 

1 

15 

Ri.,vx 

1 

0.1-5 

7 

E 

■‘-^vv—x 

2.98 

2.5  -3.5 

16 

Ri.,vvx 

1 

0.1-5 

8 

^vv 

10i° 

1012  _  1014 

17 

Rvv.,x 

1 

0.1-5 

9 

^vx 

1013 

1012  _  1014 

perspective.  Therefore,  an  alternative  method  of  accomplishing  the  sensitivity  analysis 
was  investigated. 

In  1979,  M.D.  McKay  [1]  developed  a  theory,  later  expanded  npon  by  R.L.  Iman 
and  others  [2],  that  nsed  Latin  Squares  to  select  input  parameters  to  best  cover  the 
tested  space.  A  Latin  Square  is  an  n  x  n  array  hlled  with  n  symbols,  each  occurring 
1  time  in  each  row  and  each  column.  A  3  x  3  example  is  shown  in  Figure  67.  A  Latin 
hypercube  is  the  generalization  of  this  concept  to  N  dimensions.  McKay’s  method  is 
often  called  Latin  Hypercube  Sampling  (LHS). 

To  determine  the  simulation  space,  each  variable  range  is  broken  up  into  m  equally 
probable  bins.  Points  are  then  randomly  selected  from  within  those  bins  and  placed 
in  the  simulation  space  to  satisfy  the  Latin  hypercube  requirement  (occurring  only  1 
time  in  each  row,  column,  page,  sheet,  etc.). 

The  mathematical  software  suite  Matlab  [3]  has  a  built-in  function  called  Ihsde- 
sign(n,m)  that  takes  as  arguments  the  number  of  parameters,  n,  and  the  number  of 
bins,  m.  The  output  of  this  function  is  an  n  x  m  matrix  of  numbers  between  0  and 
1  that  satisfy  the  Latin  hypercube  requirement.  By  applying  this  matrix  with  the 
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Figure  67.  This  example  of  a  Latin  Square  shows  the  three  values  A,  B,  and  C  — 
positioned  such  that  each  value  is  used  exactly  once  in  every  row  and  column. 

lower  and  upper  bounds  of  the  parameters,  the  range  can  be  sampled  efficiently.  The 
result  is  a  simulation  space  that  covers  all  ranges  without  requiring  10^^  parameter 
sets  to  sample  the  entire  range. 

For  the  first  sensitivity  study,  the  generation  rate  of  Frenkel  Pairs  (FP)  was  con¬ 
stant  for  400/is.  The  instantaneous  lU  NREF  neutron  energy  flux  was  folded  into  the 
displacement  damage  kerma  for  silicon  as  found  in  the  ASTM  Standard  E722  [4]  and 
is  plotted  in  Figure  17.  The  total  damage  was  integrated  over  the  spectrum  (plotted  in 
Figure  68),  and  the  result  was  divided  by  the  amount  of  energy  needed  for  the  libera¬ 
tion  of  a  silicon  atom  from  the  lattice  and  thus  formation  of  vacancy/interstitial  pair. 
This  was  used  as  the  generation  rate  for  interstitials  and  vacancies,  Gi  =  Gy,  and  as 
an  input  in  the  model.  The  initial  concentration  conditions  {Gy(t  =  0),  Gyy(t  =  0), 
etc.)  and  other  known  parameters,  such  as  the  temperature  and  the  irradiation  time, 
are  input  parameters  to  the  model. 

The  result  of  this  analysis  was  then  plotted,  and  an  example  of  the  results  is 
displayed  in  Figure  69.  In  this  figure,  the  vacancy  concentration  is  on  the  2:-axis  while 
the  X-  and  y-axes  are  time  and  silicon  vacancy  diffusivity,  Dy.  There  is  a  change  in 
the  vacancy  concentration  over  the  range,  which  shows  that  the  input  affected  the 
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Figure  68.  The  time-averaged  differential  neutron  flux  at  the  lU  NREF  that  was  used 
for  sensitivity  analysis. 

output.  The  same  analytical  techniques  were  applied  to  the  other  parameters.  Of 
the  initial  17  parameters  selected  for  sensitivity  analysis,  the  vacancy  diffusivity, 
presented  any  sensitivity  over  its  range.  Therefore,  when  transitioning  into  using  the 
newer  model  with  11  instead  of  6  species  and  subsequently  performing  the  parameter 
search,  only  D.u  was  added  to  the  list  of  parameters  for  this  search. 

5.2  Parameter  Search 

Following  the  sensitivity  analysis,  the  model  was  improved  from  6  species  to  11 
species,  to  include  multiple  defect  clusters.  Additionally,  the  ability  to  calculate 
device  gain  based  on  semiconductor  material  parameters  was  added.  At  this  point, 
the  model  parameters  were  analyzed  again.  This  time,  however,  the  results  predicted 
by  the  model  were  compared  against  experimental  data  collected  during  the  hrst 
experiment  at  lU  NREF.  In  this  way,  the  LHS  methodology  was  used  as  a  searching 
mechanism  for  the  correct  parameters  rather  than  an  analysis  of  the  sensitivity  of 
those  parameters  to  variation. 
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Figure  69.  The  vacancy  concentration  over  a  400  /is  interval  and  for  a  range  of  vacancy 
diffusivity  values  shows  little  change  across  the  entire  range. 

Parameters  were  selected  for  this  analysis  based  on  their  sensitivity,  on  a  lack  of 
values  from  literature,  or  because  they  were  added  after  the  sensitivity  analysis  from 
Section  5.1  was  completed.  To  this  end.  Table  23  displays  the  parameters  and  initial 
values  selected.  In  this  table.  E/p  is  the  energy  required  to  create  a  Frenkel  pair; 
is  the  diffusivity  of  vacancies;  C^o  and  C^o  are  the  initial  concentration  of  vacancies 
and  impurities,  respectively;  Wj,  and  We  are  the  base  and  emitter  widths,  respec¬ 
tively;  and  Nde,  ^de,  and  are  the  donor  concentration  in  the  emitter,  acceptor 
concentration  in  the  base,  and  donor  concentration  in  the  emitter,  respectively. 

Initially,  wide  ranges  of  values  were  selected  for  the  parameters,  which  would  be 
rehned  during  each  iteration  of  the  LHS.  The  LHS  methodology  for  this  analysis 
started  by  creating  the  simulation  space  from  the  9  parameters  and  the  number  of 
bins  for  each  range,  n.  Then,  each  of  the  n  number  of  sets  of  9  parameters  was 
solved  for  the  species’  concentrations  as  a  function  of  time.  These  concentrations 
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Table  23.  First  guess  for  upper  and  lower  bounds  for  parameter  search  using  LHS 
methodology. 


Parameters 

Lim 

Lower 

lits 

Upper 

1.00x10-5 

1.00x10-3 

D„ 

1.00x10-15 

1.00x10-5 

C-yo 

1.00x105 

1.00x1015 

Cxo 

1.00x105 

1.00x1015 

W, 

1.00x10-5 

1.00x10-3 

We 

1.00x10-5 

1.00x10-3 

^DE 

l.OOxlOi^ 

l.OOxlO^o 

Nab 

1.00x1013 

1.00x1015 

Ndb 

1.00x1013 

1.00x1015 

were  converted  to  device  gain  for  comparison  against  the  experimental  data.  Finally, 
a  regression  analysis  of  the  predictive  results  from  the  mathematical  model  and  the 
experimental  results  from  the  hrst  20-minute  dosimetric  measurement  at  lU  NREF 
(see  Section  4.1)  provided  a  measure  of  the  goodness  of  those  . 

Since  the  range  of  the  initial  parameters  were  quite  large,  a  rehnement  procedure 
was  required  to  truncate  the  range  over  which  the  search  proceeded  for  each  subse¬ 
quent  analysis  space.  After  the  initial  simulation  space  was  completed,  a  histogram 
was  generated  for  each  parameter  that  collected  the  number  of  parameters  within 
the  bin  that  met  a  threshold  value  for  the  regression  coefficient,  R^.  Then,  from  this 
histogram,  the  slope  of  the  leftmost  4  and  rightmost  4  bins  was  calculated  and  the 
y  =  0  intercept  was  calculated  for  these  two  lines.  These  two  x  values  became  the  new 
lower  and  upper  bounds,  respectively.  The  results  of  this  procedure  applied  to  the 
initial  LHS  analysis  are  tabulated  in  Table  24.  In  this  table,  the  dark  delineate  when 
a  bound  was  smaller  after  rehnement,  light  cells  delineate  when  a  bound  grew  after 
rehnement,  and  mid  cells  mean  that  there  was  no  change.  In  this  way,  the  parameter 
ranges  grow,  shrink,  or  stay  the  same  through  multiple  rehnements. 
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Table  24.  Change  to  upper  and  lower  bounds  for  the  parameter  search  after  first  cut 
using  LHS  methodology;  colors  are  only  to  guide  the  eye  to  the  parameter  bounds  that 
became  smaller  (red),  larger  (green),  or  stayed  the  same  (yellow). 


Parameters 

LB 

ALB 

UB 

AUB 

1.00x10“® 

-9.00x10“® 

1.00x10“^ 

-9.00x10“® 

1.00x10“^® 

-9.00x10“^® 

1.00x10“®® 

-9.90x10“®® 

C-yo 

1.00x10^ 

1.00x10^ 

1.00xl0®2 

-1.00x10®® 

P 

^XO 

1.00x10® 

0 

1.00x10®® 

9.90x10®^ 

Wb 

1.00x10“® 

-9.90x10“^ 

4.80x10“® 

-9.52x10“® 

We 

1.00x10“® 

-9.90x10“^ 

1.00x10“® 

-9.00x10“® 

^DE 

1.00x10^9 

9.90x10^® 

1.00x10^® 

9.00x102® 

^AB 

1.00x10^® 

9.90x10^^ 

1.00x10®^ 

9.00x10®® 

N db 

1.00x10^^ 

9.00x10^® 

1.00x10®® 

0 

After  the  second  analysis  was  completed,  it  was  obvious  that  the  initial  vacancy 
concentration  (C^o),  emitter  width  {We),  and  donor  concentrations  in  the  emitter  and 
base  {Nde  and  Neb,  respectively)  changed  very  little  over  the  two  rehnement  steps. 
Therefore,  they  were  set  to  a  parameter  within  their  range  and  no  longer  were  a  part 
of  the  search  procedure  in  order  to  reduce  sampling  time.  The  other  hve  parameters 
were  again  shifted  by  the  values  tabulated  in  Table  25. 


Table  25.  Change  to  upper  and  lower  bounds  for  the  parameter  search  after  first  cut 
using  LHS  methodology;  colors  are  only  to  guide  the  eye  to  the  parameter  bounds  that 


became  smaller  (red),  larger  (green 


or  stayed  the  same  (yellow). 


Parameters 

LB 

ALB 

UB 

AUB 

^/p 

1.00x10“^ 

-9.00x10“^ 

1.00x10“® 

-9.00x10“® 

1.00x10“®^ 

-9.00x10“®^ 

1.00xl0“®2 

-9.00x10“®2 

P 

^VO 

1.00x10® 

-1.00x10® 

1.00x10®° 

-9.90x10®® 

Wb 

1.00x10“® 

0 

1.00x10“® 

-4.70x10“® 

Nab 

1.00x10®® 

-9.00x10®® 

1.00x10®® 

-9.00x10®® 

Again,  the  results  of  the  simulation  allowed  for  three  more  parameters  to  be  set: 
initial  impurity  concentration  {Cxo),  base  width  {Wb),  and  acceptor  concentration 
in  the  base  {Nab)-  Therefore,  only  the  Frenkel  pair  creation  energy  {Efp)  and  the 
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vacancy  diffusivity  {D^)  remained. 

In  Figures  70,  71,  72,  {Efp)  and  {Dy)  were  used  and  their  ranges  modified  based  on 
the  results  of  the  >  0.9  tolerance.  Figure  70  plots  the  value  of  the  regression  coef¬ 
ficient  for  those  parameter  sets  that  met  the  0.9  threshold  value.  Figure  71  presents 
the  values  in  Figure  70  as  histograms.  The  top  two  charts  of  Figure  71  shows  that 
significantly  more  values  for  Ejp  meet  the  cutting  criterion  in  the  low  end,  so  the 
range  for  both  parameters  was  modified  in  the  second  analysis  step  (shown  in  the 
middle  two  charts).  This  time,  however,  the  range  for  D^,  was  shifted  too  far.  In  the 
final  two  plots  of  Figure  71,  the  two  parameter  ranges  are  more  stable.  While  they 
could  be  refined  further,  it  was  sufficient  for  analysis  of  the  data. 

Finally,  Figure  72  displays  a  three  dimensional  histogram  of  the  values  of  Ejp  and 
D„  that  met  the  tolerance  criterion.  The  relationship  is  linear,  which  is  expected 
since,  as  E/p  increases, that  fewer  Frenkel  pairs  would  be  produced.  Thus,  it  is  shown 
that  the  diffusivity  of  vacancies,  i.e.  their  ability  to  move  around  and  recombine  or 
annihilate  with  other  species  or  themselves,  must  decrease. 

In  conclusion,  the  parameter  search  resulted  in  the  selection  of  values  for  the  nine 
parameters  selected  for  analysis.  Only  two,  the  vacancy  diffusivity  and  the  energy 
required  to  create  a  Frenkel  pair  E/p  were  found  to  be  sensitive  across  the  range.  This 
shows  that  the  generation  and  diffusion  of  vacancies  dominates  the  prediction  of  gain 
in  the  model. 

5.3  Model/Data  Comparison 

With  reasonable  values  for  unknown  parameters,  a  comparison  to  experimental 
data  was  accomplished.  Using  the  parameters  as  determined  from  the  Latin  Hyper- 
cube  Sampling,  Efp  =  72eV  and  Dy  =  1.05  x  10“^°  cm^/s,  the  comparison  is  plotted 
in  Figure  73.  The  model  underpredicts  damage  for  the  first  ten  minutes  of  the  ex- 
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Figure  71.  The  regression  coefficient  from  Figure  70  is  binned  in  histogram  form  for 
the  (top  two)  first  try,  (middle  two)  second  try,  and  (bottom  two)  final  try. 
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Figure  72.  Three  dimensional  histogram  of  >  0.9  showing  where  both  Dy  and  Ef. 
meet  the  cut. 


Figure  73.  The  change  in  gain  as  a  function  of  time  for  the  model  as  compared  with 
experiment  for  lU  NREF  Shot  IT  (beam  frequency  of  20  Hz,  pulse  width  of  400  fis, 
and  20  mA  beam  current).  There  is  a  small  deficit  when  compared  with  the  data  until 
about  10  min,  at  which  time  it  follows  very  well  within  experimental  error. 


periment.  However,  from  this  point  on,  the  model  matches  the  experimental  data 
within  error  bounds.  The  slope  changes  for  the  model  data  between  200  and  400  sec¬ 
onds,  which  is  not  in  the  experimental  data.  The  behavior  in  the  model  is  explained 
as  the  point  in  which  the  singular  impurity  concentration  drops  below  the  vacancy 
concentration,  so  the  gain  response  begins  to  be  dominated  by  vacancy  concentration 
changes.  It  is  possible  that  this  behavior  is  being  dominated  by  another  mechanism 
not  currently  implemented  in  the  model.  Further  study  and  improvements  to  the 
model  are  warranted  to  correct  these  differences. 

Figure  74  displays  the  annealing  factor  calculation  for  the  model  as  compared 
with  the  empirical  results  from  lU  NREF  Shot  IT.  The  annealing  factor  at  time,  t, 
is  calculated  in  Equation  59. 

_ p 

AF{t)  =  ^ ^  (59) 

/3oo  /3o 
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Figure  74.  The  annealing  factor  as  a  function  of  time  shows  good  agreement  for  the 
comparison  of  the  model  against  experimental  data. 

In  Equation  59,  /3(t)  is  the  gain  at  time  t,  f^o  is  the  gain  in  the  device  before  irradiation, 
and  /doo  is  the  gain  after  a  very  long  annealing  time.  In  Figure  74,  the  gain  at  long 
annealing  time  has  been  taken  to  be  the  gain  at  the  end  of  the  annealing  period.  To 
follow  the  ASTM  Standard  more  precisely,  /3oo  should  take  the  value  of  device  gain 
after  a  period  of  two  hours  in  an  80°C  oven.  However,  since  these  devices  have  been 
misplaced  and  cannot  be  baked,  an  approximation  was  made. 

In  Figure  74,  the  model  is  being  dominated  by  the  two  parameters  Efp  and  D.^. 
The  model  response,  and  thus  the  device’s  predicted  gain,  is  strongly  dependent  on 
the  neutron  flux.  However  in  the  long  term,  where  the  flux  goes  to  zero,  the  response  is 
flux- independent.  This  shows  that  the  short-lived  species,  such  as  vacancies  with  their 
relatively  large  diffusivity,  are  more  sensitive  during  and  immediately  after  irradiation. 
However,  as  the  simulation  progresses  to  later  times,  these  species  are  almost  non¬ 
existent  and  do  not  affect  the  device  gain. 
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VI.  Summary,  Conclusions,  and  Way  Forward 


With  the  model  developed,  experimental  measurements  taken,  and  a  compari¬ 
son  made  between  the  two,  conclusions  can  be  made  regarding  the  outcome  of  the 
research.  In  this  section,  a  summary  of  the  results  of  this  research  are  presented, 
conclusions  that  can  be  drawn  from  the  results  of  this  research  are  presented,  and  the 
plans  for  future  work  are  suggested. 

6.1  Summary  of  Research 

In  the  first  part  of  this  research,  a  mathematical  model  was  developed  that  de¬ 
scribes  the  interaction  of  defect  species  after  they  are  formed  in  a  time-  and  energy- 
dependent  neutron  environment.  The  mathematical  model  proceeded  through  a  num¬ 
ber  of  steps  in  development,  starting  out  with  Bateman-like  equations  of  the  accu¬ 
mulation  of  defect  species  in  the  same  manner  that  radioactive  atoms  build  up  in 
irradiated  materials.  An  analytical  solution  was  sought;  but  as  the  number  of  de¬ 
fects  that  must  be  tracked  expanded,  the  subsequent  complexity  of  the  system  of 
differential  equations  made  this  impossible. 

At  this  point,  a  change  to  a  chemical  kinetics  approach  modihed  the  system  of 
equations.  Each  type  of  interaction  between  species  was  modeled  as  a  rate-dependent 
chemical  interaction  of  two  species  coming  together  to  form  a  third  or  a  single  species 
breaking  into  two  other  species.  The  rate  of  each  interaction,  designated  a,  was 
dehned  using  parameters  such  as  the  diffusivity  of  the  species  interacting,  the  radius 
of  capture  of  one  species  by  another,  and  a  rate  of  dissociation  of  the  species.  The 
system  of  equations  was  expanded  to  include  9  different  defect  species  and  the  two 
semiconductor  charge  carriers,  electrons  and  holes.  The  nine  species  modeled  were 
vacancies,  divacancies,  interstitials,  vacancy-impurity  complexes,  divacancy-impurity 
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complexes,  impurities,  2-vacancy  clusters,  3-vacancy  clusters,  and  annihilations.  The 
result  of  the  development  was  a  system  of  equations  displayed  in  Equation  Set  46.  In 
each  equation,  C  designated  the  concentration  of  a  particular  species  {v  for  vacancy, 
X  for  impurity,  etc.). 

In  Equation  Set  46  was  developed  the  idea  of  the  generation  rate  of  vacancies/in¬ 
terstitials,  and  G*  respectively.  These  generation  rates  took  into  account  the 
displacement  damage  kerma  when  calculating  the  amount  of  energy  deposited  into 
the  material.  Additionally,  the  gamma  generation  rate,  G.y,  was  set  as  the  number  of 
electron/hole  pairs  created  in  the  material  by  the  gamma  dose. 

The  solution  of  this  system  of  coupled,  non-linear  ordinary  differential  equations 
is  a  system  of  defect  species  concentrations  as  a  function  of  time,  and  thus  fluence  of 
neutrons.  An  example  of  the  resultant  solution  is  plotted  in  Figure  25. 

This  system  of  defect  species  concentrations  was  then  nsed  to  calculate  the  BJT 
current  gain  in  the  device.  The  vacancy,  2-vacancy  clusters,  and  3-vacancy  cluster 
concentrations  contributed  to  the  acceptor  concentration  in  a  region,  since  the  va¬ 
cancy  lies  very  close  to  the  valance  band  within  the  silicon  bandgap  structure.  The 
divacancies  acted  as  thermal  electron/hole  pair  generators,  as  they  sit  around  the 
mid-point  in  the  silicon  bandgap.  The  Messenger- Spratt  behavior  was  investigated 
by  plotting  the  change  in  inverse  current  gain  as  a  fnnction  of  neutron  flnence.  An 
example  is  displayed  in  Figure  26  showing  that,  for  a  10  second  irradiation  at  a  1-MeV 
Equivalent  (Si)  nentron  flux  of  (pi-MeV  =  5  x  10^°  n/cm^/s,  the  hypothetical  silicon 
device  responds  in  a  linear  (thus  Messenger-Spratt  like)  manner. 

Following  the  development  of  the  model,  a  series  of  two  experiment  sets  were 
accomplished  at  the  lU  NREF  nentron  source  and  the  OSURR,  which  is  a  thermal 
nentron  reactor-based  system. 

At  lU,  two  types  of  experiments  were  conducted.  A  Multiple  Sequential  Interval 
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(MSI)  test  looked  at  the  response  of  the  Si  BJT  to  multiple  intervals  of  irradiation  and 
annealing,  all  in  a  row.  Each  interval  of  irradiation  could  be  mapped  to  Messenger- 
Spratt  behavior. 

In  addition,  the  annealing  intervals  were  £t  to  a  power  law  function.  With  no  more 
than  30  s  of  annealing  information  for  this  particular  experiment,  the  power  law  fit 
the  data  extremely  well.  However,  it  was  noted  that  a  power  function  is  not  the  best 
representation  as  the  device  gain  would  go  to  inhnity  at  inhnite  time,  whereas  the 
actual  behavior  of  BJTs  after  irradiation  is  to  asymptotically  approach  a  value,  less 
than  /So,  at  long  time  after  irradiation.  Therefore,  it  was  verified  that  device  response 
has  more  than  one  annealing  mode. 


Figure  75.  The  A(l//3)  versus  fluence  for  the  experimental  data  collected  from  the 
second  experimental  series  at  lU  NREF  as  compared  to  data  from  various  single-pulse, 
short-duration  irradiation  facilities  such  as  SPR-III  and  WSMR  FBR. 

For  the  Single  Long  Interval  (SLI)  test,  the  BJT  was  irradiated  in  the  III  NREF  for 
one  interval  corresponding  to  a  particular  1-MeV  Equivalent  (Si)  interval,  and  current 
gain  was  measured  in-situ  during  and  after  the  irradiation.  Dosimetry,  via  neutron 
activation  analysis,  verihed  the  neutron  flux  to  be  approximately  5  x  10^°  n/cm^/s 
because  of  the  mechanical  jack  that  positioned  the  DuT  closer  to  the  beryllium  target. 
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Figure  75  showed  that  the  change  in  inverse  current  gain  (A1//3)  is  linear  with 
fluence  for  many  different  intervals  of  irradiation.  Not  only  that,  but  the  damage 
constant  during  irradiation  was  greater  than  that  found  in  literature  for  the  same 
devices.  However,  it  was  noted  that  the  literature  values  were  determined  after  the 
damaged  devices  had  been  baked  for  two  hours  in  an  oven  at  80°C. 


Figure  76.  Simulation  of  three  different  irradiation  conditions  at  the  lU  NREF  and 
OSURR. 

Next,  a  comparison  was  made  between  the  sequentially  pulsed  source  at  lU  NREF 
with  a  steady-state  source  at  the  OSURR,  undertaken  because  the  model  predicted  a 
different  gain  degradation  response.  Figure  76  displays  this  difference.  The  OSURR 
provided  a  steady-state  flux  of  neutrons,  in  contrast  to  the  multiple  pulses  at  the  lU 
NREF,  to  reach  the  desired  total  neutron  fluence.  In  addition,  the  neutron  energy 
spectrum  for  the  OSURR  is  much  softer  than  the  spectrum  at  lU  NREF.  The  semi¬ 
conductor  devices  were  measured  at  the  OSURR  at  two  power  levels,  1%  and  23.4% 
of  maximum  reactor  power,  which  corresponded  to  a  neutron  flux  of  2  x  10®  and 
4.94  X  10^®  n/cm^/s  in  the  experimental  column.  The  SLI  test  was  repeated  at  the 
OSURR;  the  BJTs  being  measured  were  irradiated  for  approximately  20  min  and  the 
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gain  measured  in-situ,  the  same  as  at  lU  NREF.  The  result  of  this  analysis  showed 
that  A1//3  remained  linear,  but  with  a  different  damage  constant  than  the  one  de¬ 
termined  during  the  lU  NREF  experiment  set.  At  the  OSURR,  the  damage  constant 
measured  during  irradiation  was  between  1.60  and  2.36  x  10“^®  cm^/n,  which  is  much 
closer  to  that  published  in  the  ASTM  standard  value  of  2  x  10“^^  cm^/n  than  that 
measured  at  the  lU  NREF,  which  was  >  3  x  10“^®  cm^/n. 

One  hnal  comparison  was  made  against  legacy  systems,  and  it  is  plotted  again  in 
Figure  77,  this  time  with  a  logarithmic  scale  on  the  x  and  y  axes.  This  allows  one  to 
observe  the  variations  in  BJT  current  gain  response  at  lower  fluences  that  may  not 
be  immediately  evident  when  plotting  on  a  linear-linear  scale.  The  curvature  in  the 
OSURR  data  may  be  a  result  of  changes  in  neutron  flux  due  to  a  rise  in  temperature 
in  the  core,  causing  the  flux  to  be  elevated.  Once  it  was  noticed  by  the  reactor 
operator  and  adjusted,  the  flux  was  depressed  to  compensate. 


Figure  77.  When  comparing  all  datasets  together,  the  devices  respond  differently 
for  different  neutron  spectra;  however,  it  is  yet  to  be  confirmed  whether  the  ASTM 
elevated  temperature  annealing  interval  will  have  a  significant  enough  impact  to  bring 
all  devices  back  to  the  literature  value  for  the  damage  constant. 


To  complete  this  research  work,  a  comparison  was  made  between  the  empirical  re- 
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suits  and  predictions  from  the  mathematical  model.  However,  the  model  parameters 
were  found  to  vary  widely  throughout  literature,  so  a  sensitivity  analysis  and  param¬ 
eter  search  were  undertaken  to  determine  parameters  that  best  £t  the  experimental 
results.  A  Latin  Hypercube  Sampling  technique  was  developed  to  more  efficiently 
search  through  the  entire  range  of  parameters’  values.  The  sensitivity  analysis,  ac¬ 
complished  when  the  model  only  contained  six  defect  species,  fonnd  that  the  vacancy 
generation  rate  parameter  had  the  greatest  affect  on  results  in  that  current  imple¬ 
mentation. 

However,  the  model  was  modified  and  additional  species  were  added.  Therefore, 
a  new  comparison  was  accomplished  using  an  eleven-species  model.  The  parameters 
with  greatest  sensitivity  were  the  Frenkel-pair  generation  energy  {Efp),  vacancy  diffn- 
sivity  {D^),  vacancy  and  impurity  initial  concentrations  (C^o  and  C^o,  respectively), 
base  and  emitter  widths  (H4  and  We,  respectively),  and  dopant  concentrations  in  the 
base  and  emitter  [N^e,  N ab,  andA’ce). 

A  series  of  range  rednctions  were  made  based  on  the  sensitivity,  and  it  was  fonnd 
that  Efp  and  Dy  were  the  most  sensitive  of  these  parameters.  After  setting  all  of 
these  parameters  to  their  most  likely  valne  based  on  this  sensitivity  analysis  and 
parameter  search,  the  model  was  used  to  predict  the  response  of  a  hypothetical  BJT 
snbjected  to  a  20  min  irradiation  at  the  lU  NREF.  The  nentron  sonrce  was  configured 
to  mimic  typical  operating  parameters  at  lU  NREF,  such  that  the  beam  frequency 
was  20  Hz,  the  nentron  pulse  width  was  400  ps,  and  the  proton  beam  current  was 
20  mA.  The  model  was  shown  to  underpredict  damage  in  the  first  20  min,  but  the 
response  matches  well  for  the  rest  of  the  20  min  irradiation  interval.  However,  the 
annealing  behavior  did  not  match  as  well.  While  it  stayed  within  experimental  error, 
the  model  responded  by  almost  instantly  relaxing  to  its  final  valne  of  gain.  On  the 
other  hand,  experimental  evidence  indicated  annealing  with  longer  time  constants.  It 
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was  surmised  that  this  is  due  to  the  need  for  a  more  complex  model  to  account  for 
the  long-time  annealing  mechanisms. 

6.2  Conclusions 

The  kinetics-based  model,  in  its  current  implementation,  is  able  to  predict  device 
gain  degradation  with  reasonable  accuracy.  However,  the  annealing  behavior  of  the 
device  does  not  match  experimental  results  as  well.  It  is  suspected  that  this  is  due 
to  the  device  behavior  being  dominated  during  irradiation  by  generation,  which  is 
directly  driven  by  the  parameter  Efp.  The  kinetic  interaction  of  species  has  little 
effect  during  irradiation,  and  this  is  observed  in  the  empirical  results,  as  the  temporal 
prohle  has  little  effect  on  the  linearity  of  device  response.  However,  during  annealing, 
defect  generation  ceases  and  only  species’  interactions  drive  the  behavior.  Therefore, 
there  must  be  some  aspect  of  these  interactions  that  is  missing  in  the  model.  This 
may  be  a  missing  defect  species  or  the  interpretation  of  how  these  species  are  used  in 
calculating  gain  may  be  incorrect. 

This  work  is  the  hrst  measurement  into  the  advantages  and  disadvantages  of  using 
an  extended  irradiation  period  to  simulate  total  displacement  damage  in  silicon  semi¬ 
conductor  devices.  It  was  revealed  that  the  devices  continue  to  follow  the  Messenger- 
Spratt  linear  behavior  in  the  change  in  inverse  gain  as  a  fnnction  of  flnence  for  both  a 
seqnentially  pulsed  source  such  as  lU  NREF,  as  well  as  a  steady-state  source  such  as 
OSURR.  Non-linearity  was  observed  in  device  response  for  the  26%  power  irradiation 
at  the  OSURR,  however  deviations  in  the  reactor  flux  leave  some  of  the  results  in 
question.  Therefore,  the  hypothesis  that  the  temporal  prohle  of  the  nentron  sonrce 
will  cause  a  deviation  from  linear  behavior  has  proven  false. 

It  continnes  to  be  extremely  important  to  fully  understand  the  time-averaged  nen¬ 
tron  energy  spectrnm  for  the  facility  using  compntational  particle  physics  techniqnes 
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combined  with  experimental  verification,  especially  if  the  device  is  being  measured 
in-situ.  This  is  due  to  the  effect  that  the  neutron  spectrum  has  on  the  types  and 
quantities  of  defect  species  as  they  are  created  and  interact  in  the  semiconductor  ma¬ 
terial  during  irradiation.  This  research  has  shown  that,  during  irradiation,  the  harder 
spectrum  at  lU  NREF  causes  damage  at  a  faster  rate  than  the  softer  spectrum  at 
the  OSURR  at  the  same  fiuences.  Unfortunately,  the  devices  irradiated  at  lU  NREF 
were  misplaced  and  cannot  be  baked.  However,  the  difference  in  damage  constant  is 
small  enough  that  an  80° C  bake  would  most  likely  anneal  enough  damage  to  bring 
the  damage  constant  back  to  the  expected  value  from  literature,  which  was  the  case 
with  the  devices  at  the  OSURR.  For  this  reason,  the  hypothesis  that  neutron  energy 
spectrum  has  a  profound  effect  on  the  damage  constant  is  not  supported  in  this  work. 
It  is  true  in  the  sense  that  the  spectrum  has  an  effect  in-situ]  however,  the  ASTM 
standard  bake  requirement  accounts  for  this.  Therefore,  the  hypothesis  is  false  in 
regards  to  long-term  effects  in  the  devices. 

This  research  has  provided  better  insight  into  why  and  how  facilities  such  as  lU 
NREF  and  OSURR  can  be  used  for  displacement  damage  studies  of  semiconductor 
devices. 

6.3  Ideas  for  Future  Research 

The  first  goal  for  any  future  research  would  be  to  refine  the  model.  As  noted 
in  Section  6.2,  the  annealing  behavior  of  devices  does  not  match  well  with  empirical 
results.  The  procedure  for  generation  of  defect  species  could  be  modified  to  better 
reflect  the  current  understanding  of  defect  species  formation.  To  do  this,  the  list  of 
defect  species  that  are  tracked  currently  would  need  to  be  modified.  In  the  current 
implementation,  2-  and  3-  vacancy  clusters  are  possible,  but  only  by  the  interaction 
of  vacancies.  These  2-  and  3-vacancy  clusters  could  be  removed  from  the  system  of 
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equations  and  could  be  replaced  by  a  species  called  a  cluster.  Now,  the  generation  of 
Frenkel-pairs,  which  is  driven  by  low  energy  neutrons,  could  be  supplemented  by  the 
generation  of  clusters,  which  is  driven  by  higher  energy  neutrons.  This  would  involve 
a  much  more  complicated  routine  for  determination  of  the  generation  rates  than  is 
currently  implemented  in  the  model.  However,  this  would  more  accurately  describe 
how  defect  species  are  created  in  devices  under  neutron  irradiation.  Figure  78  is  an 
example  of  how  the  spectrum  at  lU  NREF  might  be  broken  up  into  three  regions:  one 
that  creates  Frenkel-pairs,  one  that  creates  clusters,  and  one  that  has  the  probability 
to  create  either. 
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Figure  78.  The  differential  neutron  energy  spectrum  at  lU  NREF  with  an  example 
of  the  way  the  spectrum  may  be  broken  up  into  the  region  that  creates  Frenkel-pairs 
(left),  the  region  that  creates  clusters  (right),  and  the  region  that  has  some  probability 
of  creating  either  (middle). 


Another  idea  for  future  research  is  to  expand  the  SLI  fluence  dataset.  Many  of 
the  datasets  contain  points  that  consist  of  one  measurement  at  a  particular  fluence, 
mainly  because  of  time  and  hnancial  restrictions.  Getting  better  statistics  by  mea¬ 
suring  multiple  devices  at  that  fluence  will  lead  to  decreasing  error  and  rehning  the 
understanding  of  the  device  response.  Additionally,  the  maximum  fluence  reached 
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during  this  research  was  <  4  x  10^^  n/cm^.  It  would  be  useful  to  irradiate  a  device 
to  higher  fluences  to  see  if  there  is  a  catastrophic  threshold,  he.  a  fluence  where  the 
defects  become  so  tightly  packed  that  the  neutrons  begin  to  have  less  effect  on  device 
gain.  The  mathematical  model  in  its  current  iteration  predicts  this  will  occur  but  not 
until  much  higher  fluences  are  achieved. 

Finally,  once  silicon  devices  have  been  investigated  fully,  it  would  be  important  to 
apply  these  concepts  to  novel  semiconductor  materials  such  as  III-V  devices  (GaAs, 
etc.),  carbon  nanotubes,  and  others.  Obviously,  it  would  be  necessary  to  modify 
the  input  parameters  for  the  model,  however  that  may  be  the  only  required  change. 
Unfortunately,  as  was  noted  with  silicon,  many  of  these  parameters  are  not  well-known 
or  are  still  being  determined  via  current  research.  Therefore,  it  may  be  necessary  to  do 
calculations  or  experiments  to  better  quantify  the  unknown  parameters  or  determine 
them  through  a  LHS  technique. 
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